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Critical Flow Velocities for Collapse of 
Reactor Parallel-Plate Fuel Assemblies 


Theoretical formulas are presented for prediction of the fiow velocity at which collapse 


occurs in long parallel-plate assemblies. 


Beyond the critical velocity the pressure- 


unbalance forces, developed as a consequence of a small deflection, exceed the correspond- 


ing elastic restraining forces, and the plates collapse 


Both flat and curved-plate as- 


semblies are considered, and the applicability of the formulas to design of reactor fuel- 
plate assemblies is discussed. 


| COLLAPSE of reactor-fuel-plate mockups was ob- 
served in tests conducted at Oak Ridge ten years ago [1],? and 
this problem has been encountered more recently in the original 
fuelelements or the Engineering Test Reactor, ETR [2]. When 
the flow velocities in the Oak Ridge tests were raised to sufficiently 
high values, flat plates were deformed plastically, and curved 
plates were buckled. The flow channels between plates were re- 
ported to have been ‘“‘almost or even completely closed through 
part or all of the length of the plates.’’ The ETR reactor was 
initially operated below its rated flow velocity to avoid collapse of 
its fuel plates. 

The analysis which follows applies the Bernoulli theorem for 
incompressible flow to determine the pressure differences which 
would be developed across the plates by modification of local 
fluid velocities as a result of plate deflections. Equating these 
pressure differences to those which would be required to produce 
the corresponding deflections gives an expression for the critical 
velocity in terms of plate and fluid parameters. 


Assumptions Used in Analysis 
Unless otherwise stated, the following assumptions were made 
in derivation of the formulas for critical velocities: 


1 The plates are homogeneous, isotropic, elastic, initially flat 

1This paper is based on work sponsored by the U. 8. Atomic 
Energy Commission under Contract No. W-31-109 Eng 52 with the 
General Electric Company. 

2 Numbers in brackets designate References at end of paper. 

Contributed by the Nuclear Engineering Division and presented at 
the Annual Meeting, New York, N. Y., November 30-December 5, 
1958, of THe AMERICAN SocieTY OF MECHANICAL ENGINEERS. 

Nore: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 
the Society. Manuscript received at ASME Headquarters, April 8, 
1959. This paper was not preprinted. 
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SEC.C-C SEC. A-A 
Fig. 1 Outlines of plate assembly 


SEC. B-B 


or uniformly curved, uniform in spacing and dimensions, and free 
of unidentified sources of deformation. 

2 The coolant is incompressible; all channels have the same 
mass flow; at any cross section normal to the longitudinal axis the 
flow within any channel is uniform; and leakage between chan- 
nels is suppressed. 

3 Plates are broad enough in comparison with their thickness 
that shear deformation is negligible, and are long enough in com- 
parison with their breadth that plates can deflect locally without 
significant redistribution of flow among the coolant channels. 

4 Side plates or supports are rigid. 


Derivation of Formulas for Critical Flow Velocity 


Case (a) Flat Plates With Built-In Edges. This case represents a 
configuration which has been used or proposed for many reactors. 
On the order of ten similar plates with length large in comparison 
with width are assembled with uniform spacing and with their 
long edges restrained by attachment to side plates. Fig. 1 shows 
schematically the shape of the original transverse cross section of 
such an assembly and also shows the nature of the distortion 
which is postulated at the lowest critical flow velocity. It should 





Nomenclature 


moment tending to straighten 


= plate thickness 

= width of flat plate or chord of 
curved plate 

gravitational acceleration 

initial flow channel thickness at 
midspan 

longitudinal distance between 
plate-edge supports 

pressure difference tending to 
bend or stretch fuel plate 

pressure difference imposed 
across plates due to irregu- 
larity in configuration or 
flow distribution 
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longitudinal distance from edge 
support 

deflection of plate relative to 
supports 

deflection of curved plate rela- 
tive to midspan tangent 

area of beam cross section per 
unit width of beam 

a function of a defined by 
Equation (41) 


= Young’s modulus of elasticity 


moment of inertia of beam 
cross section per unit width 
of beam 


curved plate per unit width 
of beam 

membrane tensile load in 
curved plate per unit 
width of beam 

membrane compressive load in 
flat plate per unit axial 
length of plate 

critical value of P in absence of 
other loads 

original cross-sectional area of 
coolant channel 


(Continued on next page) 
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Fig. 2 Load-deflection diagram 


be noted that the deflections of adjacent plates are in opposite 
directions. This condition corresponds to a lower critical velocity 
than would be characteristic of the case of a single plate which 
divides a rectangular duct into two channels, with the same plate 
and channel dimensions as the multiple-plate assembly. 

At the critical velocity the modifications of local linear velocity, 
resulting from a deflection of the plates, produce across the plates 
a pressure difference sufficient to maintain the deflection. If the 
plates are deflected over a longitudinal distance several times the 
width of the plates, then near the center of the deflected portion 
of a plate the longitudinal variations in deflection can be ne- 
glected, and the deflection curve of the plate shown in Fig. 2, 
modified from the narrow-beam deflection curve from [3] by use 
of the factor (1 — v?) for a wide beam is 


px 1 — y*)(xz? — 2bz + b?) 


y= = (1) 
’ 24ET 


Integrating (1) and multiplying by two to determine the total 
change in channel cross-sectional area gives 


(1 — v?)pbs 


360K] 


(2) 


AQ = 


Substitution of the moment of inertia of a strip of unit width in 
equation (2) and division by bh gives the per-unit change in chan- 
nel cross section, 

AQ (1 — v?)pb! 


= (3) 
Qo 30Fa*h 


The pressure differential developed across the plate away from 
the upstream and downstream ends of the deformed regions, 
derived from Bernoulli’s theorem, is equal to the difference in 
velocity heads on the two sides of the plate multiplied by the 
weight density of the fluid 
Vo Vo 

AQ AQ 
a 4 


1 
Qo Qo 


which, as AQ/Q, approaches zero, reduces to 


2pV.27AQ 
p= —— 
gQo 

Substitution of (3) in (4) determines the critical velocity, 


l5gEa%h °]'/2 5) 
= — (5 
pb(1 — v?) 


The configura- 


c 


Case (b) Flat Plates With Simply Supported Edges. 
tion considered here is the same as for Case (a) except that the 
long edges of the plates are simply supported. Means are assumed 
to be provided to prevent leakage from one channel to another 
without additional restraint on the plates. 

Use of the expression for the deflection curve of a uniformly 
loaded wide beam with simply supported ends and the procedures 
of Case (a) lead to the following expression for the critical 
velocity 


y 5gEa*h ae 
© | 2pb41 — vy?) 


Fig. 3 Mode of distortion with longitudinally spaced edge supports 


Case (c) Flat Plates With Longitudinally Spaced Supports. Fuel 
plates might be supported only at uniform intervals along their 
long edges. This arrangement would introduce the possibility 
that the lowest critical velocity will correspond to deflections 
which vary along the length instead of across the width of the 
plate. The case of this type which is most easily analyzed is that 
in which the plate deflection varies only in the axial direction, and 
the simple supports only provide local control of spacing. This 
mode of collapse would be approximated if the supports were 
spaced at longitudinal intervals long in comparison with the 
plate width, so that transverse curvatures would be small in 
comparison with longitudinal curvature. Fig. 3 shows a sketch 
of the plates and their deflection curves. In this case all plates 
would tend to deform in a similar manner over their full length 
with adjacent plates deflecting in opposite directions at any 
longitudinal position. 

The velocity at any cross section in an enlarged channel is 





Nomenclature 


change in Qo 

initial mean radius of curva- 
ture of curved plate 

linear velocity of coolant in de- 
formed channels 

linear velocity of coolant in 
undeformed channel 

value of Vo for which collapse 
occurs due to flow-induced 
effects alone 

= collapse velocity with superim- 
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posed transverse compression 
= ratio of V, for curved-plate as- 
sembly to V, for flat-plate 
assembly for the case of fixed 
longitudinal edges 
same as V,, except applicable 
to hinged longitudinal edges 
ratio of V, for curved-plate as- 
sembly with fixed longitu- 
dinal edges to V, for curved- 
plate assembly with hinged 
longitudinal edges. 


shear force at edge of curved 
plate per unit width of beam 

curved-plate are between sup- 
ports 

angles of curved-plate are 
measured from midspan 

membrane stress and strain in 
curved plate in direction of 
its width 

weight density of coolant 


Poisson’s ratio 
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V. 1 Voh - 


h + 2y 
and in a constricted channel is 


Voh 
ee 8 
ee (8) 
The pressure difference across the plates, from Bernoulli’s 
theorem and equations (7) and (8) is 
Vo? 1 s 
pa Hs — —_—— (9) 
29 2y 2 
1— 
h 


In the limit as y/h approaches zero, (9) becomes 


tpVo7y 


10 
gh a 


| 


The equation for the deflection curve of one of the plates is 


KI d‘y 
-—— oe (11 
(1 — v*) drt P ) 


Equating (10) and (11) gives 


kI d‘y tpVo*y 

sit Tae gg eee (12) 

(1 — vy?) dx‘ gh 

The solution of (12) corresponding to the simplest mode of de- 
flection and the lowest critical velocity is 


Tr 
y = Bsin ~ (13) 


where B is an arbitrary constant. 
Substitution of (13) and the expression for J in terms of plate 
thickness in (12) gives the lowest critical velocity 


4qEa3h 1/2 
V, = T*gharh ; (14) 
48pl{1 — v?) 


It is of interest to note that for this case the local differential- 
pressure loading on the beam is proportional to the local deflection 
of the beam, just as is true for free vibration of a beam. Because 
of this similarity, the critical velocities for various deflection 
modes have the same ratios as the corresponding natural fre- 
quencies of beam vibrations. Using natural frequency formulas 
from Den Hartog [6], the critical velocity for a cantilever and for 
a beam on fixed end supports can be obtained by multiplying 
equation (14) by 3.52/9.87 and 22.0/9.87, respectively. 

Case (d) Curved Plates With Fixed Edges. The assembly considered 
here consists of several plates whose long edges are fixed and 
whose cross sections normal to the direction of flow initially have 
a uniforin curvature. Cross sections of such an assembly are il- 
lustrated in Fig. 4 in the original and in the deformed conditions. 

Before deriving an expression for the deflection of the curved 
plate under the action of a uniform pressure differential, an ex- 
pression is needed for the bending moment per unit of beam width. 
The latter can be derived by using Castigliano’s first theorem, fol- 
lowing the procedure given by Timoshenko [4] for the case of 
centrifugal loading of a flywheel. The present case is analogous 
to that of a flywheel with rigid spokes and an elastic rim. Fig. 5 
shows a representation of half of a plate cross section and identi- 
fies some of the nomenclature used in the analysis. 

The following expressions were obtained by the author: 


pk 


hui 
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SEC. B-B Sec: C=C SEC: A 
Fig. 4 Outlines of plate assembly 
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Fig. 5 Diagram for curved-plate deflection analysi 
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No = ph 1 — 22, 
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M, = £ ( — cot a 
2p; a 


° ; cot @ 
m= or (1 = «2 


pR? sin? : 


2 ~~ pR? ( <3 
Bisina 26, \sina a 


N = pR (1 --— =) 


26; sin a 
where 
AR? 


B; = 7 fla) + fila) 


1 a sin 2a 
Ie) = osin? a (; sd -*) 
1 
fla) = fila) — 2a 
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The solution shown by Timoshenko includes the coefficient 2/; 
multiplying the right-hand side of equation (15). The author’s 
solution checks the flat-beam expression for X where @ approaches 
zero, whereas the solution shown by Timoshenko does not. 

The displacements of the curved plate in a direction normal to 
a midspan tangent define an approximate deflection curve from 
which the change in flow-channel cross section can be determined. 
This procedure neglects the relatively small displacements 
parallel to the specified tangent. Neglecting deformation due to 
shear and assuming the plate thickness to be very small in com- 
parison with the radius of curvature, it can be shown that the 
following defines the displacements normal to and relative to the 
midspan tangent to the plate as shown in Fig. 5 


Y 2R?M(1 — v?) ; — 
y= deve - ait. sin (; ? cos ddd 
/0 EI a 
¥ RN(1 — v?) 
T st sin yd@ 
0 EA 


For small angles the second integral may be neglected and the 


(25) 


first can be simplified to 


8 M(1 — v?) 
_ : - x'ds 
Jo EI 


which represents the moment-area method [5], modified by use 
of ds instead of dz as suggested by Roark [7]. Equation (26) 


(26) 


was used here to simplify the analysis. To obtain the displace- 
ment function y below, the y’-function was shifted as shown in 
Fig. 5. 

Thus 

sin? 


2 sin @ 


sin @ cos @ my 
(2% 
2 a 


Integrating equation (27) with respect to distance along the 
midspan tangent over the distance between the side plates 
and multiplying by two gives the change in coolant-channel 


6pR(1 — v?) | y sin y cos Y 
y= eens + — 
B, Ea* a a 


cross section. Dividing by bh, which for small @ and normal 
plate-thickness to coolant-channel-thickness ratios is a close ap- 


proximation to the initial coolant-channel area. gives 


eT r I sin 2a cos 2a 
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Substitution of equation (28) in (4) gives for the critical 
velocity A 
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Fig. 6 Critical velocity versus a for curved- 
plate assembly relative to that of flat-plate 
assembly; edges fixed 
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Fig. 7 Critical velocity versus b/a for 
curved-plate assembly relative to that of 
flat-plate assembly; edges fixed 
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The ratio of the afore-mentioned critical velocity to that for the 
flat plate with fixed edges is 
48; sin’ @ "h 


: (30) 
sin 2@ 
45 (. -—3 — + 2 cos 2a) 
a 


Vi, = 


Using A = aandI = a'/12 in equation (22) and expanding 
trigonometric functions in power series of @ with sufficient terms 
to give good accuracy for a < 0.5 gives 

9a?]'/2 
l —_ - 
14 | 


b?a? 5a 
Vi, = va 1 — —]|+¢+ 
15a? 42 


This critical-velocity ratio is plotted in Figs. 6 and 7. 

Case (e) Curved Plates With Hinged Edges. The assembly con- 
sidered here is identical with that of Case (d) except that the 
longitudinal edges of the plates are hinged to rigid side plates. 
It is assumed that flow channels are isolated from each other be- 
tween the ends of the assembly. 

Following the procedure used for the fixed-edge curved plates 
gives the following equations: 


(31) 


pk 


X = Bs 


(32) 
(33) 


c 
- (ese a — cot a) 


(34) 


a 
M = — (sin ? — sin? — 
B, sin a 2 2 


P E cos@ 
N = pR (1 _ 2B, ae 


where 
a Ak? 


By 7 


fla) + fila) (38) 


1 3. a cos 2a] a 
f(a) = ——— | a — = sin 2a + ———— (39) 
2 sin? @ 4 2 
and f;(@) is defined by equation (23). 

Use of equation (36) and the procedure used for the fixed-edge 
assembly gives the following expression for the displacements of 
the plate relative to the supports and in the direction normal to 
the midspan tangent 

pRML — v*)C 


1 eee aa 


4 
2EIp, sin a (o) 


where 


2 3 
—— + sin 2a|—cota—a 
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4 
* cos 2a + a cota ( 


C= 


Determining the fractional change in coolant-channel cross- 
sectional area as for the fixed-edge curved-plate assembly gives 
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Substitution of equation (42) in (4) gives the following expres- 
sion for the critical velocity 


vs | eee sin’ =|" 2 
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The ratio of this critical velocity to that for the simply supported 
flat-plate assembly is 


; 1/3 
Van = (Aa) (44) 


15C 


Expanding equation (44) in power series in @ with sufficient 
terms for good accuracy for a < 0.5 gives 


2b%a? a? lla? ]'/2 
Vv, = | — (1-—) 4+1-— 
* 5a? =) ia 21 


This critical-velocity ratio is plotted in Figs. 8 and 9. 

Dividing equation (29) by (43) gives the following ratio of the 
critical velocity of a fixed-edge curved-plate assembly to that of 
an assembly which is identical except for having hinged edges 


ba? 5a? 9a?]\*/2 
0 | (1 ir =) die -*| 
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This critical-velocity ratio is plotted in Fig. 10. 
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Effect of Unbalanced Channel Pressures 


If pressure differences are imposed in opposite directions across 
adjacent plates due to unbalanced conditions, the deflections due 
to the imposed pressure difference will be magnified by an ap- 
proach to critical flow velocity. The significance of unbalanced 
pressures can be illustrated by assuming that in the plate as- 
sembly of Case (a) at velocities below critical there are superim- 
posed pressure differences (yo across each plate) which cause 
small deflections of the form shown in Fig. 2. Adding po to equa- 
tion (4) gives the total pressure differential acting on the plates 
in the deformed region 

2pV.27AQ 


P = Po (47) 
Combining equations (3) and (47) to eliminate p gives an equa- 
tion which reduces to 
a Po 
Qo _ 30Ea*h Vo?pb(1 — v?) 
(1 — v?)b4 15gEa*h 





which on substitution of equation (4) becomes 
40. Po 
Qo 30Eath — , Vo \? 
(1 — v?)b4 V. 
Equation (48) shows that deflections due to initial pressure dif- 
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(48) 


Aprit 1960 / 87 

















CRITICAL VELOCITY 











RELATIVE 




















= oo * 1 | 
0.| 0.2 03 04 06 
a, RADIANS 








Fig. 8 Critical velocity versus a for curved-plate assembly relative to that of flat-plate 
assembly; edges hinged 














VELOCITY 








CRITICAL 


se amma 
| 





> 
wi 
a 
x 





—E——EE 

















100 


oT 


on i 
60 
b/a 


Fig. 9 Critical velocity versus b/a for curved-plate assembly relative to that of flat-plate 
assembly; edges hinged 
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Fig. 10 Critical velocity of fixed-edge curved-plate assembly relative to that of hinged- 
edge curved-plate assembly 
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ferentials across the plates are increasingly amplified as the flow 
velocity approaches the critical value. Note that the term 
[1 — (Vo/V,)?] is the ratio of the net stiffness due to combined 
fluid and elastic effects to the elastic stiffness, and the term 
(Vo/V,)? can be visualized as a normalized negative spring con- 
stant or stiffness factor. Actually, the assumption of alternately 
high and low pressures in successive channels is academic. If only 
a single plate in the assembly is subject to an imposed pressure 
differential, the deflection will be amplified to an extent smaller 
than that indicated by equation (48). 


Effect of Membrane Compressive Stresses 

If flat fuel plates are subjected to transverse loadings which 
cause membrane compressive stresses in the plates, collapse can 
occur due to these alone. Such a condition can result from pres- 
sures on the outside of the side plates exceeding the pressure in 
the channels between plates. The existence of a transverse com- 
pressive loading in combination with the effects of coolant flow 
allows collapse of the plates with less extreme conditions than 
those which individually would cause collapse. The analysis of 
this problem can be simplified by assuming that the local pressure 
difference across a plate is proportional to the local deflection in- 
stead of the previous assumption of uniform velocity at any 
transverse cross section of a channel. An expression for the 
critical velocity of an assembly of flat plates with continuous longi- 
tudinal-edge supports which impose a uniform compressive load 
on all plates was obtained by equating the bending-strain energy 
of the plate for a particular deflection to the sum of: (a) The 
work done by the fluid-pressure difference produced by the de- 
flection, and (b) the work done by the compressive load. The ex- 
pression for critical velocity as a function of compressive loading 


was then converted to dimensionless form to give 


(% ae 
P. -) 5 


(49) 


where P/P, is the ratio of the actual transverse compressive load 
to that which alone would cause elastic instability, and V,'/V, is 
the ratio of critical coolant velocity to that which corresponds to 
flow-induced collapse alone. The critical velocity V, for the 
velocity distribution assumed for this analysis is 10 per cent lower 
than that predicted for flat plates with uniform local velocity in 
the deformed channels. 

Equation (49) can also be deduced by visualizing that the two 
instability phenomena each produce a component of negative 
bending stiffness, just as an axial compression of a column in- 
creases the transverse deflection of a column under the action of 
a transverse force. When the sum of the negative stiffness due to 
fluid effects (proportional to Vo?) plus the negative stiffness due 
to transverse compression (proportional to P) is equal and op- 
posite to the positive stiffness of the plate, the system is unstable. 

If only one plate is in compression due to thermal expansion or 
welding residual stresses, the situation is approximately equivalent 
to the case of a single plate centrally located within a rectangular 
duct. The critical velocity of the single plate without the com- 
pressive load is easily seen to be higher by the factor +/2 than 
the critical velocity for the assembly of several plates with the 
same coolant-channel dimensions. This factor is appropriate 
because, based on the assumptions specified previously, a single 
plate centrally located within a rigid duct would be equivalent to 
any plate in a multiple-plate assembly in which the coolant- 
channel thickness is equal to twice the clearance between the 
single plate and the inside wall of the duct. The condition for in- 
cipient collapse of the single plate under compression is then 
approximately as follows in the region where the plate is near 
collapse due to membrane compression alone 
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where V.is computed for the assembly without consideration of 
the compressive load on a single plate. Fig. 11 shows plots of 
equations (49) and (50) and a curve which shows a plot based on 
judgment to represent the stability of the single plate in the 
region where both phenomena contribute to instability. 


Fig. 11 
plates 


Discussion 

The assumptions made for derivation of the critical velocities 
are considered not to impair the applicability of the critical- 
velocity formulas to practical fuel-plate assemblies with the ex- 
ception of factors discussed ‘below. 

The assumptions of initial flatness of plates, uniformity of plate 
thickness, spacing, mass flow, and pressure distribution lead to an 
overestimate of the collapse velocity if nominal parameters are 
used for calculation. The assumption for Cases (b) and (c) that 
there is no leakage from,one channel to another may lead to an 
underestimate of the collapse velocity, but it is not obvious that 
such leakage will have a strong effect on the critical velocity. 

The assumption of incompressibility of the coolant does not 
impair the applicability of the derivations to systems using liquid 
coolants and perfect gases insofar as the inception of instability is 
concerned. 
elastic deflections which are developed by velocities which are 
slightly above that at which instability begins. 

As the plate length-to-width ratio is reduced from large values, 
it appears that the critical velocity will increase. For many 
nuclear reactor applications this ratio is likely to be so large as 
not to have a large effect on the critical velocity. 

It may be possible to have a plate assembly in which the col- 
lapse modes of Case (b) and Case (c) would occur at velocities on 
the same order. Under this condition it appears that the lower 
of these two critical velocities, predicted on the basis of formulas 
presented here, would be a substantial overestimate of the actual 


Compressibility in a gaseous coolant may limit the 


lowest critical velocity. 

A small error was introduced in the curved-plate critical- 
velocity formulas by eliminating the second integral in equation 
(25). For a = 0.5, the fixed-edge curved-plate critical velocity 
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indicated by the formulas derived is overestimated by about 6 
per cent, and the hinged-edge curved-plate critical velocity is 
overestimated by about 4 per cent by using this approximation. 
The error is smaller with smaller values of a. 

It has been shown that transverse compressive loads on flat 
plates reduce the critical velocity. It is also clear that transverse 
tensile loads on all plates would increase the critical velocity. 

There may be other sources of instability which may interact 
with column instability and flow-induced instability. One of 
these is thermoelastic instability which has been analyzed by the 
author [8, 9] and independently by Iliffe [10]. Thermoelastic 
instability may arise due to the effects of fuel-plate deflection 
upon flow distribution and heat transfer, which, in turn, can cause 
temperature gradients which tend to maintain the deflection. It 
may also be possible that high neutron-flux gradients introduce 
an additional contribution to thermoelastic instability. In the 
idealized case a straight beam is unstable thermoelastically if, due 
to its temperature distribution, it tends to acquire a sufficiently 
high local curvature proportional to the local deflection. If this 
local curvature per unit of deflection of a beam of length L is de- 
noted K, the critical value of K is (#/L)? for a beam with simple 
end supports and (27/Z)? for a beam with fixed ends [9]. 

In the case of curved plates the presence of unbalanced pres- 
sures imposed across the plates introduces the possibility of 
buckling independently of the effects of plate deflection on the 
pressure difference existing across the plates. There is also an 
asymmetrical mode of flow-induced collapse which for large values 
of @ may have lower critical velocities than those presented here. 
In this mode of collapse the two half-widths of plate would de- 
flect in opposite directions. This mode of collapse at worst would 
limit the minimum critical velocity of the simply supported 
curved-plate assembly to about four times the critical velocity of a 
similar flat-plate assembly, because the more complex mode loses 
the membrane restraints which are present in the simple mode of 
distortion. The second mode of instability at worst would limit 
the critical velocity of a fixed-edge curved-plate assembly to some- 
what more than four times the critical velocity of a similar simply 
supported flat-plate assembly or somewhat less than four times 
the critical velocity of a similar fixed-edge flat-plate assembly. 
Actually the effect of this mode of collapse is probably much less 
than is indicated by the ‘‘at-worst’’ limits indicated in the fore- 
going. 

In view of the dependence of critical velocity on Young’s modu- 
lus of elasticity, it is apparent that the occurrence of yielding or 
creep of fuel plates due to unbalanced pressure, thermal stresses, 
and growth of fuel meat reduces the plate stiffness and the critical 
velocity. Roberts and Cottrell [11] have shown that alpha- 
uranium creeps during exposure to neutrons under external loads 
which without neutron flux would not produce significant creep. 
Willis [12] has suggested that alloys and dispersions of fissionable 
material will have a similar loss of resistance to creep deformation 
under external load. Thus it appears likely that in plate assem- 
blies, which are not afforded stability against collapse by the 
stiffness of cladding alone (neither shear nor bending stiffness 
attributed to the fuel meat), collapse may proceed slowly over a 
period of months or years instead of occurring almost instantane- 
ously. If creep collapse is a possibility, it would be prudent to 
remove a fuel-plate assembly from the reactor at intervals for 
inspection until it is confirmed that creep collapse does not occur. 
Tf the fuel-meat were assumed to have no shear or bending stiff- 
ness, the critical velocity of the assembly of Case (a) would be 


y 30gEt,th “]'/2 
xii pb4(1 — pv?) 


where ¢, is the thickness of cladding on each surface. 
There would probably be a slight loss of critical velocity of 
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fixed-edge plates due to local yielding in cladding at the longi- 
tudinal boundary of the fuel meat due to growth of the meat with 
burnup of fissionable material. The zone of yielding would 
probably be so localized as to make a slight deviation of the 
fixed-plate critical velocity toward that for simply supported 
plates. It is interesting to note that the significance of this 
problem is smaller for curved plates than for flat plates, since for 
the former at large values of @ the critical velocities coincide. 

The assumption of uniform coolant velocity in a channel across 
a plane normal to the flow is slightly optimistic for Cases (a), (b), 
(d), and (e). This is indicated by the analysis leading to equa- 
tion (49), which assumed that the local velocity was inversely 
proportional to the local spacing between plates. This assump- 
tion gave a critical velocity for flat plates 10 per cent lower than 
did the assumption of uniform velocity. Actually, the velocity at 
the region of minimum plate separation may be smaller than the 
average velocity in the channel. This condition would make 
the critical velocities somewhat higher than the derived expres- 
sions indicated. 

In some fuel-plate assemblies collapse has occurred at velocities 
on the order of one half of that predicted by the formulas pre- 
sented here. The calculated critical velocity for the ETR fuel- 
plate assembly is about 54 fps for water at room temperature. 
Observed collapse velocities are reported by Doan [2] to be as 
low as 30 fps. The fuel element consists of 19 aluminum-clad 
uranium-aluminum alloy plates 0.050 in. thick, about 2.62-in. 
span, with water-channel thickness of 0.108 in. and length of 36 
in. Doan cites lack of flatness of some plates as contributing to 
collapse and also points out that the plates would deform plas- 
tically with pressure differences across the plates of one and one- 
half to two psi for the outermost plates and about three psi for 
the inner plates. He also reports that use of work-hardened 
plates or perforation of the side plates allowed application of 
flow at 25 per cent or more above the design flow rate of 35 fps 
with little or no distortion. With the 25 per cent variation in flow 
distribution indicated by Doan, the work-hardened plates have 
attained, without collapse, velocities almost as high as the 
theoretical collapse velocity. 

The data on the Oak Ridge collapse observations [1] are not 
sufficiently complete to provide a check on the derived collapse- 
velocity formulas. The MTR plates are made of aluminum-clad 
urarium-aluminum alloy, 0.059-in. plate thickness, 0.117-in. 
water channel, 2.67-in. free span, 5!/:-in. radius of curvature. 
The plates are brazed into side plates. The theoretical collapse 
velocity for these plates with water at room temperature is about 
240 fps. With an assembly with the same dimensions except for 
plates being flat the theoretical critical velocity would be about 80 
fps. Reference [1] reports collapse of 1.0-rnm (0.039-in.) and 
1.5-mm (0.059-in.) curved plates and 1.0-mm flat plates, 2.68 in. 
wide. Water-channel thicknesses were not identified for specific 
assemblies, but a normal dimension of 3 mm was specified. 
Collapse velocities were not specified, but the test series, including 
assemblies which did not collapse, was run up to velocities of 
about 45 fps. The theoretical critical velocity for the 1.0-mm 
flat-plate assembly is about 44 fps so it is not surprising that it 
collapsed. Collapse of the other plates may have been influenced 
by deviations from ideal conditions. 

These experimental data lend some support to the formulas 
presented here. It would be of considerable interest to have tests 
conducted to determine the validity of the critical-velocity 
formulas for idealized assemblies and to investigate the effect of 
derivations from conditions assumed for derivation of the formu- 
las. It is also essential that the effect of irradiation upon the re- 
sistance of the fuel meat to distortion be evaluated. 

In development of fuel-plate assemblies it is appropriate to 
conduct flow tests at velocities substantially above rated flow to 
confirm the existence of a substantial margin against collapse at 
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rated flow. This margin should be made large enough to allew 
for differences between test conditions and in-pile conditions in- 
cluding fuel burnup. 

In view of the reduced plate stiffness introduced by the ap- 
proach to collapse velocity, it would be of interest to know the 
effect of this reduced stiffness upon the natural frequencies of vi- 
bration of the plates. In the lowest-frequency mode of bending 
vibration of the plates the flow-induced negative stiffness would 
be a factor only for the case where adjacent plates vibrate in 
phase opposition. 
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Irwin Beretsky® 


Correction Factors to Permit Application of Critical Velocity Deter- 
With Initial Local Deviations in Passage Area 


tion to P 9 





Consider as before a series of parallel plates that are spaced 
uniformly along their axial position except for a local position 
AA’ where they are initially deflected as shown: 
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Wr , _ . 
Assume that y’ ~ D sin b which satisfies the end and middle 


conditions of the deflected plate as shown. 

With this in mind, proceed as in the paper to postulate a de- 
flection curve for some initially flat plate with a particular type of 
end fixity. 

As an example, we can consider the simply supported flat plate 
subject to a uniform pressure loading. The deflection y can be 
written as 

Px(1 — v?)(b? — 2bx? + x5) 


(51 
y 24EI1 ~ 


If we take equation (51), integrate it from 0 to b, and multiply 

it by two, we will determine the change in area of the channel, 8. 
P(L — v?)b§ 

Pa sans (52) 

60EI 


and since J per foot of plate = a8/12 


Pi — v?)b* 


5Ea* 


8 (53) 


Applying Bernoulli’s theorem, the pressure differential P is 


-? = \(j : A\ 43 (; a a 


Differentiating equation (53) and (54) to obtain dB/dP we have 


from (53) 


dB (1 — v2)b5 


dP —s*OBEa® 
dB - g (A? — B)8 


(55) 


= - (56) 
dP 2A%V*%A2 +4362) 13 


and (54) from 


Since the variation of 8 with pressure in equation (55) is equal 
to that of equation (56) we have 


7” 5Ea'g (A? 7 ey 57) 
Emme Ray Yall, Seer ms (O4 
2po(1 — v2)b§ A® | A? + 362 

Multiplying the right-hand side of equation (57) by bh/A = 
we find 
5Ea*h 1(A?2 — 6?) 
V2 = Bnd or B ) t (58) 

2p(1 — v2)b AA? + 36?) 

This equation will reduce to the author’s equation for areas (8) 

very small compared to A. 

However, equation (58) can be written as 


V=eVC (59) 


Where V, is critical velocity as previously determined and C is 
the correction factor 


1 | (A? — Bs 


_- —— (60) 
A? ¥ (A? + 36?) 


oC = 


Since in the beginning of this section we assumed an initial de- 
flection y’ ~ D sin m2x/b, we can calculate the @ that initially 


exists. 
b 
T 46D 
2 D sin jie = (61) 
0 


4] T 


Noting also that A = bh, and substituting into equation (60) 
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Fig. 12 Correction factor C versus D/H 
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If we let D = ah where 0 < a@ < !/2 we have finally 


1 | (a? — 16a?)3 
Tr? T? + 48a? 


(64) 


This is the correction factor which is plotted in Fig. 12. 
Nore: All nomenclature used except where noted, is Mr. 
Miller’s original nomenclature. 


E. B. Johansson* 

In this paper it is assumed that the plates are long enough com- 
pared with their width that they can deflect locally without sig- 
nificant redistribution of flow among the coolant channels. 

It. is also assumed that the difference between the friction drops 
in the constricted and opened channels, in the deflected region, is 
small. In many reactor plate designs the length of the plates is 
of the order of 6 to 15 times the plate width. Any local deflection 
of the plates must extend over a length of the order of two plate 
widths, or there would be a significant increase in the stiffness 
of the plates for this deformation because of axial bending. 
Therefore in many practical cases the deflected region is an 
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appreciable portion of the total length, and some redistribution of 
flow is to be expected. 

An analysis has been made to determine the effect of this flow 
redistribution on the critical velocity. A region of length equal 
to two plate widths, at a distance (z) from the inlet end of the 
plates, was assumed to be deflected in the same manner as in 
Miller’s paper. In this region the friction pressure drop per unit 
length in the constricted channels is higher than in the unde- 
flected portion of the channels. This occurs because of the local 
increase in velocity and the decrease in hydraulic diameter of the 
channel. In the channels which are expanded locally, the situa- 
tion is reversed and the local friction drop per unit length is lower. 
Therefore the flow through the constricted channels is decreased 
and the flow through the opened channels is increased. 

The flow redistribution was determined by imposing the con- 
dition that the pressure drops through all channels are equal. This 
pressure drop is the sum of the entrance loss, the friction drop in 
the deflected portion of the channel, the friction drop in the unde- 
flected portions, and the exit loss. With the velocities deter- 
mined and the variation of pressure drop over the channel length 
known, the static pressure difference acting on the deflected 
plates is readily found. This static pressure difference is given by 
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static pressure difference acting on deflected portion 
of plate 

pressure difference from inlet to exit of plate assembly 

inlet loss coefficient 

exit loss coefficient 

diffusion loss coefficient for deflected region 

length of deflected region 

original flow area of channel 

change in flow area of channel 

length of plates 

friction factor 

distance from inlet to beginning of deflected region 


Q 
AQ 
” BD 


(65) 


where 


Ap, 


This relation is valid for large deflections provided that the plates 
do not yield. For small deflections, Ap, increases approximately 
linearly with AQ/Q, but as AQ/Q increases the curves of Ap, 
versus AQ/Q, tend to flatten out. This is illustrated in Fig. 13 
where the results of calculations for a particular plate geometry 
are presented. The plates were assumed to be clamped along 
their edges; this defines the deflection curve and the relation be- 
tween AQ/Q) and the deflection. The effect of the axial position 
of the deflected region should be noted; a deflection downstream 
from the inlet generates a smaller pressure difference than a de- 
flection at the inlet. 

A family of straight lines which represent the stiffness charac- 
teristics of the various thicknesses of plates are also plotted on the 


same figure. These lines are simply plots of 


_ 1 Bp — *) 
384 EI 


and represent the pressure difference that a plate can equilibrate 
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(1) STATIC PRESSURE DIFFERENCE ACTING ON 
DEFLECTED PLATES (SOLID CURVES) 
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LIGRATS PLATE SPRING FORCE (DASHED 
CURVES) FOR SEVERAL PLATE THICKNESSES 
VERSUS LOCAL PLATE DEFLECTION 
Fig 13 
at a given deflection. For a given plate thickness, undisturbed 
channel velocity V> and position of deflected region, the plate is 
stable if the plate characteristic lies above the hydraulic force 
curve [equation (65)]. It is possible that the plate characteristic 
may intersect the hydraulic force curve. In this case the analysis 
predicts a stable configuration with the plates partially collapsed 
and having a deflection corresponding to the point of intersection 
of the two curves. 

If one defines the critical velocity by the condition that tbe 
plate characteristic and the hydraulie force curves have the same 
initial slope, a fairly simple expression can be obtained for the 
This velocity is that at which instability is 
In equation (65) the terms 


critical velocity. 
initiated in the region in question. 
in AQ/Q) are expanded in a series and higher order terms are 
dropped. After some algebra one obtains 

Ap AC 3fl 

DP Ww Jy 14 

ae “S 
pV? 


6 l [ +14 f(x + */el) ]\ 66) 
== 4 (00 
Yo 2 5 lee tee 


ky + ky 


(fl 
2h 
Miller obtains, neglecting flow redistribution 
Ap, 
eg 


“? o 
9 PV" 7 


where I 
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Comparing these two relations, the ratio of critical velocity 
considering flow redistribution to the Miller critical velocity is 


Ver’ fy etl 


\: PY a fest 
4 2h 21+ RR) L 


, 


rae oh 


(67) 


This equation is plotted in Fig. 14 for the same plate geometry 
used for Fig. 13. It is seen that, for this particular geometry, the 
critical velocity considering flow redistribution is within 10 per 
cent of the Miller critical velocity for entrance collapse. For 
collapse downstream from the inlet, the critical velocity is con- 
siderably higher. 














1.0 


Fig. 14 Ratio of critical velocity considering flow redistribution to 
critical velocity without flow redistribution versus distance from inlet to 
beginning of collapse region 





It should be noted that as //L — 0, the critical velocity ap- 
proaches a value smaller than the Miller critical velocity. This 


occurs because of the term - The friction drop in the con- 
stricted channel is higher than in the opened channel in the de- 
flected region. This generates a pressure difference which adds to 
the venturi effect. 

Fig. 15 shows schematically the variation of total pressure 
along the channel for various positions of the deflected region. 
The static pressure difference acting on the plates is given by the 
difference in the total pressures minus the difference in dynamic 
pressures. It can be seen from the plots that the static pressure 
difference is highest at the inlet end of the plate assembly. In 
addition, for deflection at intermediate points, the static pressure 
difference is higher at the downstream end of the deflected region 
than at the upstream end. Therefore there is a tendency to 
move the deflected region downstream. 

The close correspondence of the critical velocity considering 
flow redistributions and Miller’s critical velocity at the inlet oc- 
curs because there are two compensating factors at the inlet: 


1 The flow in the constricted channel is reduced, thus re- 


ducing the venturi effect. 
2 The friction drop is more rapid at the entrance in the con- 


strieted channels than in the opened channels (see Fig. 15). This 
means that the total pressure (static plus dynamic) is lower in 


the constricted channels. 
The net result is that the static pressure difference acting on 


the plates for the geometry considered is not greatly changed by 
the flow redistribution. In other geometries these two effects 
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may not so nearly cancel. The length of the collapsed region was 
taken as two plate widths. If a shorter distance were assumed 
the plates would be stiffer because of axial bending restraint, with 
1 consequent increase in calculated critical velocity. If a longer 
collapse region were assumed, the calculated critical velocity 
would be higher because of the effect of*/ in the equation for 
Ver'/Ver. Thus it is considered that the assumption of two plate 
spans gives to a very good approximation the minimum and cor- 
rect value for Vor’. 

If a support comb is used at the entrance, midspan entrance 
collapse cannot occur, and the collapsed region must start some 
distance downstream, perhaps one plate width. The following 
remarks are speculative. Now suppose that the velocity is 
slightly above the critical velocity for collapse at this position. 
The plates begin to deflect locally. The static pressure difference 
is highest on the downstream end of the collapsed region. This 
tends to extend the collapse region downstream. However, this 
increase in the extent of the collapsed region increases the critical 
velocity, perhaps to the point where the static pressure difference 
at the upstream end of the collapsed region is not sufficient to 
maintain the defiection. In this case the deflected region would 
move downstream. As it moves downstream the critical velocity 
increases, and the deflection may vanish before it reaches the end 
of the plates. 1f this model is correct, one would predict possible 
vibration of the plates with an inlet comb at velocities above the 
critical, and collapse would not occur until much higher velocities 
are reached. 

Experimental data have confirmed the collapse velocity pre- 
dicted by the author’s formulas in the case of no entrance sup- 
port comb. However, collapse has not been observed with an 
entrance comb at velocities as high or higher than twice the 
theoretical velocity. 


94 / APRIL 1960 


The following proportions and flow coefficients were used to 
construct Figs. 13 and 14: 


l/h = 11.4 


L/h = 50 UL = 0.23 
0.05 


f = 0.0156 


Author's Closure 

The author appreciates the interest shown by Mr. Beretsky 
and by Mr. Johansson. Their discussions deal with important 
aspects of the problem of fuel plate collapse. 

The author has analyzed the situation discussed by Mr. 
Beretsky and disagrees with his The error in 
Mr. Beretsky’s analysis lies in his lack of justification for equating 
his two expressions for d8/dP. Only when the plates are initially 
flat does Beretsky’s procedure correctly determine the conditions 
for incipient collapse. Actually, to represent a condition of static 
equilibrium, the pressure differential developed by flow with an 
elastic deformation superimposed on the initial deformation 
should be equated to the pressure differential which is required 
to maintain the elastic deformation. 

Consider an assembly of plates with long edges simply sup- 
ported in which at some local region along the length of the as- 
sembly there is a deformation such that the area of the coolant 
channels is changed an amount AQ,, adjacent plates being de- 
formed in opposite directions. If the nomenclature and assump- 
tions previously used by the author are retained and if AQ, de- 
notes the change in local channel cross-sectional area due to flow- 
induced pressure difference, the following relation can be derived: 


Vo? 1 . | ‘ 
iS ead ere) 
29 l—n l+n 


_ AQ + AQ 
7 Qo 


results. 


(67) 


where 


n (68) 
Note that Qo is the channel cross-sectional area in the undeformed 


regions. By substituting AQ. for 8 in equation (53) we have 


5Ea® AQ, 


p= ; (69) 
ti - p?)bs 

The pressure differentials defined by equations (67), (68), and 

(69) must be equal at velocities lower than that corresponding to 

instability or collapse. Thus for velocities up to the critical 

value we obtain the following, after combining equations (6), 


(67), and (69) 
AQA1 — n?)? (P)" 
bhn ~ ee 


c/ 


(70) 


Substituting bh = Q> in equation (69), and substituting the re- 
sult in equation (70) gives 


AQ: E E (2 + Ae)" vy 
bh ("") 


oss 7 


“AQ. + AQ, (71) 


Equation (71) can be solved by trial and error to find values of 
AQ, for any value of Vo less than that for which AQ» goes to in- 


finity. If the initial deformation of the plate had a form cor- 
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responding to a uniform pressure differential, the value of (AQ 
+ AQ.) would be 16/25 when flow was sufficient to cause adjacent 
plates to touch. Under this condition the term in brackets in 
equation (71) reduces to 


E _ (A — ae)" = 0.59 
bh 


Note that equation (71) indicates that with no initial deforma- 
tion the flow required to maintain a deflection becomes smaller 
as the deflection increases. This effect was eliminated in the 
presentation in the paper by the assumption of small deflections. 
It is also offset by other factors which were not considered in the 
analysis here and in the paper. These factors include frictional 
effects and nonuniformity of flow velocity within the deformed 
channel. In a plate assembly with long edges fixed the ratio 
of channel area change to maximum deflection of plate is smaller 
than for the assembly with simply supported plates, and this 
factor would increase the value of the term in the brackets in 
equation (71) for a particular maximum deflection. Finally, 
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with deformations limited to values which are tolerable in a reac- 
tor the term in brackets in equation (71) is essentially unity and 
the equation reduces to 


A+ A 1 
ane) 
i 


Equation (72) with V, evaluated by using an appropriate equa- 
tion in the paper should give a satisfactory estimate of deforma- 
tions in the idealized cases considered by the author. 

Mr. Johansson’s discussion makes a valuable contribution to the 
theory of flow-induced deformation and collapse of reactor fuel 
plates. He has extended the work upon which he based his dis- 
cussion of the author’s paper and has prepared a paper [13] which 
has been submitted to ASME. 





(72) 


Additional Reference 


13. E. B. Johansson, ‘‘Hydraulic Instability of Reactor Parallel- 
Plate Fuels Assemblies.”’ 
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The Propeller Turbines of the 
St. Lawrence Power Project 


Although the St. Lawrence River has long been recognized as a great economical source 
of water power, its development has been delayed until now for many reasons. The 
major obstacle has been the 125 miles of river which forms the common boundary be- 
tween the United States and Canada. In this distance lies the International Rapids 
section of the St. Lawrence River, which necessitated a joint undertaking by both coun- 
tries if power was to be developed. The St. Lawrence Power Project is the fulfillment of 
this joint undertaking. This paper describes the control structure required to harness 
the river and the design features of the hydraulic turbines which convert the power of this 


G. G. FINTAK 


Application Engineer, 
Hydraulic Department, 
Allis-Chalmers Manufacturing 
Company, Milwaukee, Wis. 


mighty river to useful, dependable electrical energy. 


Bini OF THE MOST DEPENDABLE RIVERS in the world, 
the St. Lawrence River, whose maximum flow is only slightiy 
more than twice its minimum flow, will soon be harnessed to pro- 
duce dependable hydroelectric power. 

The Power Authority of the State of New York and its Cana- 
dian partner, the Hydro Electric Power Commission of Ontario 
are now completing the gigantic task of developing the power 
potential of the International Rapids Section of the St. Lawrence 
River. This section stretches from Ogdensburg, N. Y., to 
Massena, N. Y., a distance of approximately 48 miles. 

In order to develop the potential power of this section of the St. 
Lawrence River, it is necessary to use the drop in water level from 


[ owranio | 


ternational Boundary divides the structure into an American 
powerhouse and a Canadian powerhouse. It will act as a gravity 
dam 3300 ft long and 167 ft high, with an enclosed erection bay 
located at each end of the structure. Six ice sluices are incor- 
porated in the structure, two at the center astride the Inter- 
national Boundary, and two each under the erection bays. 

The exterior appearance of both powerhouses is identical with 
the cranes and service equipment being interchangeable. 

Upstream a few miles, the Long Sault Dam stretches from 
Barnhart Island to the United States mainland. It is a con- 
crete, gravity curved-axis, spillway dam, 2960 ft long and 114 ft 
high, having 30 spillway gates each 50 ft wide and 30 ft high. 


iNcvesioe (TN) O LONG SAULT 


INTERNATIONAL #4 


~ BOUNDARY 9 


| NEW YORK | 


St. Lawrence Power Project, control structures 


Fig. 1 


the eastern end of Lake Ontario to the powerhouse site near 
This drop is spread over 125 miles of river and is of 
By effective 


use of control structures, a head pond is created above the power- 


Massena. 


little value unless it is concentrated at one point. 


house which will produce an average head of 81 ft to operate the 
32 hydraulic turbines installed in the powerhouse. 

The control structures necessary to provide this operating head 
consist of a powerhouse and two dams, all of which are being con- 
structed without interfering with navigation on the present 14-ft 
Cornwall ship canal. Fig. 1 shows the control structures of the 
St. Lawrence Power Project. 

The powerhouse structure, called the St. Lawrence Power Dam, 
stretching from Barnhart Island, U.S. A., to the Canadian main- 
land, consists of two adjoining powerhouses in one, since the In- 

Contributed by the Hydraulic Prime Movers Committee and pre- 
sented at the Hydraulic Conference, Ann Arbor, Mich., April 13-15, 
1959, of THe AMERICAN Soctety OF MECHANICAL ENGINEERS. 

Nore: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those 
of the Society. Manuscript received at ASME Headquarters, 
November 21, 1958. Paper No. 59—-Hyd-11. 
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LO * 
SAULT DAM 


GRASS RIVER 


This dam and the St. Lawrence Power Dam combine to maintain 
the head-water level on the hydraulic turbines. 

Twenty five miles upstream another control dam, Iroquois 
Dam, spans 1980 ft from the United States to Canada. This 67-ft- 
high dam will regulate the outflow of water from Lake Ontario 
into the forebay produced by Long Sault Dam and the St. 
Lawrence Power Dam. It will have 32 gate-controlled sluiceways 
each 50 ft wide and 48 ft high. 

Two navigation locks near Massena, N. Y., and one lock at 
Iroquois, Canada, will allow river navigation to pass the afore- 
mentioned control structures. 

The St. Lawrence Power Dam is the largest single structure of 
the project and will house 32 hydroelectric units that will produce 
over 1,800,000 kw of electrical energy. Known as a modified out- 
door type, the powerhouse will not have the conventional super- 
structure over the generating room, but instead the units will be 
protected by removable covers. Fig. 2 is a typical cross section 
of the American powerhouse which shows the location of various 
powerhouse equipment. 

The vertical-shaft synchronous generators will have a name- 
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Fig. 2 American powerhouse, typical cross section 


plate rating of 57,000 kw or 60,000 kva at 0.95 power factor, 
13,800 volts, 3 phase, 60 cycles, operating at 94.7 rpm. 

The hydraulic turbines are of the vertical-shaft, fixed-blade- 
propeller type, with concrete semispiral casings. Of the 16 hy- 
draulic turbines in the American powerhouse, 8 are being manu- 
factured by the Allis-Chalmers Manufacturing Company, and 8 
are being manufactured by the Baldwin-Lima-Hamilton Cor- 
poration, The 16 hydraulic turbines in the Canadian powerhouse 
are being manufactured by the English Electric Company, Ltd., 
of Canada. Table 1 gives comparative ratings and data for the 
hydraulic turbines being installed in the St. Lawrence Power 
Dam. 

Because several manufacturers are supplying the hydraulic 
turbines under two distinct specifications, variations of turbine 
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design exist. The following description of the mechanical features 
of the hydraulic turbines will be that of the turbines supplied by 
the author’s company. 

The shape of the concrete semispiral casing which conveys and 
directs the water to the turbine propeller is shown clearly by 
means of a model, Fig. 3. The model has been sectionalized at the 
horizontal center line of the turbine distributor and only the 
lower half is shown. The gentle upward slope of the casing floor 
to the stay-ring cone can be seen in the water passage on the 
right. The actual inlet is 65 ft wide including the two seven-ft- 
wide concrete piers. 

Generally, the main parts of the St. Lawrence hydraulic tur- 
bines are of what might be called normal design for large fixed- 
blade-propeller-type turbines. Fig. 4 is a cross section through 
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Table 1 


Item 


Manufacturer Allis-Chalmers 


NO OE UNA. si05 sac eaws 8 
Guaranteed Performance: 
81 ft rated net hd, max 
IMIG DD «6 o'ne waco 79000 
710002 
87 ft max net hd, max 
output, hp.. a5 
75 ft min net hd, max 
output, hp 
Turbine speed, rpm..... 
Runner diameter, in... . 
Elevation of center line of 
distributor a 
Diameter of outer-head 
cover, in 
Servomotor capacity, ft-lb 


85000 


72000 
94.7 
er. 240.0 


154.0 


578000 


* Best efficiency points. 


Fig. 3 Model of semispiral casing 


the turbine, which identifies the various turbine parts described 
in the following paragraphs. 

The stay ring, which is the main structural member of the tur- 
bine, is made of cast steel in five sections, designed to be rigidly 
bolted together. This ring has ten ribs to transmit the super- 
imposed loading of the powerhouse to the foundation and to re- 
sist the internal pressure of the water-filled casing. The stay ring 
is approximately 32 ft diam with a distributor height of 103 in. 
Fig. 5 shows the stay ring assembled during manufacture. The 
crane-mounted radial drill press being used was designed for 
machining large pieces, such as this stay ring on the erection 
floor. 

A cast-steel bottom ring, containing the bronze-bushed 
bearings for the lower trunnions of the wicket gates, is shop- 
welded to a fabricated plate-steel discharge ring to form an inte- 
gral ring. The combination ring is sectionalized in halves and 
rigidly bolted together. The ring is straight bored to provide a 
running clearance with the propeller runner of 0.125 in. Re 
newable carbon-steel facing plates are fitted to the ring adjacent 
to the lower ends of the wicket gates. 

The turbine head cover consists of a combination of three 
annular covers; namely, the outer, intermediate, and inner-head 
covers. With this combination of covers the propeller, runner, 
main shaft, guide bearings, inner, and intermediate-head covers 
can be removed through the generator-stator bore without re- 
moving the outer-head cover and wicket gates. The outer-head 
cover is made of cast steel in four sections rigidly bolted together. 


1960 
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American powerhouse 


St. Lawrence power project, turbine data 


Canadian 
powerhouse 
Baldwin-Lima English Electric Co. 
Hamilton 


8 16 


77500 
710002 
82400 
70500 


94.7 
229.0 


154.0 


308.0 
625000 


It is 27 ft diam and is provided with two bronze-bushed bearings 
for the upper shafts of each wicket gate. Renewable carbon-steel 
facing plates are fitted to the outer-head cover adjacent to the 
upper ends of the wicket gates. The intermediate-head cover is 
formed and fabricated of plate steel. The main guide-bearing 
housing and the gate-operating ring are supported by this cover. 
The inner-head cover is fabricated of plate steel and supports the 
main-shaft water seal. 

The 20 wicket gates, which regulate the flow of water to the 
propeller runner, are formed and fabricated of plate steel welded 
to a steel shaft. Longitudinal machined surfaces are provided on 
each wicket gate along the areas that will make contact when the 
gates are closed. An adjustable cap at the top of each wicket-gate 
shaft allows the gates to be centered in the turbine distributor to 
maintain a 0.0625-in. top and bottom clearance between the 
outer-head cover and the discharge ring. 

A two-piece lever is keyed to the upper shaft of the wicket gate 
with overload protection provided in each lever by means of 
shearing pins. The levers are connected to the gate-operating 
ring by means of adjustable links. 

The gate-operating ring is made of fabricated plate steel in one 
piece. The gate-operating ring moves on grease-lubricated bronze 
ways located on the intermediate-head cover. The force to 
operate the wicket gates is obtained by two 26-in-diam oil-pres- 
sure servomotors which are mounted in the turbine-pit liner. 
Steel connecting rods transmit this force to the gate-operating 
ring, which in turn, by means of the links and levers operates the 
20 wicket gates simultaneously. 

The main shaft is a hollow-bored, steel forging, 37 in. diam, 
with integral forged flanges at each end for bolting to the runner 
hub on the lower end and to the generator shaft flange on the 
upper end. A renewable stainless-steel sleeve is attached to the 
shaft where it passes through the water seal. 

The cast-steel propeller runner consisting of six blades and hub, 
is of particular interest as it is too large to be shipped as an in- 
tegral casting. Therefore, the blades of the 240-in-diam propeller 
runner were cast separately and arranged for bolting to the hub 
in the field. Ten studs in each blade trunnion permit an accurate 
adjustment of the 25-deg blade tilt in the field prior to the final 
doweling of the blades to the hub, Fig. 6. To provide sufficient 
bolt diameter to carry the large thrust and torque loading on the 
overhanging blades and to prevent nut interference inside the 
hub, two bolts were omitted from the bolt circle at a point where 
the blade trunnion is in compression under load. 

Each propeller blade was prewelded on the bottom side with 
stainless steel '/, in. thick on an area (approximately 15 sq ft per 
blade) most likely to be subjected to pitting due to cavitation. 
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Fig. 4 Turbine cross section, Allis-Chalmers turbine 


This area is clearly shown in Fig. 6. A stainless-steel strip '/2 in. 
thick was welded to the periphery of the blade to protect the 
running clearance between the blades and the discharge ring. 

In order to obtain a uniform deposition of stainless-steel over- 
lay and to reduce the cost of grinding the blade surface after weld- 
ing, the overlay was automatically placed by means of a welding 
machine and positioner which automatically maintains the point 
of weld level in respect to the welding head. Fig. 7 shows the 
automatic welding machine and positioner in operation, 

The propeller runner assembly was statically and dynamically 
balanced in the shop by means of a balancing machine, Fig. 8. 
The balancing is done in three steps: (1) Each blade is weighed 
before assembly and its moment about the center line of the main 
shaft determined. The blades are so arranged in the hub as to 
give the least possible unbalance. (2) The runner assembly is 
statically balanced on the balancing machine and the necessary 
corrections made. (3) The assembly is revolved at 10 rpm for 
dynamic balancing. 

The main guide bearing, the only main-shaft bearing on the 
turbine, is of the babbitt-lined, oil-lubricated, pressure-feed type. 
Two motor-driven gear pumps, one ac and one stand-by de, lo- 
cated on the intermediate-head cover, circulate the oil from an 
annular reservoir below the bearing to the bearing surface. 
Circulating-water-type oil coolers are provided in the oil system. 
Thermostatic-controlled heaters are provided in the oil reservoir 
for maintaining the oil temperature when the turbine is idle. 

Below the main bearing and adjacent to the stainless-steel 
sleeve on the main shaft is a conventional packing box to seal off 
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leakage water along the shaft. Standard square packing is used 
with grease and clear water lubrication. 

Below the packing box is a hollow rubber sealing tube surround- 
ing the main shaft and normally clear of the shaft. When the 
turbine is at rest, this sealing tube can be inflated to seal against 
the shaft, permitting the packing box, which is below low tail- 
water level, to be serviced without unwatering the turbine. 

To improve the turbine performance when operating at small 
wicket-gate openings, a cam-operated turbovent allows free air to 
be admitted through the runner hub to the turbine draft tube. 

Should it be desired to motor the unit, provisions have been 
made to admit compressed air below the outer-head cover. With 
the wicket gates closed, it will be possible to depress the tail-water 
level in the draft tube below the level of the runner blades and 
allow the runner to revolve in air, reducing the power require- 
ments for motoring the unit. 

A centralized automatic forced-grease lubrication system is pro- 
vided to lubricate the moving parts of the turbine. The proper 
amount of grease is supplied to each of 135 lubrication points by 
means of measuring valves. Two piping systems are provided, 
each connected to separate time-controlled grease-distributing 
valves. One system, operating on a short cycle, lubricates the 

lower and center wicket-gate shaft bearings (40 points) which 
are subject to water pressure. The second system, operating on 
a long cycle, lubricates the remaining points of lubrication. The 
length of the lubricating cycles will be adjusted to suit field 
operation. 

An elbow-type draft tube is formed in concrete with the upper 
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Fig.7 Automatic we!ding machine and positioner for applying stainless- 
steel surface protection 
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Fig. 9 Shop assembly of turbine looking down into the turbine pit 


portion, which is below the discharge ring, lined with plate steel. 
The draft-tube top liner contains a mandoor to allow access to 
the turbine below the runner when the unit is dewatered. 

The turbine pit opening above the headcover is lined with plate 
steel. Bases to support the wicket-gate operating servomotors are 
located in the turbine pit liner. 

Fig. 9 is a shop erection view of the St. Lawrence turbine look- 
ing down into the assembled turbine. The gate-operating ser- 
vomotors, the yt pr ring, the intermediate-head cover 
with the main-guide bearing oil pumps, the main-shaft guide 
bearing, turbine pit liner, and stay ring are clearly shown. 

Fig. 10 is another shop assembly view, looking into the water 
passages of the turbine. The stay ring, pit liner, intermediate- 
head cover, and several of the wicket gates are shown. 

The 16 hydraulic turbines in the American powerhouse are 
controlled by Allis-Chalmers twin system, oil pressure, cabinet- 
type actuator governors, Fig. 11. The twin governor cabinets 
are located on the generator floor between each pair of units. 
Each governor cabinet contains two actuators, two motor-driven, 
gear-type oil pumps, oil sump tank, and all control mechanism 
and devices. Two oil-pressure and two air-pressure tanks are 
located behind each cabinet at the same level. The normal oper- 
ating oil pressure is 300 psi and each 390-gpm oil pump is capable 
of supplying the required oil for the two actuators with the second 
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Fig. 11 Twin cabinet governor, front view 
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Fig. 12 Model performance curves at various blade tilts for the 
St. Lawrence fixed-blade propelier turbines 


pump standing by. The governor is of sufficient capacity to move 
the wicket gates a full-operating or full-closing stroke in 8 sec. 

The governors are equipped with the necessary apparatus for 
adjusting the load limit, speed drop, and speed level, and for 
automatic shutdown, remote control, remote indication, and auto- 
matic intermittent brake application. 

A permanent magnet generator mounted above the generator 
exciter and driven by the main shaft of the unit, electrically 
transmits turbine speed to the flyball motor and thence to the 
flyballs, the speed-sensing element. 

Before design and manufacture of the turbines began, extensive 
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Fig. 13 Expected performance curves for the St. Lawrence fixed-blade- 

propeller turbines 
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Fig. 14 Field performance curves for the St. Lawrence fixed-blade 
propeller turbines 


laboratory research was made on an homologous turbine model 
which included intake, scroll case, turbine, and draft tube water 
passages. Model runners, approximately 12!/2 in. diam were 
tested at heads of 68 to 87 ft and showed a maximum model ef- 
ficiency of 89.5 per cent. Although the prototype runner has a 
fixed-blade tilt angle, the model runner was constructed so that 
the blades were movable but could be fixed at various tilts and 
complete tests run at the tilt selected. Fig. 12 is a curve of model 
unit horsepower at various blade tilts plotted against efficiency. 
This curve shows that the 25-deg blade tilt produced optimum 
efficiency. It is apparent that a flatter or steeper tilt would result 
in lower efficiency values. 

Therefore, the 25-deg blade tilt was selected for the prototype 
runner together with a runner diameter of 240 in. which would 
produce in the field 71,000 hp at 81 ft net head at point of best 
efficiency. On the model tests this point represents a unit model 
horsepower of 0.262 at optimum efficiency with a 25-deg blade 
tilt, Fig. 12. 

The guaranteed full-load power of 79,000 hp at 81 ft net head 
is represented on the model at a unit model horsepower of 0.292 
with a 25-deg blade tilt. 

Exhaustive cavitation and runaway tests were performed to 
determine the operating limits of the turbine. Fig. 13 is a curve 


aprit 1960 / 101 





of expected performance of the turbines with output plotted 
against efficiency based on model tests. 

Due to the fact that the powerhouse forms part of the dam 
structure, the water intakes are short and of the semispiral de- 
sign which results in the water flow being unequally distributed 
across the inlet, making it practically impossible to conduct an 
accurate field acceptance test in accordance with the ASME Test 
Code for Hydraulic Prime Movers. 

The measurement of effective head and horsepower output 
presents no special problems but accurate measurement of the 
large quantities of water involved is practically impossible. Con- 
sequently accurately conducted laboratory testing was important 
as contract guarantees were based on model performance. 

Preliminary power-gate tests indicate a substantial increase in 
field performance over the model results. Fig. 14 is a curve show- 
ing the field performance of the turbines compared to the model 


APRIL 1960 


102 


performance. The close conformity of the shape of the power- 
gate curves indicates that the field performance verifies the model 
test performance. The increase in capacity of the prototype over 
the model is also an indication of an increase in the efficiency of 
the prototype over the efficiency of the model. 

The first St. Lawrence hydraulic turbine was started in July, 
1958. Successive units were placed in operation through 1958 
with the Power Project to be in complete commercial operation 
in 1959. The St. Lawrence Power Dam will be second only to 
Grand Coulee Dam as the largest power project in the world. 
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Influence of Trailing-Edge Geometry 
on Hydraulic-Turbine-Blade Vibration 
Resulting From Vortex Excitation 


A study was made to determine effects of trailing-edge geometry on the vortex-induced 
vibrations of a model blade designed to simulate the conditions at the trailing edge of a 
hydraulic-turbine blade. 


For the type of tratiling-edge flow encountered, characterized 


by a thick boundary layer relative to the blade thickness, the vortex-shedding frequency 
could not be represented by any modification of the Strouhal formula. The amplitude 
of the induced vibrations increased with the strength of a vortex in the von Karman 
vortex street of the wake; one exception was provided by a grooved edge, which ts dts- 
cussed in some detail. For a particular approach velocity, the vortex strength is pri- 
marily a function of the ratio of distance between separation points to boundary-layer 
thickness, the degree of ‘‘shielding’’ between regions of vortex growth, and frequency of 


vortex shedding. 


— MANUFACTURERS are faced with 
the problem of turbine-runner vibration, the excitation force 
being provided by the vortexes which regularly shed from the 
trailing edge of the blade and which arrange themselves in the 
von Karman vortex street downstream. The situation becomes 
serious when the excitation frequency coincides with one of the 
resonant frequencies of the runner. Parmakian and Jacobsen 
[1]! demonstrated that the detrimental resonant vibrations in 
turbines could be eliminated by a modification of the trailing-edge 
geometry. 

Any object in relative motion to a real fluid generates vorticity 
in the boundary layers; i.e., a rotating motion to an infinitesi- 
mally small fluid element. Somewhere along the upper and lower 
surfaces of the object (for two-dimensional flow) the boundary 
layers must separate from the object. Beyond the separation 
points the vorticity in the boundary layers causes them to roll 
up into vortex spirals at regular time intervals, each vortex being 
made up of a section of a boundary layer. The vortexes originate 
alternately from the two separation points at a certain frequency. 
For every vortex which forms, a transverse force is exerted on the 
object; the direction of the force changes with the alternately 
shedding vortexes. 


1 Numbers in brackets designate References at end of paper. 

Contributed by the Hydraulic Prime Movers Committee of the 
Hydraulic Division and presented at the Hydraulic Conference, 
Ann Arbor, Mich., April 13-15, 1959, of Tae American Society 
ro MECHANICAL ENGINEERS. 

Nore: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 
the Society. Manuscript received at ASME Headquarters, November 
24, 1958. Paper No. 59—Hyd-7. 


The frequency of vortex formation for bluff? cylinders has been 
found to follow Strouhal’s experimentally determined relation 


[2) 


fa s— 
d 


where 


V = velocity of approach 

d = largest projected cross-sectional dimension on approach 
flow 

S = Strouhal number, experimentally determined to be 
about 0.18 


Strouhal’s formula is valid for flow across bluff objects, but is not 
applicable to flow across slender objects such as a turbine blade. 
For slender objects Gongwer [3] found that Strouhal’s formula 
had to be modified in the dimension d, which was chosen as the 
sum of the trailing-edge thickness and an experimentally deter- 
mined fraction of the boundary-layer displacement thickness.’ 
His considerations excluded the effects of trailing-edge geometry. 

The magnitudes of vortex-induced forces were investigated by 
Donaldsen [4]. From experimentation on a wide variety of trail- 
ing edges on the same blade, he found that some edges induced 
great vibrations while others induced negligible or no vibration. 


2 A bluff cylinder may be thought of as one which produces a wake 
width about equal to its largest cross-sectional dimension. 


i (U — u)/U}dy, 


where U = free-stream velocity, u = velocity in boundary layer, y = 
perpendicular distance from wall. 


3 Displacement thickness defined as D = 





Nomenclature 


= velocity of approach, fps 
characteristic dimension of object, 
ft 
: sec 
= frequency of vortex shedding, cps 
fd/V, Strouhal number, dimension- 
less 
fluid density, lb sec*/ft* edge, fps 
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= distance between vortex rows, ft } = 


longitudinal vortex spacing, ft 
= vortex strength or circulation, ft?/ 


velocity of vortex centers relative 
to velocity outside wake, fps F = 
free-stream velocity at trailing 


velocity outside wake, fps 

slope of frequency versus velocity 
(U,,) curve, eps/fps 

ratio of actual vortex strength to 
potential vortex strength 

force per unit span, lb per ft 
Subscripts z refer to conditions at the 

location z in the wake. 


aprit 1960 / 103 





The objective of this investigation was to study the effect of 
trailing-edge geometry on the flow about and downstream of 
the trailing edge and to correlate these effects with the frequency 
and magnitude of the vortex-induced forces. In this initial ap- 
proach, only a qualitative correlation was sought; the results 
to be derived would be valuable in establishing the proper pa- 
rameters for a subsequent numerical correlation. 
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Fig. 2 Test section 


Investigation 


Experimental Setup. A closed-loop water tunnel, as shown in 
Fig. 1, The steel frame of the 
test section, having approximate inside dimensions of 4!/2 x 
4'/, X 18'/ in., was fitted with heavy plexiglas walls to permit 
visual observation and photographic recording of the flow. A 


was used for the investigation. 
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photograph of the test section and model blade mounting is shown 
in Fig. 2. A single-suction centrifugal pump, driven by a 25-hp, 
d-e electric motor, was capable of producing water velocities up 
to 60 fps in the test section 

The blade, designed to simulate the conditions near the trailing 
edge of a hydraulic-turbine blade, is sketched in Fig 3. Fig. 4 
illustrates the mounting and waterproofing of the strain gages 
which were used to sense the blade vibration. Vibration meas- 
urements were made with Consolidated Electrodynamics Corpora- 
tion strain-measuring equipment. 

Velocities at two locations were measured. 
velocity near the trailing edge (defined later) was measured with 
a shaft-reinforced '/,-in-OD total-pressure probe and wall 
static-pressure tap. The velocity outside the wake (defined later) 
was measured with a shaft-reinforced '/,.-in-OD pitot tube. A 
water-mercury manometer was used for pressure indications. 

A high-speed motion picture camera, operated at approxi- 
mately was used for the photographic re- 
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3600 frames per sec, 


FLOW 
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Fig. 3 Sketch of test blade 
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Fig. 4 Method of waterproofing strain gages 
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cording of the flow near the trailing edge. The paths of the flow 
were made visible by a proper amount of air bubbles entrapped in 
the water. 

Experimental Procedure. Ten trailing-edge configurations were 
chosen for this investigation. The various shapes are shown in 
Fig. 5. The edges were divided into three edge sequences, A, B, 
and C, as shown, according to the systematic changes occurring 
in each of the sequences. The trailing edge under consideration 
will be referred to by the number designated in Fig. 5. 

For each trailing edge, the following relations were determined: 


1 Vibration (=frequency of vortex shedding) 
versus velocity. 
2 Vibration amplitude versus velocity. 


3 Vortex strength versus velocity. 


frequency 


In addition to these measurements, observations of the ap- 
proximate locations of the separation points were made. 

The free-stream velocity near the trailing edge U,, was chosen 
as the reference velocity. This velocity was measured outside 
the boundary layer of the blade, about three inches upstream of 
the separation points. 

The velocity was varied in the range 10—60 fps in small steps, 
and an oscillograph recording of the blade vibration (as sensed by 
the strain gages) was made at each incremental velocity. Thus 
data were obtained for the frequency and amplitude relations. 
For the trailing edges where only slight or no vibration was re- 
corded on the oscillograph, the frequencies were measured from 
the motion pictures which were recorded at several velocities for 
each edge, by counting the vortexes shed in a certain time interval. 
For the cases where no vortexes were distinguishable, the fre- 
quency of periodic wave motion in the wake was considered 
equivalent to the vortex-shedding frequency. 

Vortex strength calculations were based on von Karman’s cri- 
teria for a stable configuration of two infinitely long rows of 
staggered vortexes in inviscid flow. If a represents the spacing 
between two vortexes in the same row and h represents the dis- 
tance between the rows, then [5], 


h/a = 0.281 
T =2V/2 ua 
where 


I’ = vortex strength = 2mrv, where v is the tangential ve- 
locity at any radius r of a circular vortex 
u = velocity of vortex centers relative to velocity outside the 


wake 
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Fig. 5 Trailing edges selected for study 
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From a number of measurements of vortex spacing, as revealed 
by the high-speed motion pictures, the average ratio h/a was 
calculated 


h/a = 0.29 


(Measurements were made in the region between the second and 
sixth vortex appearing behind the trailing edge.) Because of the 
close agreement between the measured value and equation (1), it 
was assumed that equation (2) was justified for calculations of 
vortex strength. 

Observing that for some location z in the wake 


U 
u; = Wee a, fa, _ Us ( ye) “wi ka 
, Use 


U2 = velocity outside wake at location z 
f frequency of vortex shedding 
+ = longitudinal vortex spacing at location x 


k = slope of frequency versus velocity (U,,) curve 


where 


a 


equation (2) can be written for the vortex strength at z as, 


: U.. 
Us U., — ka, | a, 


Location .c was chosen as three inches downstream of the trailing 
edge. The ratio (U,,/U,,) was found to be 0.93 for all edges and 
velocities, so that 


Tr, = 


/9 


Tr. =2V7 2 


2+/2 U,,{0.93 — ka,]a, (4) 

The longitudinal vortex spacing a was measured from the 
motion pictures; it was found to be essentially invariant in the 
region between the second vortex behind the trailing edge and 
location xz and did not vary with velocity (neglecting conditions 
of resonance). The vortex strength caused by each edge was 
therefore calculated from equation (4) as a linear function of U,, 


for nonresonance. 


Table 1 


Vel, 
fps 


Vel, 
fps 


Freq, 
cps 


Vel, 


Ips 


Edge 
no. 

3 740 
780 
720 
790 
820 
840 
850 
850 
850 
850 


Frec ly 
cps 
350 
389 46.° 
400 7 
450 

460 

560 

510 

515 

550 

590 

655 

608 

667 


Freq, Edge 
eps no. 


145 2 36 
ae 200 38 
21.9 200 40 
25.9 237 42 
28.6 280 44 
30.6 298 46: 
32.9 308 47 
34.8 324 49 
36.3 326 50 
38.3 350 51 
39 .¢ 366 54 
41 357 56 
42 375 
44 394 
46 384 
aT a 436 
432 
500 
532 
580 


Edge 


ho. 


] 16.5 


>= OIW OW OD 


250* 
370* 
190* 
520* 
633* 
700* 


795* 


Gt ho 00 00 


259 
283 
330 
284 
340 
323 
400 


590* 
738* 


950* 


1? oe ota | 
“I 0 OC & Ow 


465 
497 
582 
596 
667 


250 35 .é 

250 36 

280 39 
32.6 300 41.1 
34 325 43.2 


100 
111 
196 


* Measured from motion pictures. 
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Results 


In general, the relationship between frequency of vibration, or 
frequency of vortex shedding, and velocity followed the ex- 
pected linear trend. It should be stated here that this was true 
only for the rigid (nonresonant) blade condition. During and 
near resonance, the elastic properties of the blade affected the 
rate of vortex shedding; an example of this effect will be shown 
later. 

The frequency-velocity data are listed in Table 1. Graphs of 
frequency versus velocity for all trailing edges are shown in Fig. 6. 
The linear frequency versus velocity curves are not completely 
valid for edges Nos. 2 and 3, which exhibited increasing slope 
trends in the upper velocity ranges. In these cases the slopes 
formed by experimental points in the lower ranges (0-40 fps) 
were used for Fig. 6. 

For several edges during the first series of test runs one or more 
screws holding the blade in the test section came loose and even 
fatigue failed when the blade was resonating moderately to 
violently. Thus, the rigidity of the blade mounting was altered, 
which in turn changed the resonant frequencies. In order that 
the vibration amplitudes caused by the various edges could be 
compared, it was necessary that the comparison be made at the 
same resonant frequency for all edges. Because the requirement 
of a constant blade-mounting rigidity was not satisfied in the first 
series of tests, a second series of tests was run, incorporating a 
modified blade mounting. The ten mounting screws were re- 
placed by six tapered dowels, which were driven through the 
mounting plates into the top and bottom of the test blade. In 
the interest of not disturbing the blade mounting once the blade 
had been installed in the test section, interchangeable trailing 
edges, which could be fitted to the blade by screws, were ma- 
chined. 

The second series of tests gave satisfactory amplitude data, ob- 
tained at a resonant frequency which varied between the extremes 
of 245 and 268 eps for the two most disagreeing tests. (Most 
likely this frequency spread was caused by a varying degree of 
fluid damping rather than a varying mounting rigidity.) A 
sample of frequency and amplitude data obtained in this test 


Tabulation of frequency versus velocity data (recorded in first test series) 


Edge 
no. 


Freq, 
cps fps 


208 8 37.8 
255 39. 
256 41 
257 42 
258 44 
258 55 
269 
273 
265 
369 
357 
386 
385 
470 


Freq, Edge Vel, 

fps eps no. fps 

6 22.9 217 7 33 
25.5 230 33 
27 233 36 

236 38 

247 ‘ 

277 

283 

286 

293 

310 


Edge _ Vel, 


no. 


Vel, 


CO WKNQO wm 


103 
102 
106 
112 
114 
115 


118 
394 
404 
390* 
400 


> g + 
PO WNHKHOUERPRONNWE 
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series appears in Fig. 7 (for edge No. 1); frequencies measured 
in the first test series are included to demonstrate agreement be- 
tween the two tests. 

Refer to the frequency points of second test series in Fig. 7. 
If a point-to-point connection is made in the velocity range 22-52 
fps, an S-shaped curve results. Similar more and less pronounced 
S-curves occurred for the other edges in the second test series 
and were concluded to be caused by the particular elastic proper- 


1000, T T T T T 
NOTE: 1-DATA POINTS OBTAINED DURING 
RESONANT VIBRATION NEGLECTED 
IN DETERMINATION OF 
+ CURVES. + + 


2-DATA FOR CURVES 
| LISTED IN TABLE | 
| | 
t + 
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FREQUENCY ~-CPS. 


50 60 
VELOCITY -FPS 


Fig. 6 Frequency of vibration versus velocity for all trailing edges, 
recorded in first test series 
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Fig. 7 Frequency and relative amplitude of vibration versus velocity 
for edge No. 1 


ties of the vibratory system. The first test series, incorporating 
a different blade mounting, gave deviations of frequency points 
from the straight-line curve to a smaller extent, thus rendering it 
more reliable for the frequency versus velocity determinations 
during nonresonance operation. 

Table 2 lists the results of the (a) frequency determinations 
(as the slopes of the frequency versus velocity curves); (6) ampli- 
tude determinations (as relative amplitude of vibration at the 
resonant frequency near 260 cps, referred to the amplitude = 100 
caused by edge No. 1); (c) vortex strength calculations made 
from equation (4), valid for nonresonance operation; (d) ob- 


EDGE No.7 


EDGE No.8 


te Me 
EDGE No.9 


EDGE No.5 EDGE No lO 


Fig.8 Appearance of vortex streets as produced by the various trailing 
edges at the 260-cps r t freq y, recorded in second test series 





Table 2 Summary of results 


Frequency 
vs. velocity, 
cps/fps 
9.2 
10.3 
12.1 


Trailing edge 
Sequence no. 
A 


_ 
w 
~I 


-_— 


CWONOKorwhde 
— 
COmSIDOS 
CIN AOwWNnothds-~] 


C 1 
10 


° From first test series, for nonresonance. 
+ From second test series, for resonance at about 260 cps. 


Journal of Engineering for Power 


Distance between separation points 


Equal to blade thickness 100 
( 
Decreasing 


Equal to blade thickness 100 
Equal to blade thickness 190 
Equal to blade thickness 

Decreasing 43 


Equal to blade thickness 
Equal to blade thickness 31 


Relative? 
amplitude of 
vibration 


Vortex® strength, [t?/sec 


35 & 1078U ce 
38 26 X 10-3U oe 
3 7X 107U oe 
No vortexes visible 
No vortexes visible 
35 X 10-*U., 
43 X 10-3U., 
73 X 1073U 4 
25 X 10-3U 4: 
No vortexes visible 
35 XK 10-2Uo: 
35 X 10-3U. 


No vibration 
No vibration 


380 


No vibration 
100 
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served trend of distance between separation points in going from 
one edge to the next. 

Some frames from the motion pictures have been selected, 
Fig. 8, to show the appearance of the vortex streets at the 260-cps 
resonant frequency. 


Discussion 


Frequency Versus Velocity. The great variations in the slopes 
of the frequency versus velocity curves may seem surprising in 
view of the results published earlier. Gongwer’s modified 
Strouhal formula‘ [3] predicts identical frequency versus velocity 
relationships for all edges. Gongwer experimented with dif- 
ferent trailing-edge thicknesses while the trailing-edge geometry 
remained the same, semicircular. Thus an important parame- 
ter remained constant, the ratio of trailing-edge thickness to dis- 
tance between separation points. The present study shows that 
this parameter greatly affects the frequency of vortex shedding, 
which is in agreement with basic considerations. 

The data have been analyzed for a possibility of linking the 
vortex-shedding frequency with a characteristic dimension of the 
flow over the trailing edge into a further modified Strouhal rela- 
tion. It was soon realized, however, that such a dimension would 
have to be different for the flow over edges Nos. 1, 6, 7, and 10, 
which all produced different vortex-shedding frequencies. All four 
edges appeared to produce equal distances between the separa- 
tion points, and their boundary layers were assumed to be identi- 
cal up to the separation points. Thus the desired dimension 
could not be derived from the locations of the separation points 
or the characteristics of the boundary layers. It would either 
have to be determined from other features of the trailing-edge 
flow, or the assumptions made in the foregoing should be checked 
with careful measurements. Therefore it was felt that the data 
available were not sufficiently complete and accurate for the de- 
termination of a modified Strouhal relation. 

Correlation of Vibration Characteristics. The summary of results, 
as listed in Table 2, serves as a basis for the following discussion. 

Going through edge sequences A and B, it is apparent that the 
amplitude of vibration during resonance increases with increasing 
vortex strength (during nonresonance), a fact which may not 
seem surprising. (U,, varied only slightly between the various 
edges at resonance.) It is important to establish this relationship, 
however. Edge sequence C seems to violate the foregoing con- 
clusion; while the vortex strength for edge No. 10 is the same as 
that for edge No. 1, the vibration amplitude caused by edge No. 
10 is less than one third of the amplitude caused by edge No. 1. 
Further discussion on this behavior will be given later. 

It is believed that the following factors affect the strength of 
the vortexes forming near the trailing edge: 


1 Overlap of velocity fields of oppositely shedding vortexes, 
thereby allowing the vortexes to destroy each other partly or 
completely. 

2 Intensity of turbulence in the region occupied by the vor- 
texes while “rolling up.” 

3 Viscous effects. 

4 Frequency of vortex shedding. 


Of these four factors, only Nos. 1 and 4 are believed to be of 
major importance; the other factors will be neglected in this 
discussion. 

4 According to Gongwer, the frequency of vortex formation should 
be represented by 

S 
f 6 
where S = Strouhal number = 0.18; 6 = d + 6’; d = trailing-edge 
thickness; 6’ = 0.0297[z/(V2/v)'/s]; 2 = blade length; v = kine- 
matic viscosity. 
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Overlap of velocity fields is a common occurrence in flow across 
slender objects; it may at times destroy or prevent all vortex 
formation. (In the case of conventional airfoils, no vortexes were 
observed by Chuan and Magnus [6] for angles of attack less than 
six deg.) A criterion for the degree of vortex-strength attenuation 
caused by this factor seems to be the ratio of boundary-layer 
thickness at the separation points to the distance (normal to the 
approach flow) between the separation points. The greater value 
this ratio assumes, the more the vortex-velocity fields overlap, 
and the smaller the fraction of the potential vortex strength which 
survives. It is conceivable that this effect could be diminished, if 
desired, by introducing a shield between the two regions behind 
the two separation points where the vortexes mature, provided the 
original locations of the separation points were retained. The 
effect. of shielding appears to be of importance in the present 
study. 

If the quantity € represents the ratio between the actual vortex 
strength and the potential vortex strength (which would exist if 
no viscous, turbulent, and vortex-velocity-field overlap effects 
were present) then it can be shown that [7] 


where I’, = strength of vortex at location z in wake and € is 
evaluated for location z. Thus, by assuming that € and U,, are 
constant, vortex strength is inversely proportional to vortex- 
shedding frequency. 

The following explanation can now be offered to account for the 
behavior of the vortex-induced excitation forces. 

First, consider edge sequence A. Going from edge No. 1 
through edges Nos. 2, 3, 4, 5, the separation points approach 
each other. One separation point is always located at the lower 
corner of the edge while the other stepwise approaches the 
lower separation point. This effects an increasing degree of 
vortex-velocity-field overlap, which should favor a stepwise re- 
duction in vortex strength. The frequency of vortex shedding 
increases, which also favors a decrease in vortex strength. The 
total effect of these factors is a diminishing vortex strength and 
excitation force through ‘the sequence. 

Next, consider sequence B. For the first three edges of the 
sequence (Nos. 1, 6, 7) the distance between the separation points 
at the sharp corners is the same. For the remaining edges of the 
sequence (Nos. 8 and 9) the separation points move increasingly 
farther down the beveled surfaces and therefore approach each 
other. <A certain degree of vortex-velocity-field overlap is 
present for edge No. 1, but the two regions occupied by the op- 
positely shedding vortexes become increasingly more isolated by 
the introduction of the angular “shields’’ of edges Nos. 6 and 7. 
This reduces the velocity-field overlap, and the vortex strength 
increases in the edge order Nos. 1-6-7. Simultaneously the fre- 
quency of vortex shedding is reduced which further favors an 
increasing vortex strength and induced excitation force. The two 
remaining edges (Nos. 8 and 9) of the sequence produced separa- 
tion points successively closer to the vertex of the triangular edge, 
thus diminishing the distance between separation points. The 
shields became increasingly less effective, the degree of vortex- 
velocity-field overlap increased, and the shedding frequency in- 
The net effect was a decrease in vortex strength and 


’ 


creased. 
excitation force. 

Edge sequence C proved very interesting. The vortex strength 
for the grooved edge, edge No. 10, was the same as for edge No. 1, 
yet its vibration amplitude was about one third of the amplitude 
“aused by edge No. 1. Because this behavior could not be ex- 
plained on the basis of previous discussion, it was suspected that 
a new flow phenomenon was associated with the grooved edge. 
A study of the motion pictures of edge No. 10 revealed the pres- 
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Fig. 9 Flow in groove of edge No. 10 


ence of a stationary slug of fluid, bounded by the edge walls and 
the shedding vortexes, which was in rotational motion and was 
oscillating with the frequency of vortex shedding. The following 
hypothesis may now be justified. Refer to Fig. 9. 

The oscillation of the fluid in the groove was activated by the 
vortexes shedding from the two separation points. The sense and 
magnitude of the oscillatory motion would at any instant be a 
function of the magnitude and direction of the net velocity in- 
duced at the trailing edge by the vortexes in the process of rolling 
up. Because the velocity induced by vortexes in one row is 
opposite in direction and phase to the velocity induced by the 
other row, and the magnitude of the induced velocity falls off 
rapidly as a vortex recedes from the point of vortex maturity, a 
periodically reversing induced velocity forces the fluid slug into 
oscillation. 

It is generally recognized [3] that the excitation forces are 
derived from the constantly changing “bound’’ circulation in- 
duced about the blade by the shedding vortexes, through the 
Kutta-Joukowsky relation: 


Ver 
where 


F = force per unit span 
V= 
T = 


approach velocity 
bound circulation 


The bound circulation is ideally always equal to the net circula- 
tion (or strength) of the shedding vortexes in magnitude, but has 
opposite sense of rotation. 

The oscillating fluid slug represents a periodically changing 
circulation which will very nearly be in opposite phase to the 
periodically changing net circulation of the shedding vortexes 
Therefore, the net ciréulation behind the trailing edge will be 
reduced by an amount determined by the magnitude of the fluid- 
slug circulation. This in turn requires a reduced bound circula- 
tion, and the excitation force is reduced. 


Remarks 

The importance of the boundary-layer characteristics at sepa- 
ation and their effect on vortex shedding has been pointed out 
For the particular conditions of this investigation, the boundary- 
layer thickness was approximately equal to the blade thickness. 
A calculation made for a representative hydraulic-turbine blade 
shows the same order of magnitude for the ratio of boundary-layer 
thickness to blade thickness at the trailing edge. It appears that 
the greater thickness the boundary layer has relative to other 
dimensions in the trailing-edge flow, the more it influences shed- 
ding frequency and vortex strength. 


Conclusions 


1 The frequency of vortex shedding caused by the various 
trailing edges cannot be determined accurately from the classical 
Strouhal formula or any known modification thereof. Most of 
the edges investigated yielded linear frequency versus velocity 
relationships with slopes ranging from 7.6 to 16.7 eps per fps. 

2 Trailing edges yielding the greatest vortex strengths during 
nonresonance induced the greatest vibration amplitudes during 


Journal of Engineering for Power 


resonance and, therefore, also greatest excitation forces. One 
important exception was provided by the grooved edge which 
produced relatively great vortex strength but only slight vibra- 
tion; the explanation for this is found in the rotationally oscillat- 
ing slug of ‘“‘dead’’ water in the groove which produced a circula- 
tion in opposite phase to the net circulation of the shedding 
vortexes. 

3 In general, the criteria determining the vortex strength at 
the point of vortex maturity are: 


(a) The ratio of distance between separation points to bound- 
ary-layer thickness. The greater this ratio, the greater the 
strength of the mature vortex. 

(b) The solid shielding between regions of vortex growth. 
Thus an extension of the blade proper (the shield) beyond the 
separation points increases the strength of the mature vortex. 

(c) The frequency of vortex shedding. The lower the frequency 
the greater the vortex strength. 

(d) Turbulence and viscosity. 


The foregoing conclusions are based on velocity as a fixed 
parameter. 
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DISCUSSION 
P. S. Eagleson’ and G. H. Toebes® 


The authors are to be commended for shedding another ray of 
light on what is quite a complex hydroelastic problem. 

Recent unpublished experiments of a similar nature have been 
made at M.I.T. A few observations from these tests may help 
in interpreting the authors’ findings. 

As the free-stream velocity is increased gradually from zero the 
spectrum of plate vibrations (as determined by an harmonic ana- 
lyser) grows broader. This growth takes place at the higher 
frequency (right-hand) end of the spectrum through response of 
the plate to the current (for that particular free-stream velocity) 
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“Strouhal”’ frequency of eddy shedding. As the free stream 
velocity is increased still more the upper limit of the spectrum is 
of course raised but the spectral components which originated at 
lower velocities do not vanish. In the absence of structural 
resonance they gradually decrease in intensity until they dis- 
appear in the noise at the left-hand end of the spectrum. 

In short, a given frequency component does not appear in the 
spectrum of plate vibrations (and hence presumably in the wake) 
until it has first appeared as the eddy shedding frequency. 

This broad wake spectrum permits simultaneous resonant 
plate excitation in torsional as well as flexural modes. 

It is unfortunate that the authors did not present the fre- 
quency-velocity data of either the first or second test series in 
graphical form. The straight lines of Fig. 6 may obscure what 
the writers have found to be quite significant aspects of the plate 
behavior. 

The M.I.T. experiments showed the character of the near- 
resonant plate oscillations to become self-excited for certain 
trailing edge geometries. This can also be observed from the 
breadth of the constant frequency resonant ‘‘steps’’ in the au- 
thors’ data particularly for plates No. 7 and 8. 

It would be of interest to examine the spectrum of vibrations 
for the authors’ plate No. 7. With two natural frequencies so 
close together (i.e., 200 and 260 eps) it is likely either that two 
predominant vibratory frequencies were present or that vibra- 


tion alternated between the two values. This latter behavior has 
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been observed by the writers for plates with the same trailing 


edge configuration as plate No. 7. 


Authors’ Closure 

The observation made by Eagleson and Toebes on the behavior 
of the spectrum of blade vibrations for increasing free-stream 
velocities is of great interest and introduces speculation on the 
possibility of blade resonance at frequencies other than the vortex- 
shedding frequencies. In no case of the authors’ experiments 
was the major (in amplitude) frequency component of blade vi- 
bration different from the vortex-shedding frequency, except for 
low free-stream velocities where extraneous frequency compo- 
nents (noise) dominated. 

In response to the last paragraph of the discussion the following 
comments are made. As the water velocity (for edge No. 7) was 
increased beyond the value at which the first resonant frequency 
(200 cps) was reached, the vibration amplitude would first in- 
crease, reach a maximum, then decrease until the velocity 32 fps 
was reached. At this velocity the second resonant frequency 
(260 cps) began interrupting the first resonant frequency for short 
time intervals. Slight increases in velocity caused the frequency 
to alter between the two resonant frequencies, the second be- 
coming increasingly more predominant. Beyond a velocity of 
34 fps the second resonant frequency only would be present, and 
the amplitude would increase, reach a maximum, then gradually 
decrease until resonance ceased at 43 fps. 
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Water-Hammer Damage to 


C. C. BONIN 


Engineering Manager, 
Ebasco Services Incorporated, 
New York, N. Y. 


Oigawa Power Station 


In June, 1950, errors in operation and malfunctioning of equipment resulted in water- 


hammer surges which burst one 9-ft-diam penstock of the three-unit 68,200-kw Oigawa 
Power Station in Japan. The accident resulted in the death of three workmen and about 


one half million dollars in damage to equipment and powerhouse. 


This paper reviews 


the events leading to the penstock failure and the resulting damage. 


The Oigawa Power Station 


Tux O1gaAwa PowER SraTIon is a hydroelectric in- 
stallation located on the Oi River in Shizuoka Prefecture, about 
100 miles southwest of Tokyo, Japan. The power plant was con- 
structed in 1936 and was operated in 1950 by the Japan Power 
Transmission and Generation Company, Inc. Three Francis 
turbines operate at a head of 403 ft and develop 32,000 hp maxi- 
mum each. In addition, there are two smaller turbines, one 
driving a 600-kva house generator and the other coupled to a 
stand-by oil pump. Three steel penstocks approximately 9 ft in 
diameter carry water 800 ft from a differential surge tank to the 
powerhouse. 

Flow through the penstocks is controlled by the turbine wicket 
gates for normal operation and by butterfly valves for main- 
tenance shutdowns. Head gates are provided at the forebay in- 
take for complete shutdown of the penstocks. A pressure regu- 
lator for each unit is connected between the wicket gates and 
butterfly valve and serves as a water bypass valve to compensate 
for sudden load changes on the turbines. It is operated through 
the turbine-control mechanism. The two small turbines, whose 
failure to operate contributed substantially to the accident de- 
scribed here, can be supplied from any of the three penstocks 
through appropriate valves. 

Fig. 1 is a schematic diagram of the plant hydraulic features 
described in the foregoing. Fig. 2 shows a profile of the penstocks. 
Fig. 3 is a cross section through the powerhouse. 


Failure of Penstock 

On June 18, 1950, about one half hour after the start of load- 
reduction operations preceding a scheduled maintenance shut- 
down of the entire plant, the butterfly valve in No. 3 penstock 
slammed shut due to a series of operating errors and mechanical 
failures. The resulting water-hammer surge upstream of the 
valve caused a section of steel penstock about 25 ft in length to 
split along its crown and lay out nearly flat against the hillside, 
as shown in Figs. 4 and 5. Subsequent excessive discharge from 
the broken penstock created a vacuum upstream aud about 175 
ft of penstock collapsed, resulting in the distorted shapes shown in 
Figs. 6, 7, 8, and 9. The penstock profile, Fig. 2, indicates the 
location and extent of damaged pipe. 

For about two hours, until the head gate could be closed, water 
poured from the penstock breech through and around the power- 
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house, washing out about 16,000 cu yd of earth. Three workmen 
were killed and one seriously injured as they tried to escape from 
the butterfly valve chamber. Turbine, generator, and control 
rooms were flooded and partially buried by about 1600 cu yd of 
sediment. Figs. 10 and 11 are photographs showing wrecked 
auxiliary buildings and the sediment-filled tailrace. Fig. 12 
shows the smoking No. 3 transformer which was burned out. 

The approximate cost of repair was $445,000 of which $280,000 
was for civil engineering work and $165,000 for electrical and 
mechanical repair and replacement. The total estimated power 
loss was 90,000,000 kwhr. 


Events Leading to Penstock Failure 


Sand accumulation in the forebay over a period of six months 
had been such that a one-day plant shutdown was requested and 
approved. Repairs on runner No. 1 and minor maintenance work 
were planned to coincide with the sand-clearance work. Complete 
draining of the lower end of penstock No. 1 was therefore neces- 
sary. 

Just prior to shutdown all three units were operating at full 
load and the station-use generator was in operation, taking its 
flow from penstock No. 1. The small turbine coupled to the 
stand-by oil pump also was connected to penstock No. 1. The 
plant-shutdown procedure called for transferring the station-use 
load first to Nos. 2 and 3 generators, then to an outside loop 
connection. At midnight, according to the prearranged schedule, 
unit No. 1 and the station-use generator were shut down, the No. 1 
penstock butterfly valve was closed, and draining of the penstock 
was begun through the pressure regulator. No thought was given 
at the time to transferring the station-use and oil-pump turbines 
to the No. 2 or No. 3 penstocks. 

A few minutes after midnight the electrical loop connection be- 
tween the Oigawa Station and another local power plant was 
made as planned and load was dropped from both unit No. 2 and 
unit No. 3 generators. The turbines were permitted to operate 
at no load since there were not enough operators immediately 
available to secure them. 

Operating problems now began to compound. 
voltage could not be maintained by the outside loop and dropped 
from 3400 to 2800 volts, tripping undervoltage relays which shut 
down the oil pumps supplying oil and pressure to the operating 
cylinders of governors, wicket gates, butterfly valves, and pressure 
regulators. Approximately 22 min had elapsed since the plant- 
shutdown operation had begun. Efforts were immediately made 
to put the station-use generator back on the line, but unit No. 1 
penstock downstream of its butterfly valve had substantially 
drained and no water was available. The stand-by oil-pump tur- 
bine apparently was designed to start automatically, but it, too, 
had no water supply. Efforts to open the valve connection be- 
tween penstocks No. 1 and No. 2 were unsuccessful because of un- 
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balanced water pressure on the two sides of the valve and an ac- 
cumulation of silt. By this time the operating personnel were 
confused and disorganized. 

The switchboard attendants threw the units Nos. 2 
bine master switches to the “stop” position but the turbine-room 


and 3 tur- 


crew wedged the turbine shutdown solenoids in place and reset 
controls to the starting position. Oil supply and pressure were 
rapidly being depleted in the governors and valve controls, and 
a considerable amount of air had mixed with the oil in the operat- 
ing evlinders. As model and prototype tests later verified, the No. 
3 turbine governor started “hunting” which completely opened 
Since the wicket gates no longer had 

moved to a hydraulically 
Water was thus flowing un- 


the pressure regulator. 
sufficient oil to control them, they 
balanced partially open position. 

controlled through the wide-open pressure-regulating valve and 


partially open wicket gates. At this point the butterfly valve, 
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slightly closed and with po oil in its operating cylinder, was 
slammed shut by the flowing water. Normally, this valve could 
not be closed in less than 28 sec as long as its operating cylinder 
was full of oil. Closure now occurred in an estimated 0.3 sec. 
A typical and severe valve-closure surge developed upstream and 
the penstock failed. Since the communication facilities between 
the powerhouse and the forebay were among the immediate 
casualties, the forebay gates were not closed for nearly two hours, 
during which time the rampaging water took three lives and in- 
flicted the severe damage described. 

Fortunately the timing between shutdown operations for units 
Nos. 2 and 3 was such that unit No. 2 turbine governor did not 
tend to hunt, the unit No. 2 pressure regulator did not open, and 
the final butterfly-valve closure was from the nearly closed posi- 
tion apparently reached while oil pressure was still available. 
The final closure surges were therefore insignificant. 
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Fig. 4 Water hammer caused about 25 ft of penstock to split 
and flatten out 


Fig. 6 Nearly 175 ft of p 


Fig. 5 General view of burst penstock 


d as shown in this view 





Water-Hammer Conditions 


As noted previously, calculations have indicated that the 
butterfly-valve closure may have occurred in as little as 0.3 sec. 
The period of surge-wave propagation, u, was determined to be 
0.49 sec where uw = 2L/a and a = 3280 fps, an experimental 
value obtained from tests on the unit No. 1 penstock. For pur- 
poses of water-hammer computation the valve closure may be con- 
sidered instantaneous since the closure time was less than the 
period of wave propagation. From Joukovsky’s law 
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Hp = — Vo 
g 


where 
Hp = head rise over static head 
a = 3280 fps 
Vo = velocity of water prior to valve closure 
The velocity Vs was determined by calculation to have been 
between 10.7 and 13.5 fps as a result of flow through the wide- 
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Fig. 7 Close-up of section of ¢ 
upstream of break 





Fig.9 Another view of collapsed penstock 
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Fig. 8 View of distortion resulting from failure 


open pressure regulator and the partially open wicket gates. The 
head rise was therefore between 1100 and 1380 ft and the total 
pressure at the upstream side of the valve, with static head in- 
cluded was between 1500 and 1780 ft. The penstock was 
originally designed for a maximum pressure of 611 ft. 

After the penstock had burst, water discharged down the steeper 
grade below expansion joint No. 2 faster than it could be sup- 
plied from the gentler grade upstream. Pressures less than at- 
mospheric developed and the penstock collapsed. 


Steel Failure 

Japanese engineers made a thorough study of the failure of the 
steel to determine where the break first occurred and whether 
deficiencies in the steel or in fabrication might have been at fault. 
There was also concern that penstocks Nos. 1 and 2 might not be 
safe. However, the investigation revealed no evidence of un- 
usual weakness in the penstock. The steel was originally de- 
signed for a working stress of 18,000 psi, with yield point 58,000 
psi. After extensive testing it was concluded that the steel could 
be stressed to at least 42,500 psi before failure. 

The penstock is made up in 25-ft sections joined together by 
straps with a single row of rivets holding the strap to each pen- 
stock section. The sections are constructed of welded steel plates 
staggered in such a manner that the welded joints are not in line. 
The fabrication detail can be seen in Fig. 4. Plate thickness at 
the ruptured section was about 0.83 in. There were no stiffeners. 
The sections of penstock which collapsed under negative pressure 
had plate thicknesses varying from 0.35 to 0.55 in. with stiffeners 
on sections of 0.47 in. and under. 


aprit 1960 / 115 





Fig. 12 Transformer No. 3 burning out 


At the impact of the surge the steel started to tear at the down- 
stream top edge of section No. 31 and ripped open along an almost 
straight line to the upstream edge of the section as shown in Fig. 4. 
About half the break was along welded seams and the other half 
through the center of plates making up the sections. 


Conclusion 

It is apparent from this brief review that the primary cause of 
the accident was the operator's lack of knowledge and judgment. 
The accident report of the Special Mechanical Committee of In- 
vestigation gave as its first two recommendations the improve- 
ment of technical training for station employees and the clarifica- 
There 
were, of course, a number of recommendations for mechanical and 


tion of normal and emergency operating instructions. 


electrical changes and for proper maintenance. 

This accident is a stern reminder that destructive water ham- 
mer is not just a theoretical concept, but a very real problem in 
the design and operation of water-control facilities. 
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Fig. 11 Debris and earth in tailrace 


DISCUSSION 
G. S. Bingham? 


In March, 1955, errors in operation resulted in surges which 
damaged the intake gate controi equipment and gate house at the 
Ladhon Hydroelectric Project in Greece. The amount of damage 
was small but the incident required a plant shutdown of more than 
a week for complete inspection of the gate and power tunnel. 
This discussion reviews the events leading to the occurrence. 

The writer’s firm was engaged, from 1950 to 1955, to manage 
and direct the setting up of the Public Power Corporation, a 
complete, integrated public utility enterprise covering a large 
part of the mainland of the Kingdom of Greece. Constructing 
and operating a project such as this meant working with many 
different nationalities and a babel of tongues. Drawings were 
in German, American, and Greek; Italian, American, and Greek; 
or French, American, and Greek. In several instances British 
had to be translated into American and Greek—technical British 
being different from technical American in practically every term. 

Although Greek is the foundation for many technical terms, 
modern-day Greek has few equivalents for many of the words 
we find commonplace. As an example, the tachometers for the 
generating units, in a version the station operators could under- 
stand, were labelled ‘“strofometron-strovilov,’’ which literally 
means “‘the meter of revolutions of a rotating body.”’ Naturally, 
there were misunderstandings which resulted in some unusual 
occurrences, one of which will be described here. 

The Ladhon Hydroelectric Project, in the heart of the Pelo- 
ponnesus about 8 hours drive from Athens, was one of three 
hydros constructed. The engineers were Italian, the workmen 
mixed Italian and Greek, the supervision American and Italian. 
The project consists essentially of a buttress-type dam 115 feet 
high, a power tunnel 13 feet in diameter and over 5 miles long 
with tunnelled surge tank, a steep penstock, and a powerhouse 
with two 25,000-kw hydro units which operate under a gross head 
of 750 feet. The tunnel and penstock profile is shown on Fig. 
13. There is also a short tunnel 14 feet in diameter around the 
dam, which was used for river diversion during construction and 
is now a reservoir bottom outlet. This is equipped with two 
vertical slide gates in tandem, the upstream one normally fully 
open and used only as a spare, the downstream normally closed. 

At the time of the incident the powerhouse had been in opera- 
tion for several weeks and the reservoir was full. In the morning 
the downstream diversion tunnel gate was opened about 2 feet to 
reduce the reservoir level somewhat for witnessing tests on the 
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spillway gates. During the noon lunch hour it was decided to 
reduce the discharge through the diversion tunnel because of 
known difficulties with scour downstream from the gates. In- 
structions were given via field telephone to the assistant gate 
tender at the dam to close the diversion tunnel gate to a one- 
foot opening. Apparently, because of language difficulties, the 
assistant gate tender did not lower the diversion tunnel gate but 
proceeded instead to almost close the intake gate. Having some 
misgivings about this, he then tried to call the construction 
office for confirmation but instead reached another part of the 
project, where he was understood to have closed the diversion 
tunnel gate and was told to do nothing further. 

‘The intake gate, with a full reservoir, was under a considerable 
head, so with the one-foot opening that had been left it was 
passing about two thirds of the water required for the one unit 
then in service. Of course the tunnel immediately started to 
drain and the unit head to fall off, rapidly at first as the gate 
shaft drained, then very gradually as the level dropped in the 
sloping power tunnel. 

At the powerhouse the Chief Operator was away at lunch, as 
was the American Plant Superintendent. The Assistant Operator 
first noted that the unit output had fallen when he took his hourly 
instrument readings. He called the Dispatcher in Athens to ask 
if possibly some other units on the system had picked up more 
load and backed the Ladhon unit off but learned this was not the 
case. Soon after, the man on the turbine floor advised the oper- 
ator that there had been a large reduction in penstock pressure 
and that it was slowly continuing to drop. The Plant Super- 
intendent, returning a short time later from lunch, called the 
dam and learned what the trouble was. After a hurried discus- 
sion with the Italian start-up man, it was decided the intake gate 
should be re-opened and the dam gate tender was so advised. 
The unit had then been in service somewhat less than one and 
one-half hours from the time the gate was lowered. It was con- 
tinued in service. 

Flow calculations and the pressure readings taken indicate 
that the tunnel was about half empty at the time the gate started 
to re-open. From Fig. 14 it will be seen that there is a combined 
air inlet and access shaft about 4 feet in diameter leading to the 
intake house, where it is closed off by a steel plate in the floor. 
A supplementary opening about 18 inches in diameter is connected 
to this shaft and terminates outside the intake house, making two 
right angle bends in the process. Soon after the gate started to 
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open, a series of rumbling noises were heard, coupled with blasts 
of air and some water coming from the access shaft. The steel 
plate cover was removed to relieve the air pressure, whereupon 
the water wet the gate hoist motor. Fearing damage to the mo- 
tor, the operator stopped opening the gate and removed the 
motor fuses. Soon the air blasts became so strong that he ordered 
all workmen in the vicinity to flee and ran away himself, leaving 
one half of the large gate house door open. Shortly thereafter 
the action became more violent, blowing out most of the windows 
in the building. The other half of the large door was then blown 
bodily into the lake and the grand finale saw the lean-to type roof 
of the intake house lift and slide completely off, coming to rest 
leaning against the front of the building. Great spouts of water 
then rose as much as 200 feet in the air, if some eyewitness re- 
ports are to be believed. 

The whole incident took about 10 minutes, which by calcula- 
tion was about the length of time required to refill the tunnel. 
The powerhouse crew was told what had happened and the unit 
shut down. 

The gate-house damage also included shattering of a part of the 
reinforced concrete floor, which fell into the gate shaft. A num- 
ber of timber stop logs at the top of the shaft and the steel hatch 
covers also were dislodged and dropped down the shaft. 

An attempt was made to close the gate by hand but it could 
not be lowered all the way because of some obstruction. A call 
for divers was sent to Athens and a team secured with some dif- 
ficulty, since the depth was well beyond their usual experience. 
They found that some of the concrete and other debris prevented 
gate closure but it was soon removed, the gate closed, and the 
tunnel drained, 

The Ladhon tunnel was constructed through much limestone 
rock of questionable character and a considerable length through 
clay, requiring very heavy timbering and a very thick concrete 
lining. The engineers had recommended reinforcing a large part 
of its length but it was finally decided that the thick lining, 
coupled with a thorough grouting program, would obtain a safe, 
sound construction. After the lining and grouting were com- 
pleted the whole tunnel was pressure tested, in sections, to the 
maximum expected surge pressure and was found satisfactory. 
The incident just described, however, apparently had subjected 
the tunnel to pressures in excess of the test conditions and it was 
considered advisable to make a thorough inspection. No evi- 
dence of any damage was found. 
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Fig. 14 Public Power Corporation, Ladhon Hydroelectric Project, intake 
gate shaft and house 


The tunnel was then refilled and the plant put back in service. 
The total shutdown was of about 12-days’ duration. 


C. C. Crawford? 


Three persons were killed; another seriously injured. Re- 
pairs cost one-half million dollars. Service was interrupted. 
Prestige, confidence, and goodwill were 
Destructive water hammer 

Since we learn from our 


Power revenue was lost. 
lost. 
is indeed a very real design problem. 
failures, Mr. Bonin has performed a valuable service by publica- 


What a tragic series of events. 


tion of this article. 

An important cause of the failure at this Power Station was a 
design in which a central hydraulic system was used to control 
and actuate wicket gate, pressure regulator, and butterfly valve. 
In the powerplants of the United States Bureau of Reclamation, 
each of these water control mechanisms has an independent hy- 
draulic system. This essentially precludes the type of accident 
which occurred at the Oigawa Power Station. Penstock butter- 
fly valves have their own operating oil supply, operating motors, 
and controls. The hydraulic control circuit contains, in series 
with the main actuating cylinder, a pressure-operated sequence 
valve which will open and allow the butterfly valve to move only 


when the system is under pressure from the oil pump. A failure 


? Engineer, Technical Engineering Analysis Section, Bureau of 
Reclamation, Denver, Colo. 
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of power to the oil pump motor while the butterfly valve is being 
closed or opened will result in the leaf being hydraulically locked 
in the position attained at the time of power interruption. Fail- 
ure of the operating linkage or the hydraulic cylinder is essen- 
tially the only way the butterfly valve leaf could slam closed. 

The pressure regulators, where used by the Bureau, operate 
essentially independently of the governor oil supply for operation 
of the turbine wicket gates. The pressure regulator valve is 
opened by a linkage to the servomotor system which closes the 
wicket gates. The linkage is strong enough to stall the servo- 
motors if the regulator valve is stuck shut. When used as a 
water-saving valve, the closing force for the pressure regulator 
valve is supplied by turbine scroll case pressure. The rate of 
closing of the regulator valve is controlled by an oil filled dashpot. 
Thus, it is virtually impossible for the wicket gates and the pres- 
sure regulating valve to be open at the same time except for short 
periods of time. Loss of oil to the servomotor system could 
leave the wicket gates partly open, but the pressure regulating 
valve, being self-contained, would still close at a controlled rate. 


Robert A. Sutherland® 

The author and the Japanese Operating Company are to be 
congratulated on making available the information contained 
in this paper. Too frequently, failures are “‘soft-pedalled’’ or 
not adequately described in the technical press. 

As the author points out, a proper operating procedure would 
have avoided the accident. A system of interlocks between the 
manifold connection valves and the main penstock valves might 
have achieved the same result by ensuring correct operation. 
It is doubtful whether more adequate stiffeners on the penstock 
would have avoided the collapse in sections 15 to 21 since the 
collapse carried through the anchor block at section 20. 

The events described in this paper emphasize the importance 
of an unfailing source of oil pressure where oil-operated valves 
are concerned. 


Author's Closure 

Mr. Bingham’s discussion again emphasizes the need for proper 
operator training, the lack of which contributed substantially to 
the Oigawa failure and the gate operation accident which he 
describes. Design engineers rightfully concern themselves with 
providing controls which will safeguard the plant against human 
However, the provision of foolproof design features is 
A careless 


errors. 
obviously uneconomical beyond a certain point. 
operator can put any plant out of commission in spite of the de- 
signer’s efforts. An essential feature of safe, efficient operation is 
proper personnel training and organization, both of which had 
been neglected in the plants under discussion here. 

Many mechanical and electrical features were also obviously at 
fault in the Oigawa design. Mr. Crawford has pointed out a 
number of Bureau of Reclamation features which would make 
Oigawa type malfunctionings impossible. In addition to those 
described there are float valves at modern oil pressure tank out- 
lets which will close on oil supply failure and keep air out of 
cylinders of oil operated valves and controls. The Oigawa pres- 
sure tanks were apparently not equipped with such valves. The 
Bureau of Reclamation features and variations have been used in 
the many hydropower plants designed by the author’s firm. 
The possibility of a series of malfunctions is, of course, kept in 
mind and every effort is made to eliminate combinations which 
may result in damage to equipment or plant. 

Mr. Sutherland has mentioned the ‘soft pedalling’’ of serious 
failures, free discussion of which would greatly advance engineer- 
’ Hydraulic Engineer, Ebasco Services, Inc., New York, N. Y. 
Mem, ASME. 
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ing science. Weof the ASME and ASCE who are concerned with 
water hammer are indebted to Mr. Sakihara and his associates for 
permitting publication of this particular failure. 

It should be noted in this connection that there are many in- 
teresting and informative reports of failures readily available in 
foreign technical literature. Because of language difficulties, 
American engineers do not keep abreast of much of this valuable 
information. Moreover, the foreign technical publications carry 
extremely interesting articles on more general subjects than oper- 
ating failures. The writer has had the good fortune of spending a 
number of years in foreign countries, many times a good distance 
away from American technical publications. At these times the 
foreign publications provided an excellent source of enjoyable 
leisure-time reading. 

It is my opinion that we have much to learn from those who 
write in foreign technical journals; their methods of presentation 
in many cases are quite superior to ours. The obvious difficulty, 
of course, is translations. It appears we could profit considerably 
by the publication in this country of the best of the foreign papers. 
Just how the engineering societies could go about obtaining 
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translations of these articles is a difficult problem. Perhaps those 
in the various societies who do read articles in other languages 
could be encouraged to translate them—or at least recommend to 
the Publications Committees of the Society the best articles for 
translation; then, it might be possible that the authors would be 
willing to furnish a copy of the paper in English for publication. 
In addition, foreign authors might be encouraged to submit papers 
in English to our societies. An interesting practice in several of 
the foreign countries is that of paying the author of a technical 
paper when it is accepted for publication. While this could 
hardly be done in this country, there might be very little difficulty 
in paying the relatively small amounts normally called for in other 
countries to authors of such papers. It is interesting to note that 
a great number of our technical papers are translated into Japa- 
nese, Chinese, French, Italian, and Spanish by the engineers in 
those countries. 

In any event, I firmly believe that a greater flow of technical 
information between the countries of the world will assist us all in 
getting along together; however, we certainly do not need any 
more conventions or committees for this purpose. 
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and are discussed. 


Introduction 


iat THE PROBLEMS with which a centrifugal 
pump designer must cope is the unbalanced radial] load on the 
impeller caused by the nonuniformity of the static pressures and 
velocities within the volute, particularly at off-design conditions. 
Unless shaft and bearing designs are modified to compensate for 
higher loads, the resultant shaft deflection due to this radial] load 
can sometimes be of sufficient magnitude to cause rubbing at 
close running clearances with the consequence of premature wear 
on the parts in contact. Furthermore, the load is cyclic because 
of rotation and can lead to shaft failure due to fatigue. Hence 
it is important that the radial forces acting on the impeller be 
given consideration in the mechanical design of a centrifuga] 
pump. 

It is surprising that there is relatively little published technical 
information on radial thrust, since its presence in centrifugal 
pumps and its significance on pump performance have long been 
known. In 1913 Lawaczek [1],' while experimenting with a high 
head centrifugal pump, noted excessive bearing temperatures 
after a short running period. He theorized, from pressure meas- 
urements made around the volute, that the fault was due to a 
radial thrust on the impeller caused by the unequal pressures 
around its periphery. The installation of a vane within the 
volute so as to cause the impeller to discharge into similar geo- 
metric regions corrected the difficulty. This was, perhaps, the 
first use of the double volute (Fig. 1) as a means of reducing the 
unbalanced radial thrust on a centrifugal pump impeller. 

Binder and Knapp [2] conducted an investigation on the flow 


sf LG 


DS 
DOUBLE VOLUTE 


VOLUTE DIFFUSER 
Volute designs for minimizing radially unbalanced forces 


Fig. 1 
' Numbers in brackets designate References at end of paper. 
Contributed by the Hydraulic Division and presented at the Hy- 

draulic Conference, Ann Arbor, Mich., April 13-15, 1959, of THe 
AMERICAN Society OF MECHANICAL ENGINEERS, 

Norte: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those 
of the Society. Manuscript received at ASME Headquarters, August 
8, 1958. Paper No. 59—Hyd-2. 
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An Experimental Investigation of Radial 
Thrust in Centrifugal Pumps 


An experimental investigation has been conducted to determine the magnitudes and 
directions of the unbalanced radial forces on centrifugal pump impellers. 
covers single volutes for a wide specific speed range, duuble volutes, concentric casings, 
and modifications of the concentric casing 


The work 


The results are presented in graphical form 


A method, making use of strain gages, was devised for determining the magnitudes and 
directions of the resultant radial forces and is described. 


characteristics within the volute of a centrifugal pump and 
demonstrated by actual measurement the nonuniformity of the 
static pressure and velocity distribution around the impeller at 
off-design conditions. They attributed the radial force on the 
impeller to a combination of a static pressure force and a force 
due to momentum which, added vectorially, would give the re- 
sultant radial thrust on the impeller. Their calculations of the 
magnitude and directions of the radial forces in this manner 
agreed with actual measurements of the shaft deflection. Fur- 
thermore, they found that even at the normal design condition, 
where the variation in static pressure and velocity was slight, an 
unbalanced radial thrust on the impeller was still evident. 

Knapp [3] measured the static pressure distribution around the 
same impeller installed in a conventional volute, a double volute, 
and a volute-diffuser arrangement (Fig. 1). He found that the 
resultant radial forces were greatest for the conventional volute 
at all capacities other than design. The double volute almost 
eliminated the radial thrust whereas the volute-diffuser combina- 
tion gave less thrust than the volute alone but more than the 
double volute. The radial unbalance found with the diffuser- 
volute combination is attributed to the fact that the diffuser still 
discharges into a volute arrangement. The magnitude of the 
thrust is less because the water is discharged into the volute 
casing at a lower velocity than is the case in a conventional vo- 
lute. 

Stepanoff [4] presents relationships which enable the pump de- 
signer to estimate the radial thrust on a centrifugal pump im- 
For a volute casing 

F, = Kpd,W, (1) 


peller. 


where 
F, 
p pump total head, lbp per in.? 
d impeller diameter, in. 
W; impeller width including the shrouds, in. 
K = a parameter which varies with capacity 


radial resultant force, lbp 


' 


To select the proper value of K he gives the following equation 
arrived at experimentally for volute pumps: 


ee Pe Fm 
x ~ 0a [1 - (2) 


where Q is any capacity and Q, is the capacity at the best ef- 
ficiency point of the pump. Equation (2) is for a particular pump 
and assumes that the maximum value of radial thrust occurs at 
shutoff. At the best efficiency point the value of K becomes 
zero, 
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For concentric casing pumps Stepanoff expresses the parameter 


K = 0.36 (2) (3) 


Equation (3) implies that at zero capacity K is zero and hence the 
radial force is also zero. As Q increases K increases as will 
the radial force since it is directly proportional to the constant K. 

The work reported in this paper was confined to single-stage 
units pumping water. The objectives of the experimental pro- 
gram were to determine the effect of specific speed and casing 
geometry on the developed radia] forces for the complete operat- 
ing range of the pumps. Measurements of the magnitude and 
direction of the resultant radial thrust were made on conventionally 
designed volute pumps covering a wide range of specific speeds, 
on double volute collectors, concentric casings, and on modifica- 
tions of the concentric design. 

A radial force measuring system, making us¢ of strain gages, was 
developed and employed for the test work. This force measuring 
system has so greatly facilitated the determination of unbalanced 
radial forces that it has become a part of the rating tests for new 
pump designs. A description of the technique and apparatus 
used is presented under Instrumentation. 


as: 


Instrumentation 


The problem was that of measuring the magnitude and direc- 
tion of radial impeller forces in a centrifugal pump over the nor- 
mal range of operating conditions. The apparatus designed to 
determine these forces had to simulate the actual conditions of 
rotor support in the standard pump allowing only negligible radial 
deflection so that tests could be made with normal] wearing ring 
clearances. 

For earlier test work the mechanical probe method was used. 
This involved the use of four micrometer probes mounted on the 
casing in such a manner that they could be advanced inward until 
slight rubbing against the shaft or impeller occurred. This 
rubbing was detected by electrical contacts or by a listening 
probe. Deflection at the plane of the probes was translated into 
force at the impeller center line by actual weight calibration. 
This method was time consuming, subject to vibration effects, and 
it was somewhat awkward in practice. 

In order to avoid the disadvantages of the probes, a method 





was devised using strain gages to detect bending stresses. In 
this method the resultant radial force acting on one of the two 
shaft support ball bearings is measured. This force is separated 
into two transverse components (horizontal and vertical) each of 
which is measured separately. A special bearing mounting was 
designed for one of the two bearings as shown in Fig. 2 and is es- 
sentially a tubular cantilever beam. One end of the tube was 
rigidly attached to the bearing frame and the other end contained 
the movable bearing. A section of the tube was of reduced wall 
thickness to raise the bending strains to a reasonable and meas- 
ureable range. Eight strain gages were mounted axially in pairs 
at 90-deg spacing and wired into two Wheatstone bridges having 
separated horizontal and vertical planes of sensitivity. By 
measuring the magnitude of the horizontal and vertical com- 
ponents separately, the magnitude and direction of the resultant, 
radial force are thus determined by vector addition. 

The axial spacing of the two bearings and impeller remained 
constant so that the device could be calibrated in terms of 
weights applied at the impeller center line. Not only was a static 
calibration made but a specially designed fixture permitted a 
dead weight calibration at test speed also. Periodic calibration 
checks made during the test work attested to the reliability of the 
force measuring system. 

A mechanical seal was used in place of the conventional stuffing 
box to prevent any support of the shaft by the packing and care 
was taken with the suction piping to insure a uniform velocity 
distribution at the pump inlet. 


Experimental Results and Discussion 


Volute Casings. Experimental data were obtained on sixteen 
conventional design constant velocity volute casing pumps cover- 
ing a range of specific speeds from 25 to 165. Figs. 3 to 6, inclusive, 
are the results of the test work. 

In general, the experimental work showed that the resultant 
radial forces were a maximum at shutoff, decreased with increas- 
ing capacity and reached a minimum near the best efficiency 
point. Further increase of the flow beyond the design condition 
resulted in an increase in the force. 

The shape of the radial force curves appeared to be a function of 
specific speed, Fig. 3. In order to permit comparisons of pumps 
of different sizes, both abscissa and ordinate are in terms of per- 
centage of design capacity and percentage of force at shutoff, 
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Fig. 2 Force measuring system for end suction pumps 
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Fig. 3. Radial force as a function of capacity and specific speed 
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Fig. 5 The effect of speed on radial thrust 
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respectively. Note that for the lower specific speeds there is a 
relatively large capacity range over which the minimum force is 
essentially constant. For the higher specific speeds the forces are 
a minimum almost at a single point and on either side of this 
capacity the rate of increase of the forces is high. 

From the foregoing characteristics it is evident that if the re- 
sultant radial thrust is an indication of the degree of hydraulic 
mismatch between impeller and casing then for the higher 
specific speeds there is a very small capacity range over which the 
two are compatible. For the lower specific speeds the opposite 
appears to be true since there is a relatively large range of ca- 
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180° 
pacity over which the impeller and case exhibit a favorable rela- 
tionship. 

The curves of Fig. 3 further show that the point of minimum 
radial thrust does not always occur at the normal design capacity 
and that the radial force is not always zero at this point. It was 
generally found that for the higher specific speeds (NV, > 100) the 
point of minimum thrust was at a larger capacity than the best 
efficiency point, and for the lower specific speeds (NV, < 50) the 
minimum force occurred at a lower capacity than design. For 
the intermediate specific speeds there was no fixed position of the 
minimum thrust and its location was sometimes at the best ef- 
ficiency point as well as either side of this capacity. 
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As was mentioned under Instrumentation the magnitude and 
direction of the resultant radial force were determined by direct 
measurement of its vertical and horizontal components. Fig. 4 
illustrates the method used for plotting the test results for four 
pumps of different specific speeds. The vertical and horizontal 
axes are in terms of per cent of force at shutoff while the percent- 
age values at the test points are per cent of design capacity. A 
line drawn from the origin to any point on the curves would indi- 
cate the direction and magnitude of the resultant radia] force for 
that capacity. 

While it has been known generally that the direction of the 
radial thrust varies with flow, the experimental results have also 
shown that the directions are not the same at similar flows for 
pumps of different specific speeds. This is readily seen from Fig. 
4. For the higher specific speeds (NV, = 100) the force at shutoff 
is approximately toward the casing cut-water (or casing tongue) 
and the curve for VN, = 112 is representative of the direction of the 
radial thrust for the specific speeds in this range. For the lower 
specific speeds the force at shutoff shifts in angular position in a 
direction corresponding to the impeller rotation. 

That the magnitude of the radial thrust is proportional to the 
head developed by the pump and, consequently, should also vary 
approximately with the square of the speed ratio has been verified 
experimentally and is shown in Fig. 5. Although the results of 
tests on only one pump are shown, similar results were observed 
on other units. 

The radia] thrust factor K has been computed from equation 
(1) for all pumps tested. A plot showing the variation of the 
parameter with specific speed is shown in Fig. 6 for 0, 50, and 190 
per cent of normal design capacity. Although there appears to be 
considerable scatter of the points there is, nevertheless, an evident 
trend for each curve. The dashed portion of each curve repre- 
sents an extrapolation to the highest specific speed pump tested. 

Referring to Fig. 6 again it can be seen that at the best ef- 
ficiency points the K-values are close to zero for the lower specific 
speeds and tend to increase with increasing specific speeds. The 
reason for this, as has been previously mentioned, is that the 
minimum radial force does not always occur right at the point of 
best efficiency. 

Concentric Casing. A circular casing was so designed that the 
same impeller could be tested with various width and radial height 
combinations. It was possible to vary the ratio of casing width 
(W,) to impeller width (W;) from 2.13 to 3.20 and the ratio of 
casing diameter (d.) to impeller diameter (d;) from 1.145 to 1.430. 

Fig. 7 shows the results of some of the combinations tested. 
For a comparison of the magnitude of the radial forces, the results 
of the same impeller tested in a standard single volute are also 
shown. The radial forces for the circular collectors are ex- 
pressed in percentages of the forces developed at shutoff by the 
impeller in the volute casing. Although a great number of data 
were obtained for all the geometric possibilities, only three curves 
are shown since the hydraulic performance for these casings was 
equivalent to that of the conventionally designed volute casing. 
Equivalence in this instance means that not only were the ef- 
ficiencies comparable to the standard volute casing but also that 
the peak values of efficiency occurred at the same capacity. Fig. 7 
is for a constant casing width to impeller width ratio of 3.2 and 
varying ratios of casing diameter to impeller diameter. 

From the results obtained it appears that certain combina- 
tions offer an effective means for minimizing the unbalanced ra- 
dial load on the impeller. As Fig. 7 indicates, the unbalance radial 
forces for the higher diameter ratios (1.290 and 1.43) are much 
lower than the standard volute except in the region of maximum 
efficiency. Even in this range the force is only approximately 
40 per cent of the radial thrust at shutoff for the volute case. 
Note that the force curves for the two higher diameter ratios are 
practically the same, indicating that an optimum dimension has 
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Fig.7 Radial force for concentric casings for constant width and varying 
radial dimensions 


been achieved. Furthermore, this optimization was observed for 
all width ratios down to and including 2.13. Similar force reduc- 
tions were obtained for these diameter ratios also. 

For the 1.145 diameter ratio the relationship between radial 
force and capacity is quite similar to that of the volute casing 
with the exception that the magnitude of the forces is less for all 
capacities up to 110 per cent of normal. For all casing widths, 
the characteristics of the force-capacity curves for the diameter 
ratio of 1.145 were similar to those encountered with volute 
casings. 

It is apparent from the foregoing that the use of equation (3) 
to predict the radial thrust for concentric casings may lead to 
erroneous results. Whereas the radial thrust is a maximum at 
the best efficiency point, the magnitudes of the forces are much 
less than the force at shutoff for the conventional design volute. 
Also, while they do approach a minimum near shutoff, the radial 
forces were never found to be zero at this condition. Further- 
more, it appears that for “shallow” concentric casings the radial 
force characteristics are quite similar to those expected of a 
volute collector. 

The directions of the radial force at shutoff for the circular 
casings were approximately 180 deg from the casing discharge and 
shifted through 180 deg (in the direction of impeller rotation) as 
flow was increased to the maximum. 

Modified Concentric and Double Volute Casings. The double 
volute has been used as an effective means of minimizing the 
unbalanced radial thrust on an impeller. In many cases this is a 
suitable treatment of the problem. However, there are certain 
pump designs and physical size limitations where the casting of a 
vane within the volute can cause manufacturing difficulties. The 
use of a double volute casing in small sizes makes it more difficult 
to duplicate performance. Consequently, a program was under- 
taken to determine if certain casing modifications could be made 
which would result in low radial forces at both shutoff and de- 
sign conditions without adverse effects on hydraulic performance. 

Three casings were designed for impeller specific speeds of 55, 
100, and 165. The collectors are circular for a portion of- the 
angular distance from the casing cutwater before the radial ex- 
tent of the casing wall is increased to the discharge nozzle. The 
values of the ratio of casing width to impeller width for each pump 
together with the ratio of casing diameter to impeller diameter 
for the circular portion of each collector are shown on Figs. 8, 9, 
and 10. 

In order to compare the magnitude of the radial forces de- 
veloped in the modified casings the same impellers were tested in 
a double volute and in a single volute casing. Hence there are 
three tests for each specific speed for three collector designs. 

Figs. 8, 9, and 10 show the results of the radia] force measure- 
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Fig. 8 Comparison of the effect of three casing designs on radial force 
for N, = 55 
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Fig. 9 Comparison of the effect of three casing designs on radial force 
for N, = 
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ments made on the three casings for each specific speed. For 
each specific speed the lowest radial loads on the impeller were 
with the double volute casing. However, for the 


55 and 165 (Figs. 8 and 10) with the impeller 
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Fig. 10 Comparison of the effect of three casing designs on radial force 
for N, = 165 


installed in the modified concentric casing, the reduction in the 
developed radial force is comparable to that obtained by 
the double volutes. Fig. 9 indicates that the radial forces 
for the modified collector are larger than in the double volute 
as flows are reduced toward shutoff, but when compared to 
the standard volute there is still a considerable reduction in the 
magnitude of the radial thrust on the impeller. 

The directions of the thrust on the impeller for the modified 
circular casings are shown in Fig. 11. For the higher specific 
speeds of 100 and 165 the direction of the radial force at shutoff 
is practically the same as that with the single volute pumps in the 
The direction of the thrust at shutoff 
as would be expected 


same specific speed range. 
for the lowest specific speed was the same 
for a fully circular casing. Since the lowest specific speed collec- 
tor was concentric for a greater portion of the casing development, 
it probably is to be expected that it would retain some of the 
characteristics of a fully circular collector. 

It was found that the direction of the radial thrust on the im- 
pellers in the double volutes was generally toward the casing 
tongue. No results of the directions of the forces are shown for 
the double volutes since the variation of direction with capacity 
was random in nature and would show up on a plot as scattered 
points over a small sector of the polar diagram. 
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Polar plot showing the direction of the resultant radial forces for the modified concentric casings at 
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A Summary of the Test Results 


A There are certain characteristics of radia] thrust in volute 
pumps which appear to vary with specific speed. These are: 


1 Fraction of design capacity at which the minimum thrust 
occurs. 

2 Range of capacity over which the forces are close to mini- 
mum. 

3 Direction of the resultant radial forces. 

t Radial thrust factor, K. 


B Concentric casings with certain dimensions of width and 
diameter are quite effective in minimizing the unbalanced radial 
force on a centrifugal pump impeller. 

C The modifications made to the concentric collector result 
in a reduction in radial force which is comparable to the reduction 
obtained with a double volute casing. 
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DISCUSSION 
Fred F. Antunes? 


The authors and their company deserve much credit for the 
extensive work done for the completion of this paper. The con- 
tents of this article certainly add much knowledge for the pre- 
diction of radial thrust in centrifugal pumps. It is interesting 
to note that specific speed, as in many pump studies, was also an 
effective tool in this report. 

In connection with this study the writer would like to add some 
notes on the radial load tests of the Grand Coulee model pump 
which were recorded in the ASME Trans. of 1950. This par- 
ticular pump is of interest since it employed 6 diffuser vanes in 
the casing. Utilizing the test data and formula (1) of this re- 
port, the following table was compiled: 

tadial force, 
per cent ol 
force at 
shutoff for 
volute casing 
17.0 


Model 


sapacity, 


Capacity, 
per cent tadial load, 
cu ft/sec of normal Ib 
0 0 255 
1.52 20 178 11.9 
3.05 10) LO7 vel 
4.57 60 62 x 
6.10 80 35 2.3 
7.60 100 30 0 
9.15 120 36 4 

The specific speed of this unit was 99 (cu ft/sec). Comparing 
the figures in the last column with Fig. 9(N, = 100) of this paper, 
it appears that additional vanes have little effect on decreasing the 
radial load as compared with the two-vane or double volute 
casing. 

Fig. 7 of this paper shows the results of a standard volute casing 
at a specific speed 63. Plotting this same curve on Fig. 3, it is 
noted that this particular test and specific speed does not follow 


the trend as shown. Is it possible that the combinations of 


2 Hydraulic Engineer, Morris Machine Works, Baldwinsville, N. Y. 
Assoc. Mem. ASME. 
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casing and impeller, such as small diameter impeller in given 
casing, lower or larger capucity impellers in given casing, etc., 
require further investigation? 

The considerable reduction in radial thrust illustrated in this 
report for the “‘modified concent™ic casing’’ is very interesting 
since this design appears to approach casings used for pumping 
of materials of large sphere s.ze. The writer was unable to de- 
termine the criterion used for the development of the modified 
casings and wonders if the authors would care to mention this. 


André Kovats* 

Some points of the paper are not clear, such as: Ts cu ft/sec 
used as capacity in the N, or is N, the dimensionless formula? 
Further there seem to be some liscrepancies between Figs. 3 and 
4. For the N, = 37 impeller, Fig. 3 shows a lowest thrust value 
of about 8 per cent of the shutoff thrust at 80 per cent capacity, 
while in Fig. 4 the lowest thrust value is shown as 42 per cent. 

The writer wonders if the impeller-volute ratio, such as im- 
peller diameter to tongue diameter and impeller width to volute 
width was the same for both the low specific speed pumps and the 
high specific speed pumps. Or was it a mismatch between the 
impeller and the volute of the low specific speed pumps that 
caused the high radial thrust at design point? This may be an 
explanation for the surprising result, where an impeller in a cir- 
cular casing has the same efficiency as in a volute casing. Former 
experiences have always shown a higher efficiency with correctly 
designed volutes than with circular casings. 

A valuable improvement seems to be in the “modified circular 
casing’? which may also have other advantages beside the con- 
siderable reduction of the radial thrust. 


Donald S$. Ullock* 


This paper reflects a rather extensive amount of experimental 
work by the authors on the measurement of radial thrust in cen- 
trifugal pumps. Values for radial load factors over a range in 
specific speeds are presented and in some cases information on 
design parameters is provided, but it is to be regretted that the 
authors were not able to give more design information about the 
pump casings and impellers tested. In a broad sense, their data 
confirm the experimental data and empirical formulas available 
in the literature for estimating radial thrust developed by im- 
pellers. 

Since it is customary to express the specific speed for cen- 
trifugal pumps in terms of rpm, the authors should have defined 
their units for specific speed as being rps of speed, cfs of capacity, 
and ft of head. 

In Figs. &, 9, and 10, ratios are shown for the casing and im- 
peller widths and for the casing and impéller diameters for com- 
parison of the double volute and the modified concentric casings 
with the so-called ‘‘standard volute” casing. No ratios are given 
for the other important parameter, cutwater to impeller radius. 
None of these ratios are shown on Figs. 3, 4, 5, and 11. 

References are made to the early work of Lawaczek, the work 
of Binder and Knapp, a later paper by Knapp, and Dr. Stepanoff’s 
book on ‘Centrifugal and Axial Flow Pumps.’”’ No references 
are made to the paper by Ullock on ‘Evaluating the Mechanical 
Design of End-Suction Centrifugal Pumps’’® and the paper by 
Ullock, Reynolds, and Hudson on ‘‘Mechanical Design of Cen- 
trifugal Pumps From Shaft Deflection Measurements.’’® These 
two papers were written to stimulate interest in adequate mechani- 

3’ Chief Engineer, Pumps, Foster Wheeler Company, New York, 
N. Y. Mem. ASME. 

4 Staff Engineer, Union Carbide 
Charleston, W. Va. Mem. ASME. 

5 Chemical Engineering Progress, vol. 51, 1955, p. 207. 

6 Ibid., vol. 52, 1956, p. 35. 


Chemicals Company, South 
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cal design in centrifugal pumps. Examples and measurements 
were presented that demonstrated an apparent lack of attention 
to the presence of unbalanced radial loads created in the pump 
casing. 

The plot on Fig. 6 for volute casings shows that the maximum 
radial thrust factor decreases with a decrease in specific speed 
from what appears to be a constant value of about 0.36 (for spe- 
cific speed values greater than 110) down to about 0.08 at a spe- 
cific speed of 20. With these data, no values are presented for 
the parameters—cutwater to impeller radius and casing to im- 
peller width ratios. As shown on Fig. 6 and for different combi- 
nations of these ratios, can it be stated that the radial thrust 
factors will decrease as the specific speed is decreased? 

On Fig. 7, the radial forces developed in concentric casings are 
compared with the forces developed in a volute casing. It is 
presumed that for both concentric and volute casings the spe- 
cific speed is 63 and the casing to impeller width ratio is 3.2. It 
is shown that for the largest concentric casing to impeller di- 
ameter ratio of 1.430, the radial force at the best efficiency point 
is 35 per cent of that for the volute casing at zero capacity. This 
means that the radial thrust factor for this concentric casing is 
0.35 x 0.2 or 0.07 at its maximum efficiency point. 

In evaluating the mechanical design of centrifugal pumps with 
concentric casings, I have received experimental data (from 
several pump manufacturers) on radial thrust for some twenty 
pumps of different sizes. At the best efficiency point, the ca- 
pacities ranged from 75 gpm to 1700 gpm and the heads from 50 
to 190 ft. When the radial thrust factors are plotted against 
the specific speeds, a shot-gun pattern is obtained. The specific 
speeds varied from 22 to 73 and the factors from 0.1 to 0.3 while 
one manufacturer’s pump had a factor of 0.5 for a specific speed 
of 25. 

The conclusion is drawn by the authors that the use of Stepa- 
noff’s equation (K = 0.36 Q/Qn) for radial thrust factor in con- 
centric casings “‘may lead to erroneous results.’”’ However, if we 
are interested in good mechanical design, we are interested in the 
maximum loads that may develop. In general and if no ex- 
perimental data are available on the pump under consideration, 
the only error that can result from using this equation would be 
in equipping the pumps with heavier shafts and larger bearings. 
Only the extra cost will hurt the user. In these cases, mechanical 
overdesign with the extra investment cost will be much less ex- 
pensive than underdesign with multiple repair costs. 

This paper is of special value in showing that a centrifugal 
pump casing can be modified so that its design is intermediate 
with respect to typical volute and concentric casings. The au- 
thors demonstrate that radial thrust will be minimized and could 
be almost constant over the capacity ranges for pumps in the 
specific speed range of 55 to 165. This possibility was discussed 
in the Appendix to the paper by the writer.’ 
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Authors’ Closure 

The authors wish to thank Messrs. F. Antunes, A. Kovats, and 
Dr. Ullock for their comments and the interest they have shown 
in the paper. 

The Grand Coulee model pump data introduced by Mr. An- 
tunes is an important addition to the paper and it does appear 
that additional vanes have little effect on decreasing the radial 
load when compared to the double volute casing. 

The authors suggest that if Mr. Antunes had plotted the results 
of Fig. 7 for the standard volute casing on the upper curve group- 
ing in Fig. 3 he would have obtained a curve which would follow 
more closely the trend of that group. 

As far as the design of the modified casings are concerned, the 
width and diameter ratios are indicated in Figs. 8, 9, and 10. 
In addition to these parameters the collectors were made circular 
for 270 deg and the remainder of the development approximated 
a volute. 

To clarify the units of specific speed for Mr. Kovats, rpm was 
used as shaft speed, total head was in foot-pounds force per pound 
of mass, and capacity was in cubic feet per second. 

The discrepancies between Figs. 3 and 4, as mentioned by Mr. 
Kovats, are due to the fact that different pumps were used in each 
of the plots. In the lower specific speed ranges there is less 
consistency in the minimum radial force expressed as a per- 
centage of shutoff force than in the higher specific speed range. 
Fig. 3 was meant to show the relative shape of the low and high 
specific speed radial force curves. If the results of Fig. 4(N, = 
37) were plotted on Fig. 3, a curve very similar to N, = 25 would 
be obtained. 

The impeller volute ratios were not the same for pumps of dif- 
ferent specific speeds and no hydraulic mismatch was evident 
since the performance of these units in every case met the designed 
efficiencies. Furthermore, the authors take exception to the 
statement that volute pumps have always shown higher ef- 
ficiencies than circular casings. 

Dr. Ullock has mentioned the lack of design ifnormation, par- 
ticularly the absence of the cutwater to impeller radius ratios for 
the pumps tested. During the evaluation of the test results an 
attempt was made to use this ratio in the analysis and in every 
case the plots showed so much scatter that this parameter was not 
considered to be too significant. The plot of K-factor versus 
specific speed as shown in Fig. 6 yielded the most promise for the 
volute pumps. Since specific speed does classify a pumping unit 
on the basis of performance regardless of physical size, and since 
the plot of the data with this parameter definitely indicates a 
trend, it is felt that Fig. 6 is a useful design tool. 

For the concentric casings comparison with the volute pump 
(Fig. 7) was made and shown only for those circular casings whose 
hydraulic performance and not geometry was equivalent to the 
standard volute pumps. 
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Design and Test‘of Mixed-Flow and 


JAMES J. KRAMER? 
WALTER M. OSBORN? 
JOSEPH T. HAMRICK? 


Centrifugal Impellers 


Procedures for designing the hub-shroud contour and the blade shape of mixed-flow 
and centrifugal impellers are presented and evaluated in terms of experimental results 
obtained for impellers designed by these methods. 


The experimental results in general 


indicate that these design procedures are reliable. Additional experimental results with 
varying boundary-layer allowance and some vaned diffuser test results are presented. 


Introduction 


NTEREST IN THE FIELD of mixed-flow and centrifugal 
impellers resulted in the undertaking of a program to develop re- 
liable design procedures that would obviate the need for lengthy 
development work. This program was based on the philosophy 
that proper control of the fluid velocity on the wetted surfaces of 
the machine was the key to good over-all performance. The 
procedure followed, in general, was to develop a method that 
enabled the designer to control the wetted surface velocities and 
then to investigate the reliability of the design method by over-all 
performance tests of specific impellers designed by the method. 

All the design procedures were based on isentropic flow models 
so that real fluid effects were neglected in the design. To com- 
pensate for this, experiments were conducted in order to determine 
what modifications to the isentropic design shape were required 
to compensate for viscous effects. 

Blade design procedures were also used to produce diffuser 
vanes, and tests were conducted to determine the performance of 
these diffusers in combination with one of the mixed-flow im- 
pellers. 

The research program summarized in this paper was originally 
reported in [1-8].* 


Development of Design Procedure 
Blade Design With Meridional-Plane Analysis. The first design pro- 


1 Branch Consultant, Turbopump Components Branch, Lewis 
Research Center, National Aeronautics and Space Administration, 
Cleveland, Ohio. 

? Aeronautical Research Engineer, Turbopump Components 
Branch, Lewis Research Center, National Aeronautics and Space 
Administration, Cleveland, Ohio. 

*Chief Research Engineer, Thompson Products, Rocky 
Mount, Va. 

* Numbers in brackets designate references at end of paper. 

Contributed by the Compressor Subcommittee of the Hydraulic 
Division and presented at the Hydraulic Conference, Ann Arbor, 
Mich., April 13-15, 1959, of Tam AMERICAN Society OF MECHANICAL 
ENGINEERS. 

Note: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those 
of the Society. Manuscript received at ASME Headquarters, 
December 18, 1958. Paper No. 59—Hyd-20. 


Inc., 


cedure used consisted of blade element design using the method of 
[9] and the subsequent analysis of the meridional-plane flow using 
the method of [10]. The details of the aerodynamic design are 
given in [1]. Briefly, the blade design method consisted of the 
following steps: 


(1) The following were prescribed: 


(a) Mean-blade-height line shape in the meridional plane 

(b) Design flow rate 

(c) Blade thickness at inlet and outlet 

(d) Blade height distribution 

(e) Blade surface velocity distributions at mean-blade- 
height line 

(f) Impeller-outlet mean-line speed 


(2) A distribution of blade height h was prescribed. (The 
blade height at inlet and outlet were fixed by prescribed condi- 
tions). The prescribed blade-height distribution along with the 
mean-blade-height line shape resulted in a hub-shroud shape. 

(3) The method of [9] was used to determine the number of 
blades and the blade shape at the mean-blade-height line cor- 
responding to the prescribed conditions. The blade shape at 
positions other than the mean-blade-height line was determined 
by prescribing radial blade elements and blade taper through 
the blade shape at the mean-blade-height line. 


The input required by the method was adjusted and the calcula- 
tions repeated until a satisfactory blade shape was obtained. 
Approximately six to eight hours of hand computing time were re- 
quired to obtain a provisional blade shape from the prescribed 
conditions. In order to determine hub and shroud fluid velocities 
and the blade surface velocities at positions other than the mean- 
blade-height line, the stream-filament method of [10] was used 
to investigate the flow in the hub-shroud or meridional plane. 
This analysis will indicate meridional-plane velocity distributions 
and hence warn the designer of velocity distributions conducive to 
poor performance or of flow areas that are too small to pass the 
design weight flow. 

The meridional-plane analysis was extended to yield also the 
blade surface velocities. This extension was accomplished by the 
assumption of a linear pressure distribution from blade to blade. 
The meridional-plane flow analysis method is an iterative proce- 





Nomenclature 


L = ratio of distance from inlet to any 
point along meanline or stream- 
line to total distance from inlet 
to outlet 

Q = ratio of relative velocity to inlet 
stagnation speed of sound 
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actual impeller speed at outlet 
mean-blade-height line, ft/sec 

total compressor flow rate, lb/sec 

ratio of total pressure at inlet to 
NACA standard sea-level pres- 
sure (29.92 in. Hg abs) 


adiabatic temperature-rise  effi- 


ciency 


ratio of inlet stagnation tempera- 
ture to NACA standard sea- 
level temperature (518.4 R) 
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dure, and in order to obtain good accuracy approximately 500 
hours of hand computing time were required to arrive at the final 
The amount of computing time involved is the chief 
This time can be decreased by the 


design. 
drawback of this procedure. 
use of automatic computing equipment, but the procedure is not 
an easy one to program because of the large amount of engineering 
judgment required. This design procedure was used to design two 
mixed-flow impellers, designated the MFI-1 and the MFI-2, 
which are discussed in the following sections. 

Design and Test of MFI-1. The specified design conditions for 
the MFI-1 are shown in Table 1, and the prescribed blade surface 
velocity distributions are shown in Fig. 1. These specifications 
resulted in the blade shape shown in Fig. 2. The provisional and 
the modified hub-shroud profiles (based on the stream-filament 
analysis) are shown in Fig. 3. 

This impeller was designed for a 10 per cent increase in relative 
‘velocity at the mean-blade-height-line section from inlet to out- 
let, which it was hoped would result in accelerating flow along the 
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Fig. 1 Prescribed blade surface velocity distribution for design impeller 


MFI-1 


Fig. 2 Impeller passage on developed cone surface 
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hub. Velocity gradients on the blade surfaces were controlled to 
maintain low deceleration rates. 

In order to test the design procedure it was decided to build 
and test the MFI-1 impeller. It was clear, however, that some 
allowance in available flow area had to be made in order to ac- 
count for real fluid effects. The blade height was increased by a 
percentage of the outlet blade height varying linearly with dis- 
tance along the shroud from 0 at the inlet to 74 at the outlet. In 
addition, because of the low ratio of blade height to passage width 
that resulted from the prescription of accelerating flow from inlet 
to outlet, it was decided to insert splitter vanes near the outlet. In 
succeeding tests the blade-height increase was reduced. The re- 
sults are presented in the section entitled “Investigation of Al- 
lowance for Viscous Effects.’”? The resulting modified impeller is 
shown in Fig. 4. This impeller, as well as all others mentioned 
in this paper, was tested in the rig described in the appendix. 
Instrumentation and data precision are discussed also in that 
section. The results of the over-all performance are shown in 
Fig. 5. 

The performance map indicates quite good performance for this 
modified impeller, with a total-pressure ratio of 4.0 obtained at 
design speed with an efficiency of 83 per cent. Complete per- 
formance data for the MFI-1A are given in [2]. 

Design and Test of MFI-2. The next step in the program was 
the design of a second mixed-flow impeller designated the MFI-2. 
The procedure followed was the same as for the MFI-1; the 
specifications for the MFI-2 are shown in Table 2. 


Table 1 Design conditions for impeller model MFI-1 


Impeller-outlet mean-line speed (7.00-in. outlet mean 
radius), ft/sec 

Design flow rate, lb/sec 

Slip factor (from backward-curved blades at outlet) 

Axial inlet velocity, ft/sec 

Relative velocity at mean inlet radius, ft/sec 


Table 2 Design conditions for impeller model MFI-2 
Impeller-outlet mean-line speed (7.00-in. outlet mean 
radius), ft/sec 
Design flow rate, lb/sec 
Slip factor (based on blade-outlet angle equal to zero) 
Axial inlet velocity, ft/sec 
telative velocity at mean inlet radius, ft/sec 
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Fig.3 Hub-shroud profile of design impeller MFI-1 and modified impeller 
MFI-1A 
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inverse blade taper was used. In addition, this impeller had larger 
ratios of blade height to passage width than the MFI-1. The 
blade mean-line section and the hub-shroud profile of the MFI-2 
are shown in Figs. 6 and 7, respectively. The blade height was 
modified in a manner similar to that for the MFI-1; and the re- 
sulting impeller, designated the MFI-2A, was fabricated and 
tested. A photograph of the MFI-2A is shown in Fig. 8 and an 
over-all performance map in Fig. 9. The peak pressure ratio and 
maximum adiabatic efficiency at the design equivalent speed of 
1400 feet per second were 3.95 and 79 per cent, respectively. The 
resulting performance characteristics were considered adequate. 
On the basis of these tests, the design procedure could be con- 
sidered reliable but time-consuming. 

Blade Design With Rapid Approximate Method for Design of 
Hub-Shroud Profile. As noted previously, the analysis of the flow in 
the meridional plane required about 500 hours of hand-computing 
time for the MFI-1. Further, this computation could not readily 


be mechanized. Hence in order to make the design procedure a 


Fig. 4 Modified impetier MFI-1A 
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Fig. 5 Over-all performance characteristics of MFI-1A impeller at inlet 
air pressure of 14 inches of mercury absolute 


con on 
@ 
i 
wo 





CENTER LINE OF ROTATION eee 





is 
Fig.7 Hub-shroud profile of design impeller MFI-2 and modified impeller 
MFI-2A 


—— 


ANGULAR 
VELOCITY. w 


Fig. 6 MFI-2 impeller passage on developed cone surface 


The resulting design was an impeller considerably different in 
physical appearance from the MFI-1. This was done purposely 
in an attempt to show that proper control of velocity gradients in 
the impeller would result in good performance regardless of pas- 


sage shape. Because of the thick blade sections at the meanline, Fig. 8 Modified impeller MFI-2A 
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more useful engineering tool, the calculation time in the meridi- 
onal-plane design had to be decreased. This was accomplished 
by the development of the rapid approximate method of design in 
the meridional plane reported in [11]. In this method the shape 
of a base streamline (usually the hub) and the velocity along this 
streamline are prescribed. The blade shape is completely pre- 
scribed. The procedure is then to work away from this base 
streamline by solving for the flow conditions across an incre- 
mental stream tube using the basic flow equations. Having ob- 
tained the shape and the velocity along the adjacent streamline, 
the procedure is repeated using the newly found streamline as the 
base streamline. In this manner stream tubes are stacked until 
the desired weight flow is passed. 


the analysis method, the blade surface velocities can be solved 


In a manner similar to that for 


for by assuming a linear pressure variation from blade to blade. 
An alternative procedure for obtaining the blade surface velocities 
is to use the approximate analysis method of [12] based on the 
concept of circulation. 

The blade shape that is prescribed completely at the beginning 
of the design can be obtained by the method of [9], as was done 
for the MFI-1 and the MFI-2, or can be prescribed arbitrarily. 
In the latter case the velocities along the blade surface are 
checked at each streamline and undesirable trends are 
eliminated by changes in the blade or hub shapes. 

This rapid approximate method was a considerable improve- 
ment in respect to computing time required, and it was decided 
to test its value by redesigning several existing centrifugal im- 
The procedure was to prescribe the velocity along the 
This was 


any 


pellers. 
existing hub and then to redesign the shroud shape. 
done for three impellers—the parabolic-bladed, the circular- 
bladed, and skewed-paraboiic-bladed impellers. These impellers 
had been designed and tested as part of a previous program in- 
vestigating the effect of blade shape, reported in [13], and hence 
performance data were available. 
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Fig. 9 Over-all performance characteristics of MFI-2A impeller at inlet 
air pressure of 14 inches of mercury absolute 
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Fig. 10 Hub and shroud velocities for redesigned (isentropic) parabolic, 
circular, and skewed-parabolic bladed impellers 
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(a) Parabolic-bladed impeller 


Fig. 11 Modified impellers 


These impellers were redesigned with only moderate decelera- 
tions along the hub. For a given hub shape it is sometimes dif- 
ficult to prescribe accelerating velocities along the entire hub and 
at the same time to maintain desirable velocity distributions along 
the blade surfaces and shroud. The results of these compromises 
were the velocity distributions shown in Fig. 10. In addition, the 
velocity distribution for the unmodified impeller is shown for the 
parabolic-bladed impeller. The analysis of the parabolic-bladed 
impeller using the hub-shroud analysis method was reported in 
[10], and hence complete analysis information was available. The 
hub and shroud velocities for the unmodified parabolic-bladed 
impeller are shown in Fig. 10(a). From this figure it can be seen 
that a much more desirable velocity distribution was designed 
for on the hub than was indicated by the analysis for the original 
impeller. 

All original dimensions of the impellers, except the shroud con- 
tour, were retained. These modified impellers, shown in Fig. 11, 
were tested in order to determine whether the design method did 
produce an impeller with good over-all performance. Good over- 
all performance results would tend to indicate good velocity con- 
trol within the impeller. 

The over-all performance maps of the modified impellers are 
shown superimposed on the performance maps of the unmodified 
impellers in Fig. 12. In Fig. 13 the maximum adiabatic efficiency 
is shown as a function of tip speed. The modified impellers were 
investigated over a range of equivalent speeds from 900 to 1500 
fps and flow rates from maximum to the point of incipient surge. 
These maps indicate that the modified impellers had better per- 
formance characteristics than the original impellers at all speeds 
investigated, the greatest gains occurring at speeds of 1300 fps 
and higher. During the design procedure it was noted that in the 
original impellers there was probably present the tendency for for- 
mation of a potential-flow eddy on the hub and driving surfaces. 
(Eddy formation is taken as beginning when the theoretical veloc- 
ity on the blade surface becomes negative). In the redesign, flow 
conditions favorable to the formation of such eddies were avoided 
in so far as possible. The stable operating characteristics of the 
redesigned impeller as contrasted with the unstable characteristics 
of the original impellers were probably the result of eliminating 
conditions conducive to eddy formation in the modified im- 
pellers. 
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(c) Skewed-parabolic-bladed impeller 


Fig. 11 Modified impellers 


At 1300 fps, the peak pressure ratio and maximum adiabatic 
temperature-rise efficiency for the parabolic-bladed impeller were 
3.07 and 0.825, respectively. For the same condition, the circular- 
bladed impeller and the skewed-parabolic-bladed impeller had 
pressure ratios of 3.13 and 3.15 and efficiencies of 0.737 and 0.805, 
respectively. Of the three, the parabolic-bladed impeller had the 
highest maximum efficiencies (0.854 to 0.800) and the best weight- 
flow range over the speed range tested. 

On the basis of the experimental results obtained with the re- 
designed impellers, it would seem that the general design criteria 
(i.e., minimizing decelerations and avoiding potential-flow eddies) 
and the design procedure used could be considered reliable. 
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(C) SKEWED-PARABOLIC-BLADED IMPELLERS 


Fig. 12 Performance of original and modified parabolic, circular, and 
skewed-parabolic-bladed impellers with vaneless diffuser 





Investigation of Allowance for Viscous Effects 


In the design procedures discussed previously, the assumption 
was made throughout that the flow process was isentropic. In 
order to account for real fluid effects, the blade height was made 
larger than was indicated by the isentropic design procedures. 
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Just how much larger the blade height should be made at a par- 
ticular point along the axis of the impeller to allow for the bound- 
ary-layer build-up to that point is a difficult issue. There is con- 
siderable question whether increasing the blade height is a good 
way of allowing for viscous effects. This problem in general has 
not been attacked in a systematic way for mixed-flow and centrifu- 
gal impellers. In order to contribute some information in this 
general area, a particular way of allowing for viscous effects was 
used in conjunction with the design procedures previously out- 
lined. The blade height was increased by a percentage of the isen- 
tropic blade height at the outlet, varying from 0 at the inlet to 74 
at the outlet. This gradual increase in blade height was made for 
an assumed gradual build-up in boundary-layer thickness along 
the shroud. It was supposed that the boundary layer of greatest 
thickness that on the Therefore all the 
boundary-layer allowance was made at the shroud. This pro- 
cedure is an approximation to the idealized procedure, where a 


would be shroud. 


boundary-layer allowance would be made on all surfaces so that 
the resulting free-stream flow would be the same as that for the 
isentropic design case. Detailed information for predicting the 
boundary-layer behavior on all such internal flow surfaces is not 
available. 

In order to investigate how much of a boundary-layer allowance 
should be made in an impeller where all the allowance is made at 
the shroud, over-all performance was obtained for the same impel- 
The 


isentropic design impellers were specified as MFI-1 and MFI-2. 


ler with varying amounts of boundary-layer allowance. 


The impellers with a boundary-layer allowance at outlet of 74 per 
cent of blade height were called the MFI-1A and the MFI-2A. 
The performance results for the MFI-1A and MFI-2A are shown 
in Figs. 5 and 9. These impellers were then cut down so that the 
boundary-layer allowance was reduced to half its initial value. 
These impellers were called the MFI-1B and MFI-2B. The per- 


formance maps for these two impellers are shown in Fig. 14 (a 
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Fig. 13 Variation of maximum adiabatic efficiency with speed for modi- 
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and b), For the MFI-2 impeller only, the boundary-layer al- 
lowance was eliminated completely in version MFI-2C; the per- 
formance results are shown in Fig. 14(c). In addition to the over- 
all performance curves, shroud static pressures were measured 
and impeller exit surveys were made. The results of weighing 
the various data in order to determine what was happening to the 
internal flow in the impeller are presented here as digested from 
[4]. 

The maximum over-all adiabatic efficiency of the MFI-1A at 
the design speed of 1400 fps was 0.83, and of the MFI-1B was 
0.85. The maximum over-all efficiency of the MFI-2A at the de- 
sign speed of 1400 fps was 0.79; of the MFI-2B, 0.79; and of the 
MFI-2C, 0.77. 

The boundary layer appears to build up rapidly near the inlet 
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Fig. 14 Over-all performance characteristics of MFI-1 and MFI-2 im- 
pellers with blade height changes at inlet-air pressure of 14 inches of mer- 
cury absolute 
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(a) 24-Vane diffuser 


(b) 40-Vane diffuser 


Fig. 15 Impeller model MFI-1B and vaned diffuser 


of the impeller, not gradually as was expected. Thus the assumed 
flow model seems to have been invalid. Hence a means of in- 
creasing the blade height to allow for boundary-layer build-up 
cannot be recommended from the data obtained. 

The data indicated that the consequences of a thickened or 
separated boundary layer depend not only on the design velocity 
gradients but also on the shape of the passage. Changing the 
deceleration rate on the shroud by reducing the blade height did 
little to increase the over-all efficiency. The theoretical velocity 
distribution and the outlet surveys indicate that the increase in 
efficiency fromm the MFI-1A to the MFI-1B was due to a change 
from decelerating to accelerating flow along the hub rather than 
from any change along the shroud. 


Vaned- Diffuser Investigations 
All tests to determine the over-all performance of the various im- 
pellers were conducted with vaneless diffusers. In addition, some 
vaned-diffuser tests were conducted with the MFI-1B impeller. 
Experimental results are presented for two vaned diffusers of 24 
These diffusers are shown in Fig. 15 with the im- 
From the photographs it can be seen that the 


and 40 vanes. 
peller MFI-1B. 
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Fig. 16 Over-all performance characteristics of MFI-1B impeller with 
vaned diffusers at inlet-air pressure of 14 inches of mercury absolute 
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Fig. 17 Over-all performance characteristics at inlet-air pressure of 14 
inches of mercury absolute of MFI-1B impeller with vaned diffusers. Inlet 
vane angle increased 5 deg. 


vaned sections were preceded by a vaneless section in which the 
flow was decelerated primarily by reduction of the tangential 
component of velocity. Both sets of vanes had essentially the 
same velocity distribution on the blade surfaces at the mean- 
blade-height line and performed the same amount of turning. 
However, because of the difference in blade height, a smaller 
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number of blades was required with the larger blade height. These 
blades were designed by the method of [10]. 

The results of the over-all performance tests with the two sets 
of diffuser vanes are shown in Fig. 16. For both diffusers at the 
design speed of 1400 fps, the maximum pressure ratio was 3.73 
with an efficiency of 0.77; a maximum efficiency of 0.78 and a 
pressure ratio of 3.25 were obtained at 1300 fps. These results in- 
dicate a decrease in efficiency of 0.08 and a decrease in total- 
pressure ratio of 0.27 at the design speed of 1400 fps from that ob- 
tained with a vaneless diffuser. However, the vaned diffusers 
resulted in deceleration of the flow to a Mach number of 0.30 as 
compared with deceleration to 0.70 in the vaneless diffuser. 

As would be expected, the weight-flow range with the vaned 
diffuser was much shorter than that obtained with a vaneless 
diffuser at each operating speed. However, pivoting the vanes 
of the 40-vane diffuser 5 deg and bending the vanes at inlet of the 
24-vane diffuser 5 deg to match a reduced impeller weight flow at 
the vane inlet, increased the weight-flow range from choke to 
surge at design speed from 0.6 to 1.2 pounds per second for the 
40-vane diffuser and from 1.2 to 2.0 pounds per second for the 24- 
vane diffuser. The results of the tests of the vanes with inlet blade 
angles changed by 5 deg are shown in Fig. 17. 


Summary of Results 

The design procedures investigated are capable of producing 
good impeller designs with a nominal amount of calculation time 
required for the rapid approximate method and considerably more 
for the stream-filament method. The allowance for viscous effects 
is not well understood. However, the results:cf the experimental 
work indicate that a good level of performance can be obtained 
using the methods discussed herein. However, detailed studies of 
the flow in the impeller seem to indicate that the boundary-layer 
build-up does not occur in the manner supposed, so that a different 
means of allowing for the boundary-layer build-up might well re- 
sult in better performance. 

Tests with vaned diffusers indicate that, in order to avoid the 
performance penalties incurred by the use of vaned diffusers, 
more information is needed concerning viscous effects and dif- 


fuser-impeller interactions. 


APPENDIX 
Apparatus and Instrumentation 
The over-all experimental setup used in these investigations is 
illustrated in Fig. 18, which shows the installation of the MFI-1 
impeller. Air passed through a submerged adjustable orifice, 
through a filter tank, and then into a stagnation chamber that 
housed the inlet instrumentation. From the stagnation chamber, 
the air passed through a converging inlet, the impeller, a vaneless 
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Fig. 18 Experimental setup used in investigation of mixed-flow and 
centrifugal impellers 
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diffuser, and into a collector; it was then discharged into the 
laboratory exhaust system. Constant-area vaneless diffusers 
were generally used for evaluating the performance of the im- 
pellers. However, vaned diffusers were used in two series of tests. 

The inlet pipe, collector, and outlet pipe were insulated to 
minimize the effects of heat transfer, an especially important con- 
sideration during the refrigerated-air runs. 

The weight flow through the compressor was measured by 
means of a submerged adjustable orifice. Temperatures were 
measured by use of calibrated bare-wire spike-type thermocouples 
in conjunction with a calibrated potentiometer. Total pressure 
was measured by use of Kiel-type total-pressure probes, and all 
static taps were 0.030 in. in diameter. All pressures were 
measured by use of water or mercury manometers. The precision 
of the measurements is estimated to be within the following limits: 


+0.04 
+2.0 
£1 
+0.30 


Pressure, in. Hg abs 
Temperature, deg F 
Air-weight flow, per cent 
Speed, per cent of design 


A more detailed description of the instrumentation, test pro- 
cedure, and computations may be found in [2]. 
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DISCUSSION 
Robert C. Dean, Jr.° 


The authors are to be congratulated for summarizing final 
results of the NACA centrifugal compressor program. Since 
experimental verification of the usefulness of the theoretical 
work carried out under this program is particularly searce, this 
paper is a welcome addition. 

Impressive efficiency values are presented in Figs. 5, 7, 9, 13, 
and 14. However, it is not stated how these efficiencies were 
measured, whether they refer to the impeller only or to the en- 
tire stage. This criticism also applies to Figs. 16 and 17. It 
would be useful also if geometrical descriptions were presented 
for the vaneless diffusers including radius ratios, depths, and the 
shapes of meridional mean lines. 


Authors’ Closure 


The authors wish to thank Dr. Dean for his comments on their 
The efficiency values presented in Figs. 5, 9, 13, amd 14 
were impeller efficiencies. These efficiencies were measureé and 
presented in order to demonstrate the effect of changes in impeller 
geometry on impeller efficiency. For the mixed-flow im- 
peller tests (Figs. 5, 9, and 14) these impeller efficiencies were 
measured in a vaneless diffuser at a point where the radial dis- 
tance from the center of the vaneless diffuser passage to the axis 
of rotation was equal to one and one half times the impeller outlet 


paper. 


radius. 


5 Head, Advanced Engineering Department, Ingersoll-Rand Com- 
pany, Phillipsburg, N. J. Assoc. Mem. ASME. 


Journal of Engineering for Power 


The vaneless diffuser for the mixed-flow impeller tests had a 
conical rear diffuser wall with the cone half-angle (approximately 
30 deg) equal to the mean flow direction in the meridional 
plane. The front diffuser wall was contoured so that the flow 
area was reduced to 95 per cent of the diffuser inlet area in the 
first inch of length and was maintained constant at 95 per cent 
thereafter. 

For the centrifugal impeller tests (Fig. 13) impeller efficiencies 
were measured at twice the impeller outlet radius in a vaneless 
diffuser. This procedure was the same as that used in obtaining 
the original performance data for the unmodified impellers with 
which the modified impeller performance was to be compared. 
The vaneless diffuser used in the centrifugal impeller tests had a 
straight radial rear diffuser wall and an appropriate front diffuser 
wall to result in constant flow area in the radial direction. The 
over-all diffuser diameter was 25 inches. 

In the vaned-diffuser tests the efficiencies stated were efficiencies 
for the complete stage. For the 40-vane diffuser the outlet 
measuring station was at a radius of 16 in. in a vaneless section of 
constant passage height in which the flow was turned to the 
radial direction after passing through the diffuser vanes. For 
the 24-vane diffuser the outlet measuring station was located in 
the passage between blades at a radius of 13.0 in. The Mach 
number at the outlet measuring station was 0.2 for the 40-vane 
diffuser and 0.3 for the 24-vane diffuser at the design impeller 
outlet mean-line speed of 1400 feet per second. 

Detailed descriptions of instrumentation, test procedure, and 
apparatus, as well as more complete discussion of the experi- 
mental results, can be found in references [1 to 8}. 
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Force on the Impeller, and Volute Mixing 
Losses of a Radial Flow Centrifugal Pump 


A standard volute-type, radial-flow, centrifugal pump was instrumented to obtain the 
pressure distribution in the volute and also the bearing reactions from the pump hy- 
draulic force transmitted to the shaft. 
distribution was found to give a reasonable design approximation of the radial force. 
An analysis of hydraulic conditions within the volute gave pressure distributions and 
radial-force magnitudes that were comparable to those measured with certain qualitative 
interpretations about internal recirculations. 


The resultant force from the integrated pressure 


In addition, the pressure-distribution 


analysis furnished an interpretation of the effect of the volute on the pump head-capacity 


performance with corrections to the impeller head. 


The predicted head-capacity rela- 


tionship had the form of the measured pump performance. 


Introduction 


ae: IS INSUFFICIENT EVIDENCE, and no analysis, 
to enable determination of the radial force on the impeller (or 
shaft) of a centrifugal pump of any given design. Stepanoff! re- 
lated the radial force to a combination of the pump head and a 
coefficient which is a simple function of the flow rate. Presuma- 
bly, the multiplying constant in his coefficient was obtained from 
the test results on one pump. Stepanoff further stated that his 
recommended constant may not be applicable to all radial-flow 
centrifugal pumps. From tests on one pump, Binder and Knapp? 
concluded that the radial force is related to the volute pressure 
distribution together with the nonuniform momentum-rate distri- 
bution in the flow discharging around the impeller periphery. 
They had some reservations as to the reliability of their results. 


1A. J. Stepanoff, ‘‘Centrifugal and Axial Flow Pumps,” John Wiley 
& Sons, New York, N. Y., 1948. 

2 R. C. Binder and R. T. Knapp, ‘Experimental Determinations of 
the Flow Characteristics in the Volute of Centrifugal Pumps,’’ TRANs. 
ASME, vol. 58, 1936, pp. 649-661. 

Contributed by the Pumping Machinery Subcommittee of the 
Hydraulic Division and presented at the Hydraulic Conference, Ann 
Arbor, Mich., April 13-15, 1959, of THe AMERICAN Society oF 
MECHANICAL ENGINEERS. 

Nore: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 

Manuscript received at ASME Headquarters, October 
Paper No. 59—Hyd-10. 


the Society. 
27, 1958. 


at, 


Speculation on the cause of the nonuniform pressure distribu- 
tion around the impeller periphery leads to the thought that the 
pressure distribution must be due partly to the mixing of the 
fluid emitting from the impeller with the fluid already in the 
volute. If so, the fluid system may be amenable to analysis. 

Thus, the radial force problem may be considered in two parts: 
a sufficiently accurate relationship for design 


(a) Is there 


purposes between the measured pressure distribution and the 


radial force? In a practical sense, the pressure-measurement 
technique should be as simple as possible to enable determination 
of the pressure distribution quickly and conveniently. 

(b) Is there a possibility of an analytical prediction of the radial 
force? Assuming that there is a relationship between the volute 
pressure distribution and the radial force, an analysis leading to 
the volute pressure distribution would then give the force. 

There is one further consideration based upon the physical in- 
terpretation of the pressure variations in the volute. The average 
velocity at any section in the volute is, in general, different from 
the velocity of the fluid emitting from the impeller. Mixing of 
two fluid streams, with consequent accelerations or decelerations 
of the streams, produces pressure differences and also head 
losses. Thus, there is a third consideration which follows from 
(b). 

(c) Is there a relationship between the volute pressure distribu- 
tion and the pump head? 

A program was initiated to attempt to supply answers to the 





Nomenclature 


head 
volute 


volute cross-section area H, = 
volute area at tongue 

impeller-discharge area 

volute-surface area for friction 


flow-rate coefficient 


coefficients 


volute hydraulic diameter 

volute hydraulic diameter at 
tongue 

energy-loss rate 

friction factor 

gravitational constant 


pressure 


reference 


pump head 
Euler impeller head 
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loss due to 
Stodola impeller head 


volute-area variation coefficient 


volute-hydraulic-diameter varia- 
tion coefficient 
number of impeller vanes 


pressure at 

discharge section 
volute-section perimeter 
flow rate at a volute section 
impeller flow rate 


flow rate through Ay 

impeller radius 

impeller peripheral velocity 

velocity at a volute section 

radial velocity of fluid leaving 
impeller 

tangential velocity of fluid leav- 
ing impeller 

Stodola corrected V, 

volute-area variation angle 

impeller-vane discharge angle 

fluid density 

angle from volute tongue in 
direction of impeller rotation 

shear stress 
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mixing in 


volute- 








9 VANES 


23.5° 


Fig. 1 


foregoing three considerations. Pressure distributions and radial 
forces were measured in a radial-flow, centrifugal pump with a 
specific speed of L000. A force analysis and a mixing-loss analysis 
were made for the conditions in the volute to show the possibility 
of predictions of the pressure distribution, of the radial force, and 
of the head loss due to mixing. The program was not developed 
in an elaborate fashion, but rather on the naive basis that reasona- 


ble answers may be possible from relatively simple techniques. 


Experimental Program 


A 3-in. double-bearing, horizontal-shaft, single-suction, single- 
stage, radial-flow centrifugal pump of normal industrial manu- 
facture was used for the test program. The pump was selected 
from the supply of service pumps available to the laboratory 
mainly on the basis that it looked like a pump of reasonable de- 
sign and could be instrumented as needed without too much 
difficulty. 

Essential details of the pump impeller and of the volute are 
shown in Figs. 1 and 2. The impeller was modified by fastening 
rings to the shroud periphery to reduce the clearance to a nominal 
‘/ig in. between the shroud and the casing. There were two 
reasons for this modification; one was to obtain a definable volute 
cross-sectional area, and the other was to minimize flow between 
the volute and the clearance passage formed by the impeller 
shrouds and the casing. 

The arrangement of the volute pressure taps is shown in Fig. 2. 
Relative locations of the pressure taps around the volute are 
The original intent in the test program was to 
Pres- 


shown in Fig. 3. 
use the simplest possible pressure measurement technique. 
sure taps were then located as shown in the narrow section of the 
volute such that the pressure opening in the volute wall was es- 
sentially parallel to the main fluid flow issuing from the impeller. 
One series of tests was made with the simple side-wall taps. A 
second series of tests was made with the static-probe extensions as 
shown at location 5, Fig. 2, on taps 1 through 10. The nominal 
location of the original side-wall pressure taps was '/2 in. from the 
impeller periphery and the openings in the static probes were a 
nominal !/, in. from the impeller periphery. 

The radial force taken by the shaft as transmitted from the im- 
peller was measured at the bearing mounts. A section was re- 
moved from the bearing brackets to separate completely the 
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Fig. 2 Test-pump typical volute secti Numbers refer to pressure 
taps. A typical static-probe arrangement is shown at section 5. 





bearing housings from the pump casing. Each bearing was then 
mounted on a strain-beam spider, Fig. 4. These were four rectan- 
gular steel bars welded to a center ring to form four legs at 90-deg 
intervals. The beam ends were fastened rigidly to a special sup- 
port on the bearing cap, and the center ring was fastened rigjdly 
to a stand attached to the test bedplate. When assembled the 
impeller and shaft were supported only by the strain beams. 
The beams were oriented with the plane of the spider in a horizon- 
tal plane parallel to the pump shaft, and with one pair of beams 
parallel to the pump shaft and the other perpendicular to the 
axis of the pump shaft. The plane of the spider was approxi- 
mately 3 in. below the center line of the pump shaft. 
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Fig. 3 Location of pressure taps 


DIRECTION OF 
PUMP DISCHARGE 


1" 
= DIAM 
siz 0 


8-TYPE A-7 
-4 

SR GAGES _ DIRECTION OF 
c 1. ae 

6 \ i | 


PUMP SHAFT 
CENTER BOLTED 
TO BEARING STAND 


NO'S ON LEGS 
REFER TO TOP 
GAGES 


ENOS BOLTED 
TO BEARING MOUNT 


NOS BESIDE LEGS 
REFER TO BOTTOM 
GAGES 


EXCI 


Fig. 4 Strain-beam and bridge circuits for bearing support forces 


SR-4 strain gages were fixed to each strain beam, one on top and 
one at the bottom, both near the center-ring attachment. The 
strain gages were wired to a multicontact, multiposition switch 
to form Wheatstone-bridge circuits with an SR-4 strain indicator. 
Gages in the same location on the two spiders, one at each bearing, 
were wired in series. 

The strain-beam arrangement and the circuit arrangement per- 
mitted isolation of three shaft-force components; vertical, hori- 
zontal perpendicular to the shaft, and axial parallel to the shaft. 
Each foree-component-magnitude response on the strain indica- 
tor was independent of the magnitude, or changes in magnitude, of 
the other two force components. 

Some deflection of the complete shaft and bearing-cap assembly 
was expected under the pump test-loading conditions. To prevent 
rubbing at the wear rings, all running clearances in the pump were 
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enlarged to approximately 4/6 in. Internal flow recirculations 
due to these enlarged clearances did not affect the test results of 
forces since the desired comparison was that between the meas- 
ured pressure distribution and the simultaneously measured shaft 
force. 

The normal shaft-seal packing was removed and one flexible 
chevron ring installed in each stuffing box. Some shaft-seal leak- 
age was expected with this packing arrangement, but the flexible 
seal was needed to eliminate any force on the shaft except that of 
the hydraulic loading on the impeller. 

The pump was installed on a laboratory test stand with piping 
and instrumentation for a normal pump test. The net pumped 
flow rate was measured using a calibrated orifice with both mer- 
cury-water and water-air differential manometers. The pump 
head was measured by a mercury-water differential manometer 
connected to static pressure taps in the suction line and in the 
discharge line. The pump was driven by a variable-speed, d-c, 
cradled motor through a pin-and-leather-strap flexible coupling. 
The coupling arrangement did not impose any load on the shaft 
other than the direct input torque. All volute pressure taps were 
manifolded to one side of a mercury-water manometer and to an 
inverted water-air manometer. A pressure tap in the pump-dis- 
charge nozzle was used as a reference for the volute differential 
heads. 

The strain-beam and bridge-circuit response to load was ar- 
ranged for direct calibration with the pump assembled on the test 


stand. Flexible cords were attached to the bearing housings 


HORIZONTAL 


STRAIN BEAM MOUNT” LOAD 


Fig. 5 Arrangement of external loading devices on pump 
bearings 


and led over dry ball bearings which were mounted on rigid shafts 
independently supported from the pump, Fig. 5. Weights ap- 
plied to the cords, or placed on the bearing housings for vertical 
downward loading, gave direct magnitudes of the applied forces. 

Two test sequences were run as follows: 

Test I. Pressures from the volute side-wall pressure taps, and 
shaft forces from the calibration of the strain-beam system. 

Test II. Pressures from the static-pressure tubes, and shaft 
forces from the externally applied loads needed to bring the strain 
response to the no-load or null position. 

For both the test sequences, the no-load strain response was ob- 
tained with the pump running dry at the test speed of 1000 
rpm. The strain-system response to horizontal, vertical, and 
axial loads was obtained. The mutual independence of the re- 
sponse to one force component when the other two were varied 
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was checked for various combinations of loading. In addition, the 
averaging effect of the circuit of the two strain beam mounts, one 
at each bearing, was checked by loading the bearing blocks sepa- 
rately and in combination. The complete array of calibrations 
and checks showed that each force component was measurable 
within 1 Ib at all test loadings. Repeated checks were made 
during the tests of the strain calibrations and the null indications 
of the strain-response system. 

After calibration, water was admitted to the pump, the pump 
was brought up to the nominal test speed of 1000 rpm and data 
were taken for the capacity, head, power, volute pressure dis- 
tribution, and the three components of force at the bearing 
mounts. In both tests, the operating conditions were changed 
gradually to prevent sudden loads which would cause excessive 
deflections of the impeller and shaft, and possible rubbing. 


Experimental Results 

All test results were obtained at pump speeds between 990 
and 1020 rpm. At each capacity the speed varied not more than 
2 rpm during the test. Centrifugal-pump-model laws were ap- 
plied to the reduced data to obtain the final results for the speed 
at 1000 rpm. 

The pump performance is shown in Fig. 6 for the two test se- 
quences. Differences in the head-capacity results are attributed 
to the differences in internal-clearance distributions within the 
pump. The high efficiencies of Test II are not correct due to an 
undetected mechanical fault in the dynamometer scale which 
oceurred between Tests I and II. The efficiencies are reported 
not for the actual magnitudes but to show the pump capacity in 
the range of maximum efficiency. 

Typical static-head distributions in the volute are shown in 
Fig. 7. The reported static head is the manometer differential 
as read in terms of the pumped fluid and is not corrected for the 


elevation difference between taps. 
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Fig. 6 Pump performance at 1000 rpm 
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Fig. 7 Typical static-head distribution in the volute 


The net forces on the impeller as evaluated from the integrated 
head distributions around the volute are shown in Figs. 8 and 9. 
The total width of 11%/32 inch to the outside of the added rings on 
the impeller periphery was used as the area over which the pres- 
sure acted. Also shown is the angle of the action line of the re- 
sultant force. 

The measured force from the strain-system response is given 
also in Figs. 8 and 9 for the two tests. Since the ‘‘zero’’ or null 
indication of the strain response was with a dry pump, the meas- 
ured force indication was corrected by the inclusion of the buoy- 
ancy of the impeller. Thus the results are directly comparable for 
the force static head from the distribution and for the measured 
forces at the bearings. The consistency of the data is shown by 
the scatter of the results at each flow rate. 

In the small force range (around 150 gpm), a few of the angular 
location results are not shown. The resulting angles in this range 
were not considered representative since the forces were obtained 
with measured accuracies of 1 lb and variations within 1 lb of one 
force component would change the angle by as much as 100 deg. 

Results of the two tests show, in general, agreement of the 
bearing reactions in magnitude and direction. The force direc- 
tion evaluated from the head distribution with the static probes 
shows better agreement with the bearing-reaction direction than 
does the force direction from the head distribution with the 
volute side-wall taps. 

The general agreement between the forces obtained from the 
static-head distribution and the measured bearing forces, while 
not perfect, is of practical value since the results show that there is 
a correlation, both qualitatively and quantitatively to a good first 
approximation, between the volute static-head distribution and 
the force transmitted to the pump shaft. 

There are many reasons why the correlation is not perfect. 
Binder and Knapp showed nonuniform velocity distributions in 
the discharge around the periphery of the impeller depending upon 
the position in reference to the volute tongue. Thus there is 
the possibility of an appreciable radial-force component due to the 
nonuniform momentum rate as a function of angular location. 


1960 / 139 


APRIL 





uw 
2) 





] TEST I 
© SVOLUTE PRESSURE 
_.@ * BEARING REACTION 


B= 
oO 





w 
oO 


FORCE - POUNDS 
n 
°o 


3 








o 


ono Ls 
Oo O° 








ANGLE FROM TONGU 
nm 
oO 


@ 
oO 





100 150 
CAPACITY - GPM 


Fig. 8 Radial force in magnitude and direction from volute side-wall 
taps for static head and from strain-beam calibrations for bearing reac- 
tions. Closed and open circles represent results from data taken simul- 
taneously. Other points were taken from separate test sequences. 
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The pressure as measured for both Tests I and II was a single 
pressure at a volute position. A pressure variation may exist 
across the impeller width. The volute side-wall pressure taps of 
Test I were placed in the narrow section of the volute to attempt 
to obtain representative static pressures in a simple fashion. Ad- 
mittedly, this location does not give the true pressure at the im- 
peller periphery. There may be other forces on the impeller from 
fluid pressure and shear distributions along the impeller shrouds 
and at the wear rings. 

But the general agreement between the radial force from the 
volute static-head distribution and the measured radial force is 
encouraging. If acceptable, then the first of the original three 
considerations is answered in the affirmative and the second, that 
of an analytical prediction of the pressure distribution and the 
radial force, can be attempted. 


Pressure Distribution and Radial-Force Analysis 

The volute is a collection chamber for the flow from the impel- 
Each tangential section of the impeller discharges a portion 
A force 


ler. 
of the total flow into that which is already in the volute. 
balance on a volute element should give the pressure distribution. 

Consider the volute section in Fig. 10. For simplicity, and also 
through ignorance of the pressure and velocity distributions in 
the radial direction, the following assumptions are made: 


1 The pressure is constant over the area A. 
The velocity is constant over the area A. 

3 The volute radial width is small compared to the impeller 
radius such that all moment arms about the central axis may be 
taken as the impeller radius. 

4 The discharge from the impeller is uniform over the width of 
the impeller, and also over the circumference of the impeller. 


5 The fluid is incompressible. 
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A moment balance about the central axis and in the direction 8 


1S 


| 


~ 


d 
Apr — (A + dA)(p + dp)r + (, + ’) dA, (sin @)) 


— TdA, (cosa)r = p(Q + dQ\V +dV)r 

— pQVr — pV, dQr 
Equation (1), with second-order terms neglected, reduces to 

— Adp — pdA + psina dA, — tT cosa dA, 
= pVdQ + pQdV — pV, dQ 
Also sina dA, = dA (3) 
cos a & 1.0 (for reasonable volute-area variations). And, from 

continuity, 

dQ Q 
A A? 


dA (4) 


The boundary friction represented by the shear stress may be 
put in the form containing the friction factor. Also, introducing 
the perimeter P of the shear area, and the hydraulic diameter 
D = 4A/P 

dA, fpV?Prd0 — fpV*rd0 
T = ee = 


x 2D 


i. (5) 


With equations (3), (4), and (5) substituted in equation (2), 
the differential equation for the pressure distribution is found. 


, @ (? V, fQ*rdb : 
—dp =p (2 av? = 404 eal dQ + 2A*D (6) 


In order to integrate equation (6), the variation of the volute 
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The test pump had a 
measured linear variation of area in the direction 6, 


A = Ay + Ki6 (7) 


in which Ao is the clearance area between the tip of the tongue and 


cross-section area with 8 must be known. 


the impeller periphery. 

When Ap is included, then the flow rate at any section of the 
volute will be the combination of that which discharges from the 
impeller from the tongue to the position 6 and the flow through the 
area Ag, 

(8) 
where Kz = (9) 
2r 

Q; = total flow through the impeller 

The hydraulic diameter of the test pump was found also to 
have a linear variation with 6, 

D=D,+ K;6 
where Do is the hydraulic diameter at Ao. 

When equations (7), (8), and (10) are substituted appropriately 


in equation (6) to give the differential form of the pressure as a 
function of the angular location and with the condition 


(10) 


p = p, at 0 = 2r 


the integration yields 


| ene |v K» (“)' ( Ay + KO 
= —- I 
p sh? '\ Ao + 20K 


_ (Kio — xan aha ilnria’ 
| 2K? (Ao + 24Ki)? (Ao + Ki8)? 


rf ° ] l Ao + K\6 
+ C; = = = bab Cs hie 
2 Ao + K,6 Ao + 2rK, Ao a 2rK, 


Dy. + K30 
- cxin( + 4 (11) 
Do + 2rK; 


(AoK2 — QoK1)? 
K\*( Dok, = AoK3) 


where i= 


Cy = 
K,?(DoKi — AoKs)* 


y 


.! 


(QoKs — DoK)* 
K3( Dok, — A oKs 2 
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The first two terms of equation (11) are from the mixing proc- 
ess. The remainder of the terms are the frictional effects. Con- 
sider only the mixing terms. The second term is always negative. 
If the flow rate through the impeller is small (K2 is small), then 
the first term is negative. For large flow rates the first term be- 
comes positive. The pressure difference as noted may be negative 
for low pump capacities and positive for high capacities. This 
qualitative trend agrees with that found from the test results, 
Fig. 7. 

For a quantitative evaluation numerical magnitudes are needed. 
The factors Ao, Ki, r, Do, and K; are fixed by the pump geometry. 

The tangential velocity of the fluid leaving the impeller V, is 
obtained from the impeller-discharge velocity pattern. Osborne 
and Morelli? have shown that the Stodola correction to the Euler 
head is a reasonable approximation to the impeller-discharge 
conditions. The inference is that the velocity pattern at dis- 
charge is taken from the Stodola correction. The Euler head for 
the impeller is 


U2 Q,Us 


H,; =— -- cot Be (12) 
q Ag 


The Stodola corrected Euler head for the impeller is 


i 2 U2? JU 
Hs = (: -_— fs) = cot Bs 


(13) 
n g Ag 


Then 
SS own Bs (14) 
A; 

The remaining variable that must be investigated is the mag- 
nitude of the flow rate Qo through the clearance between the tip 
of the volute tongue and the impeller periphery. The clearance 
area Ao of the test pump was small compared to the final dis- 
charge area at the volute end. Since no measurements were 
made of the details of the flow at the tongue clearance, and since, 
for a small area Ao, the flow rate is small unless the velocity 
through the area is large, the clearance flow was taken as zero as 
a reasonable approximation. 


Qk 2(QoKs eke Do Ke) oe Aok2*( KiDo = AoK;) + DoKi Kf Aok2 ~ QoKi) 


’W. C. Osborne and D. A. Morelli, ‘‘“Measured Performance of 
Pump Impellers,’” ASME Paper No. 50—A-90. 
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Fig. 11 Predicted static head variation in test-pump volute; mixing and 
friction included 


The friction-factor magnitude was taken arbitrarily at the 
nominal value of 0.02. 
The measured and assumed conditions for the test pump were 


22.5 deg 

0.100 sq ft 

0.423 ft 

9 vanes 

Q;/2m cfs per radian 


Ay = 0.00244 sq ft 
K, = 0.00346 sq ft per rad 
D, = 0.107 ft 
K; = 0.010 ft per rad 
= 0.02 (assumed 


= 44.3 fps (at 1000 rpm 0 (assumed) 


The numerical evaluation of the pressure distribution is shown 
in Fig. 11 for selected values for the flow rate through the im- 
peller. The relative influence of the volute frictional effect is 
shown for the one condition of flow at 300 gpm. The frictional 
effect is appreciable and has been included at all flow rates. 

In general, the trends of the predicted pressure variations with 
position around the volute and with flow rate agree with the 
measured results. Direct comparisons are not possible due to one 
major factor—that of internal recirculation from the volute 
through the wear rings and into the suction of the impeller. Re- 
call that the test-pump wear-ring clearances were enlarged to 
nominal magnitudes of '/is in. Internal recirculation rates may 
be appreciable in terms of the net pump flow rate. 

Kadial forces from the predicted static-head distributions in 
the volute were evaluated by a numerical summation in the same 
manner as the technique used with the measured static-head dis- 
tributions. Direct analysis from equation (11) is possible, but 
the forms to integrate are rather involved and the complex final 
expressions are as difficult to evaluate as the numerical integra- 
tion. 

The resultant predicted force in magnitude and direction as a 
function of the gross flow rate through the impeller is shown in 
Fig. 12. At first glance, the prediction does not seem to have any 
validity. But the flow-rate variable must be interpreted in terms 
of the test conditions and of the predicted conditions. In the test 
there was recirculation of appreciable magnitude from the volute 
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Fig. 12 Predicted radial forces in magnitude and location for gross flow 
rates through pump impeller 


to the impeller suction. If the predicted curves of Fig. 12 are 
shifted to lower net capacities corresponding more to the meas- 
ured pump capacity, then better agreement results. The recircu- 
lation is dependent upon the pump head, or pressure difference 
between the volute and the suction, upon the pressure distribu- 
tion in the volute, and upon the wear-ring clearances. In terms 
of the head alone, the shift of the predicted curves to lower net 
capacities should be greatest at the low flows and least at the high 
flows. 

The predicted force magnitude of 45 lb at the gross capacity of 
125 gpm is in the range of the measured force at shut-off of the test 
pump. Also, the predicted force of 45 lb at the gross capacity of 
300 gpm is in the range of the measured force at the net pump 
capacity of 300 gpm. This flow rate occurs at approximately zero 
measured head—the condition approaching minimum recircula- 
tion leakage and hence an expected gross and net capacity agree- 
No attempt has been made to adjust the gross impeller 
More information is needed 


ment. 
capacity to the net pump capacity. 
on the actual recirculation magnitudes. 

The predicted static-head distribution and the predicted re- 
sultant force show qualitative agreement with the measured 
conditions, and appear to show quantitative agreement with 
some judicious interpretations of circulation and leakage condi- 
tions within the pump. 


Head-Capacity Analysis 

If the static-head-distribution and radial-force analysis are ac- 
ceptable, then the analytical system may be investigated further 
to determine the effect of the postulated mixing analysis upon 
the head-capacity relationship of the pump—the third considera- 
tion originally stated. 

The impeller head has been stated already in equation (13). 
The mechanical energy-loss rate in the mixing process in the 
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volute results from a mechanical-energy balance on an element 
of the volute as in Fig. 10. For the element, the mechanical- 


energy rate loss is the difference between the energy-rate input 
and the energy-rate output, or, 


y2 | d (V2+ V,2) 
dE, = ol +e-> | +0 (» + P) 4p Bris | 


V + dv)? 
~ (9 + 40)| (w +a) + 0 rs *] (15) 


Expanding and reducing equation (15) gives 


V2 
dE,, = dQ — Qdp — pQVaV — p> dQ (16) 


(V2 + V,2) 
9 

Substitution for the differential pressure from equation (6) 
without the friction term, for the appropriate velocity form from 
equation (4), for the volute-area variation from equation (7), and 
for the volute flow rate from equation (8), and integrating equa- 
tion (16) over the complete volute from 8 = 0 to 6 = 2m yields 
the total energy-loss rate due to mixing. 


2rk, 
x | ve . ve + 


(QoK: — K2eAo)? 
K,*Aol( Ao + 27K) 
K, V, ) Ay + 29K, 
In 


+ (KiQo — K2Ao) ( sk te 
17K 3 TK 2 
K,? 


9 K;? 


The head loss in the volute is the total mechanical-energy rate 
loss divided by the weight rate of flow, or 


En E,, 


_ = ” ( 18 ) 
Qipg  2rK2pg 


H,, = 
Thus the pump head is the impeller head from equation (13) with 
the volute head loss from equation (18) subtracted 


H = H, — H, (19) 


In addition, other losses, such as friction in the impeller and in 
the volute, should be included in the final complete representation 
of the pump head. In order to focus attention on the volute 
mixing-loss terms, frictional effects will not be included in the 
present discussion. 

As a first step in the explanation of the significance of equation 
(19), take the condition K. = 0 (or Q; = 0), Ao is small but finite, 
Qo = 0, and for a pump with a large number of impeller vanes such 
that the Stodola correction to the impeller head reduces to one. 
The pump head becomes 


U? 


V, is zero for this condition. 
H = (20) 


g 
Since, at Q; = 0, V, = U2, then 
U2? 


2g 


H = 


(21 ) 
or the shut-off head is one half the impeller head—a value which 
is in general agreement with test values of radial-flow centrifugal 
pumps. 

The evaluation of equation (19) for the test pump is shown in 
Fig. 13 together with the Euler head for the impeller, the Stodola 
corrected Euler head, and the test-pump head as measured. 
Again, the differences of the measured and predicted head-capac- 
ity values are subject to the approximations in the analysis, to 
the omission of frictional effects in this case, and to the unspeci- 
fied recirculation flow rate between the volute and the impeller 
suction. 
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Fig. 13 Pump head-capacity relationship showing effect of volute 
mixing losses 


Friction losses will lower the predicted head, and the recircula- 
tion will shift the predicted curve to smaller values of the net 
pump flow rate. Both refinements move the predicted magni- 
tudes in the proper direction to obtain better agreement between 
the predicted and test head capacity. At least, as far as this 
evaluation is carried, the shape of the predicted pump head- 
capacity curve is certainly an improvement over the acceptance 
of the impeller head as the total pump head. A partial answer to 
(c) is thus proposed. 


Summary 


The radial force in the test centrifugal pump as obtained from 
the measured pressure distribution in the volute is a reasonable 
design approximation to the measured radial force. The pressure 
measurement need not be elaborate since results for the test pump 
were obtained with simple volute side-wall pressure taps. 

The analysis of the mixing process in the volute, while subject 
to many simplifications, resulted in volute pressure distributions 
and radial forces which were on the order of the measured radial 
forces. 

The analysis also yielded an explanation of the mechanical- 
energy losses due to mixing in the volute and a consequent head- 
capacity modification to the impeller head that shows the in- 
fluence of the volute in predicting the pump performance. 


DISCUSSION 
F. C. Gilman‘ 


Worthington Corporation has been much interested in de- 
terminations of the radial force acting on impellers of centrifugal 
pumps of various configurations and over the full capacity range 
of operation. Some recent work has been reported by Agostinelli, 
Nobles, and Mockridge.’ We are gratified that the authors 
have been investigating this phenomenon and have now published 
their findings. 

The method employed for measuring bearing reactions is in- 
genious and simple. It is unfortunate that the authors chose to 

‘Senior Research Engineer, Research and Development Section, 


Worthington Corporation, Harrison, N. J. Mem. ASME. 
5 Published in this issue, pp. 120-126. 
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modify their centrifugal pump so drastically by adding cylindrical 
surfaces at the impeller periphery and by increasing the wear ring 
clearances to the extent that upwards of 30 per cent of normal 
capacity was probably being bypassed to suction. While it may 
be argued that volute pressure distribution is dependent mainly 
on the net volute flow it is this discusser’s experience that im- 
peller discharge flow exerts an important influence. It would be 
interesting if this test could be repeated using a face seal or some 
other system such as displacement of the bearings to permit 
operating with normal clearance. 

The addition of the rings at the impeller periphery has made 
the test pump nonrepresentative of normal commercial pumps of 
American design for two reasons whose effects are cumulative. 
First, by closing off communication between the volute region 
and the regions flanking the impeller the normal tendency 
toward equalization by flow across these regions was minimized. 
Secondly, the frictional drag of the added rings could be expected 
to further increase the static pressure differential at low flows. 
An over-all change in differential pressure equal to normal pump 
head is unusually high for 1000 specific speed. 

The authors present the data of Figs. 8 and 9 in support of the 
thesis that volute static pressure distribution alone is sufficient 
to predict impeller radial forces. At the maximum force condi- 
tion both figures show magnitude discrepancies of about 25 per 
cent and direction errors amounting to approximately 30 degrees. 
At higher flows, although magnitudes seem to agree more nearly, 
angle errors are of the order of 70 degrees. The addition of the 
rings at the impeller periphery could be expected to increase the 
static pressure contribution without directly affecting momen- 
tum forces. Of course if they altered the static pressure dis- 
tribution as previously suggested, this would affect the momen- 
tum distribution. It is significant that even though the test 
pump was altered to favor the influence of static pressure, a con- 
siderable discrepancy still remained. 

The analysis of pressure distribution within the volute, results 
of which are given in Fig. 11, shows completely reversed effects in 
relation to test at normal flow and zero flow. Only at maximum 
flow is there good correspondence. It is easy to understand why 
the assumptions upon which the analysis is based entirely fail to ex- 
plain the pump action at zero net flow. With uniformly zero flow 
from all parts of the impeller, uniformly zero velocity in the volute 
and no friction drag ascribed to the periphery of the impeller, 
there is no mechanism to provide a differential static head in the 
volute. In fact there would be no reason to expect a discharge 
head higher than U,?/2g. 

The assumptions of the analysis more nearly match actual 
pump conditions at normal flow. Why is it then that, while the 
pump designer has succeeded in producing nearly uniform volute 
static pressure at the normal flow condition, the present analysis 
indicates a large differential static head corresponding to the 
tested distribution at zero flow? 


Authors’ Closure 


Mr. Gilman objects to the modifications of the tested pump 
which made the pump nonstandard and to the analysis which 
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does not predict exactly the conditions in the pump volute and 
the consequences of these conditions in terms of the radial forces. 
The authors agree to both objections. However, the reported 
test and analysis were not made to give all the answers required 
for industrial pump design. Rather, the work was done on the 
selected system to see what possibilities existed to relate the 
volute pressure distribution to the radial force and to apply a 
simplified analysis as an initial attempt toward possible predic- 
tion of volute conditions. The pump modifications were neces- 
sary to define volute areas for inclusion in the analysis and to 
furnish a workable test arrangement in this initial investigation. 

While the impeller rings tended to produce larger than ‘‘nor- 
mal’’ magnitudes of the nonuniform pressure distribution in the 
test pump by “closing off the communication between the volute 
region and the regions flanking the impeller,’ there is another 
difference between the test pump and industrial pumps that 
contributes to the “high”? change in differential head. The 
volute tongue in the test pump was brought down to the impeller 
to close off, as much as possible, the recirculation past the tongue 
in the area between the tongue and impeller. With the small, 
but finite, area at the tongue there was a recirculation through 
this area other than zero as assumed for convenience in the numer- 
ical evaluation of the analysis. The recirculation flow term in 
equation (11) modifies the predictions of Fig. 11 to show a pressure 
variation in the volute other than zero at zero impeller flow. The 
shift at low pump-flow rate is in the direction of zero negative 
differential static heads and toward better agreement with the 
pressure distributions found in the tests. 

The agreement in magnitude and direction of the radial forces 
between those measured at the bearings and those computed 
from the pressure distribution is not perfect. But most engi- 
neering design is not perfect in the sense that exact magnitudes of 
all design criteria are not known at the blueprint stage. If the 
radial force is a critical design condition, then relatively simple 
tests to determine pressure distributions in the volute may furnish 
approximate magnitudes of radial forces that are better than no 
information. 

The analysis from which Fig. 11 was developed is based on 
gross flow through the impeller. The point is well taken that the 
prediction is not valid for zero gross flow (and for other flows as 
well) since there is a shear drag from the impeller acting on the 
fluid in the volute, in addition to the previously mentioned recir- 
culation at the tongue. This drag would tend to displace all 
curves downward toward higher negative heads in terms of the 
sign convention established for the pressure distribution. Other 
simplifications in the analysis also exaggerate the predicted pres- 
sure distribution as compared to the measured pressure distribu- 
tion. Then there is the wearing ring leakage factor which is 
unknown and makes direct comparison qualitative only between 
predictions based on impeller gross flow and measurements re- 
lated to net pump flow. 

Even though most of the information presented in terms of this 
one pump may be qualitative only for design purposes, there 
seems to be sufficient agreement to stimulate further work which 
may lead to the ultimate objective of all engineering practice— 
the ability to design completely from analytical formulations. 
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Relationship of Coal-Ash Viscosity 
to Chemical Composition 


This paper discusses information on coal-ash behavior in the temperature range suitable 
for liquid ash removal and shows how this is related to boiler design and fuel characterts- 
tics. A description of a method for predicting slag viscosity from ash analysis and 
fusion temperatures is given. 


Also, a method is shown for predicting the flux require- 


ments to obtain desired slag flow characteristics from an unsatisfactory fuel. 


Woes THE GROWTH of the slag tap boiler, the ability 
to predict accurately the fluid characteristics of the ash by 
laboratory means becomes increasingly important. In the past, 
because of the smaller units and the more abundant sources of 
coal, a miscalculation of coal suitability could usually be cor- 
rected by changing fuel sources. However, many modern plants 
because of location or long term contracts are often committed 
in advance to burn a specified coal. 

Furthermore, the flow characteristics of an ash in relation to 
furnace temperature have considerable bearing on the design and 
selection of firing equipment. 


Flow Characteristics 


Ash fusion temperatures have long been used as a guide in pre- 
dicting ash flow characteristics, but they have a shortcoming in 
that they do not tell how fluid the ash will be. 
having the same ash softening temperatures may have widely 
different flow characteristics, as illustrated in Fig. 1. The degree 
of fluidity is commonly expressed in terms of its viscosity or re- 


Two coal ashes 


sistance to flow at a stated temperature. 

Although the viscosity of a slag can be measured, this operation 
requires special high temperature equipment of a type not found 
in most coal laboratories. However, most commercial labora- 
tories do have at their disposal means for determining the chemi- 
cal composition of an ash. This can be done either by wet 
analysis or by more recently developed spectrographic means 
{1}.! 
istics from chemical analysis provides a very useful tool in boiler 


Hence an accurate means of predicting ash flow character- 
design. It also provides information on how the flow character- 
istics of an ash may be altered by fluxing to obtain the desired 
characteristics. 


Past Work 


Ash as it exists in coal is composed chiefly of compounds of 
silicon, aluminum, iron, and calcium, with smaller amounts of 
titanium, magnesium, sodium, potassium, and other trace ele- 
ments. 

Many attempts have been made in the boiler and in the steel 
industries to correlate the flow characteristics of slags with their 
chemical composition. Probably the most extensive investiga- 
tion of the relationship of slag flow to chemical composition for 

1 Numbers in brackets designate References at end of paper. 

Contributed by the Fuels Division of THe AMERICAN Society oF 
MECHANICAL ENGINEERS and presented at the joint AIME-ASME 
Solid Fuels Conference, Cincinnati, Ohio, October 26-29, 1959. 

Note: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 
the Society. Manuscript received at ASME Headquarters, July 20, 
1959. Paper No. 59—Fu-4. 
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coal-ash slags that has been published to date was the one made 
by the Bureau of Mines [2, 3, 4, 5]. However, the scope of this 
investigation covered a limited range of ash composition. 


Need for Additional Information 


As new sources of coal for steam generation have come into 
use and as new and improved methods of firing and furnace de- 
signs have been developed the need for more detailed knowledge 
of ash characteristics has become increasingly important to 
operators and designers of slag tap boilers. To meet this need it 
has been necessary to expand and intensify our studies of ash 
behavior. These studies include ash characteristics that influence 
(a) fluidity in the tapping range, (b) boiler fouling, and (c) metal 
corrosion. 

This paper discusses information on coal-ash behavior in the 
temperature range suitable for liquid ash removal and shows how 
this is related to boiler design and fuel characteristics. A de- 
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Fig. 1 Plot of ash softening temperatures in a reducing atmosphere 
(hemispherical point) versus temperature for a viscosity of 50, 100, and 
250 poises 
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scription of a method for predicting slag viscosity from ash 
analysis and ash fusion temperatures is shown. 


Scope of Investigation 

The scope of this investigation has included a review of availa- 
ble information from previous data and from published articles 
along with a co-ordinated two-phase test program: 


(a) Test firing a number of fuels to observe slag tapping charac- 
teristics in laboratory and commercial units. 

(b) Determining slag viscosity and chemical analysis of these 
This work included ashes altered by the 
addition of chemicals. Data on these slags are tabulated in 
Table 1. A description of the viscometer used is included in the 


Appendix. 


and many other ashes. 


The results of these studies provide a practical basis for de- 
lineating the relationship that exists between (a) slag viscosity 
and furnace temperature, and (b) chemical analysis and the vis- 
cosity of a slag. 

The influence of these relationships in predicting the slag tap- 


Table 1 


CHEMICAL ARALYSIS 
T102 Cao ID 


Fe70; 
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ASH FUSION 
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ping limitations of a power boiler will be discussed in the follow- 
ing. 


1 Relationship Between Slag Viscosity and Furnace 
Temperature 


In order to predict the suitability of a coal ash for liquid slag 
removal, it is necessary to know its flow characteristics in relation 
to temperature. Fluidity of a slag is usually described by giving 
its viscosity in poises at a stated temperature. In order to es- 
tablish these limits, tests were made on a number of units where 
the temperature of the slag was obtained by optical pyrometer, 
and its viscosity was determined by laboratory tests of a slag 
sample. Results of these tests show that the slag viscosity as its 
taps from a boiler may vary considerably, however, it is normally 
within the range of 50 to 100 poises with 250 poises as an upper 
limit. Data reported by the Bureau cf Mines serve to confirm 
these findings [3]. There are no commonly used liquids having 
viscosities in these ranges; however, to give some visual concept 
the viscosity of 


Data on coals used for viscosity determinations 


MEASURED SLAG VISCOS’TY 
TEMP °F TEMP °F FOR POLSES 
rr $0 100 250 
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0.01 poises 
1.00 poises 
15.00 poises 
120.00 poises 
1300 .00 poises 


Water at 68 F is 
SAE #10 oil at 68 F is 
Glycerine at 68 F is 
32 F is 
0 F is 

The boiler designer has some control over setting the limits for 
satisfactory slag tapping from a furnace in the arrangement of 
the water-cooled surface and in directing the flow of the furnace 
gases. However, for a given furnace design and slag composition, 
the governing factor in determining whether the slag will tap 
satisfactorily from the furnace is the furnace temperature. Fur- 
nace temperature at the tapping zone is the adiabatic or theoreti- 
cal? furnace temperature minus heat losses up to this point. Once 
the fuel analysis, excess air, and air preheat are known, the 
theoretical furnace temperature can be calculated. Calculating 
heat loss to the enclosure is quite involved; however, for predict- 
ing limits on slag tapping ability it is sufficient to know the rela- 
tionship between furnace surface and total heat release. This is 
illustrated in Fig. 2 in which curves are plotted for a typical boiler 
design to show the upper limits on satisfactory liquid slag tapping. 
This figure shows that, for a theoretical furnace temperature of 
3900 F, a slag having a viscosity of 250 poises at 2600 F would tap 
satisfactorily at rated capacity. However, if the rating were 
dropped, that is, if its heat release rate were lowered (thereby in- 
creasing the ratio of surface to heat release), slag flow would 
cease. This does not mean that such a boiler cannot operate at 
50 per cent rating, but only that it would not tap liquid ash until 
rating is increased. Further study of this figure shows that ash 
viscosity limitations are lower on a small unit of given configura- 
tion because of the higher surface to heat release ratio. 

To calculate the theoretical furnace temperature it is necessary 
to know three things: (a) The fuel analysis, (b) the amount of 
excess air, and (c) the degree of air preheat. For example, at 600 F 
air preheat and 10 per cent excess air, the normal range is from 
4100 F for a high rank eastern coal of low ash and moisture con- 
tent down to 3500 F for a midwestern coal of 20 per cent moisture 
and 20 per centash. Furthermore, once the tapping limitations for 
a unit have been established on one coal, Fig. 2 can also be used to 
predict limits on other fuels by locating the limiting point for 
temperature and viscosity and drawing a curve parallel to those 


shown on this figure. 


2 Relationship of Ash Viscosity to Chemical Composition 

The chemistry of coal ash is quite complex due to the number of 
elements involved and the reactions that take place. Several sim- 
plified methods of predicting slag behavior based on chemical 
analysis of the ash have been proposed. 

The two methods to be described are not original with the 
authors but are modifications of previously proposed methods 
with additional data to support their validity. The first of these 
is a modification to the method proposed by the Bureau of 
Mines [4], which permits plotting a reasonable approximation of 
the viscosity-temperature relationship in the tapping range. 

The second method shown relates ash viscosity to the base-to- 
acid ratio of the chemical constituents. This method is more 
applicable to predicting the effect of flux additions or the result of 
blending two or more fuels. 


Modified Bureau of Mines Method 


Measured viscosities of a large number of coal ashes show that, 
as a slag is cooled from the true liquid state, its viscosity increases 
along a logarithmic curve until some point is reached where it 

2 Theoretical furnace temperature as used here is the calculated 
temperature obtainable neglecting dissociation but allowing for the 
heat loss to the moisture and to the slag. 
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Fig. 2 Relation of theoretical furnace temperature to iemperature of slag 
required for a viscosity of 250 poises 
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Fig. 3 Viscosity-temperature plots of typical slags (refer to Table 1 for 
analyses as listed by number) 


deviates and viscosity tends to increase more rapidly as tempera- 
ture is decreased. The point at which this deviation takes place 
has been termed ‘‘the temperature of critical viscosity’’ (7..), 
and its occurrence will be discussed later. The deviation may be 
gradual or it may take place rather sharply (Fig. 3), and is 
described as the point at which solid material begins separating 
from the liquid, resulting in a mixture of both.. Measurements 
beyond this point are no longer true viscosities; however, they 
do give an apparent viscosity which can be used as a measure of 
flow characteristics. Through proper correlation of these data, 
the measured viscosity values of a slag can be related to its 
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chemical composition and to the degree of oxidation of the iron 
in the slag. 

Viscosity in the Liquid Range. There are many ways by which the 
relationship of viscosity in the liquid range to chemical analysis 
may be correlated. However, the method proposed by the Bureau 
of Mines is comparatively easy to use and gives a good correla- 
tion. This method correlates ash viscosity in the true liquid range 
to the factor 

SiO. & 100 
SiO, + Fe.O;* + CaO + MgO 


In this paper this factor is designated as the ‘‘equivalent-silica 
percentage.’’ In the true liquid range all slags having similar 
equivalent-silica percentages will have a similar viscosity-tem- 
perature relationship. Based on experimental data a family of 
logarithmic curves can be established showing the viscosity- 
temperature relationship for a range of equivalent-silica contents. 
Such a series of curves is shown in Fig. 4. 

While this method may appear to neglect the effect of alumina, 
this is not exactly true. Measured viscosities show that as 
alumina replaces silica in a slag analysis, the viscosity decreases. 
Also, as alumina replaces silica the calculated equivalent-silica 
percentage lowers, which also means a lower calculated viscosity. 
Hence this method does in reality compensate for changes in 
alumina content. This method holds reasonably well for silica- 
to-alumina ratios from 1 to 4. 

Point of Critical Viscosity. Measured viscosities of many slags 
show that, as the slags are cooled, viscosity increases uniformly 
along a line approaching that calculated as described. It will be re- 
called, however, that, as cooling proceeds, a point may be reached 
where viscosity increases more rapidly as a drop in temperature 
This is caused by a selective separation of solid 
Many authorities refer to this as the 
Limited attempts to identify material 


takes place. 
material from the liquid. 
start of crystal formation. 
in the viscometer melt have generally been unsuccessful, since 
most of the material is noncrystalline. However, Fe.O; and 4 
Ca0-Al,Os-Fe:O; have been identified by x-ray diffraction in a 
melt from a coal-ash high in calcium and iron. Nevertheless, the 
separation of the solid material is a real phenomenon and in many 
causes can actually be seen 

The point at which this separation of solid material begins is 
not always sharply defined and may occur at higher tempera- 
tures by extremely slow cooling rates. For the work covered in 
this paper, all slags were cooled in steps of 30 deg over a period of 
30 min except in zones where rapid changes in viscosity occurred. 
\t each step viscosity was observed for 5 min to determine if a 
change was taking place. This cooling rate is considered slower 
than the chilling of slags during their passage from the combus- 
tion zone to the slag tank. 

The point at which the temperature of critical viscosity occurs 
is dependent largely upon the iron content of the slag and on the 
degree of oxidation of the iron. The degree of iron oxidation is 
normally expressed as the ferric percentage where 


Feel )s 
xX 100 


Ferric percentage 
. 1.11 FeO + 1.43 Fe 


Fe.O; + 
In stating the viscosity of a slag containing iron it is of prime im- 
portance that the degree of oxidation, ferric percentage, be ex- 
Fig. 5 shows viscosity curves for a typical slag run 
under various stages of oxidation. Note that, in the tapping 
range of 50 to 100 poises, ferric percentage has a very marked 
effect on viscosity. In making viscosity determinations, there- 
fore, it is not sufficient to know the type of atmosphere in which 
the determinations were made, but samples must be periodically 
taken to determine ferric percentage, and the atmosphere must 
be adjusted to obtain the desired ferric percentage. 


pressed. 


* Total iron as FesOs. 
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Fig. 5 Viscosity-temperature plots of a typical slag showing effect of 
ferric percentage 


Samples of slags from a large number of slag tap furnaces 
operating under normal conditions at 10 to 15 per cent excess air 
show an average ferric percentage of 20. Bureau of Mines data 
[4] also show approximately a 20 per cent ferric average with a 
range of 6 to 34 per cent on 16 pulverized fuel units that they in- 
vestigated. Hence a means of predicting the temperature of 
critical viscosity and the apparent viscosity below this temperature 
at 20 per cent ferric is desirable. We find that the ash softening 
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temperature in a reducing atmosphere (hemispherical point*) 
plus 200 F gives a fair approximation of this point for the purpose 
for which it is used. Fig. 6 shows a plot of actual determined 
temperatures versus ash softening temperature in a reducing at- 
mosphere plus 200 F trom data on Table 1. There is reasonable 
agreement with the 7... Some of the discrepancy results from 
the determination of 7',, at other than 20 per cent ferric. After 
the 7., has been established it can be located on the liquid vis- 
cosity line on the plot as shown in Fig. 4. If this point falls below 
250 poises, extending a line upward at a 10-deg slope gives a 
reasonable approximation of the temperature-viscosity relation- 
ship ‘ip to 250 poises. An example of how this is applied is shown 
in Fig. 4. Fig. 7 shows an approximate method of correcting 
temperature for other than 20 ferric percentage to 20 per cent at 
50, 100, and 250 poises viscosity. 

A comparison of viscosities as obtained by this method with 
actual measured viscosity corrected to 20 ferric percentage is 
shown in Fig. 8. With very few exceptions the predicted viscosi- 
ties fall within a band of +100 F from measured values. 


Base-to-Acid-Ratio Method 


The constituents of a coal ash can be classed as either basic or 
The acidic constituents are the silica, alumina, and 
The basic constituents are the iron, calcium, magnesium, 
For coals mined in the U. 8., the acid materials al- 
The relative amounts of the 


acidic. 
titania. 
and alkalies. 
most always exceed the basic. 
basic and acidic constituents in the ash can be used as a means of 
predicting the viscosity of the slag and also as a means of deter- 
mining what substances must be added to the ash to obtain de- 
sired flow characteristics. These values can be expressed as the 


base-to-acid ratio, which is defined as 


FeO; + CaO + MgO + NasO + KO 
SiO, + Al.O; + TiO, 


The viscosity of a slag decreases as the base-to-acid ratio in- 
creases to 1. Although our investigation beyond this point is 
limited, other investigators indicate that viscosity is fairly uni- 
form at base-to-acid ratios between 1 and 8 [6]. 

The slag temperature for any given viscosity can be plotted 
Plots for 50, 100, and 250 poises 
In slags having a low base-to- 


against the base-to-acid ratio. 
are shown on Figs. 9, 10, and 11. 
acid ratio, the temperature required for a given viscosity in- 
creases more rapidly as the silica-to-alumina ratio becomes higher. 
Ranges for a silica-to-alumina ratio below 1 have not been in- 
vestigated; however, this trend is expected to reverse as the 
alumina content exceeds that of the silica. : 

A comparison of viscosities as obtained from the curves on 
Figs. 9, 10, and 11 with the actual measured viscosity when cor- 
rected to 20 ferric percentage is shown in Fig. 12. The viscosities 
predicted by the base-to-acid ratio show about the same accu- 
racy as the modified Bureau of Mines method. 

Only a limited amount of work has been conducted on slags 
where the potassium or magnesium content was of sufficient mag- 
nitude to permit evaluating their effect in relation to other basic 
constituents. Care should be used in extending these curves to 
slags containing major proportions of these two materials. 
Other authorities [7, 8] have suggested using multiplying factors 
for the individual elements. However, within the range found in 
U.S. coal-ash slags we find this does not appreciably improve the 
correlation. 


Fluxing 


Under certain conditions a slag may be altered to obtain de- 


3 The temperature at which the cone has fused down to a hemi- 
spherical lump in which the height is one-half the width of the base. 
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Fig. 7 Amount of change in temperature at a stated viscosity at other 
than 20 per cent ferric 


sired flow characteristics. Fluxing or blending of slags to lower 
viscosity is generally effective; however, an increase in viscosity 
by the addition of materials is not so readily attained. In order 
to be completely effective all materials in the base slag and the 
added material must be intimately mixed and melted. 

Figs. 9, 10, and 11 offer a means of predicting the effect of al- 
tering a slag by the addition of a flux or asa result of blending with 
other slags. The weights of the respective materials in the re- 
sulting mix are calculated and the new base-to-acid ratio ob- 
tained. Values of materials added should be expressed as 
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Fig. 8 Plot of temperatures for measured viscosity corrected to 20 per 
cent ferric versus viscosity as calculated by modified Bureau of Mines 
method at 50, 100, and 250 poises 

oxides. The resulting temperature for 50, 100, and 250 poises at 
20 ferric percentages can be read directly from the charts. 

The selection of a flux is usually dictated by economic reasons. 
In many cases the materials to be used in this way may be availa- 
ble locally in nature or as waste materials. The following is 
offered as a guide in selecting a flux. 

Calcium Compounds. The calcium compounds are the most com- 











Slag Temperature for 50 Poises Viscosity i 


0.6 
Rase-to-icid Ratio 


mon type of flux and are usually obtainable in the form of lime- 
stone, CaCO;, or as the hydroxide, Ca(OH)». 

Iron Compounds. Mill scale or other high-iron compounds may 
effectively be used as a flux; however, they usually command a 
premium price. Iron or iron compounds tend to separate from 
the slag under adverse operating conditions, hence some care 
should be exercised in their use. 

Magnesium Compounds. Magnesium compounds are very simi- 
lar to calcium in fluxing ability. Manufactured magnesium com- 
pounds are usually costly. However, many localities have 
natural deposits of dolomite or dolomitic lime, which are suitable 
fluxes. 


Sodi Cc d 
Pr 





Sodium compounds are very effective fluxes 
and are often used in the form of soda ash, Na2CQ;; or salt cake, 
Na2SO,; for emergencies. Published literature indicates that 
sodium compounds probably contribute to boiler fouling, hence 
their continuous use is not desirable. 

Potassium Compounds. The effect of potassium compounds on 
slag viscosity has not been investigated to any extent. They are 
not considered as effective as the other compounds listed. Some 
authorities [8] suggest that in certain cases they may increase 
viscosity. 

Ash Blending. A high-viscosity ash and a low-viscosity ash can 
often be blended to obtain more desirable flow characteristics. 
One practical application of this is to fire wood and bark in com- 
bination with coal. Most wood and bark ashes have a high cal- 
cium content, which will appreciably lower the viscosity of the 
coal slag. 








Conclusions 


The goal of this paper has been to show how laboratory 
analysis of a coal and its ash can be practically and reliably used 
to predict the slag behavior in the furnace of a power boiler. This 
information is of pertinent value to the boiler manufacturer in the 
design of his equipment, to the boiler operator in comparing and 
selecting the most suitable fuel sources, and to the coal producer 
in determining the suitability of his fuel for use in the specified 
equipment. 
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Fig. 9 Plot of temperature for 50 poises viscosity at 20 per cent ferric versus base-to-acid ratio 
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Fig. 10 Plot of temperature for 100 poises viscosity at 20 per cent ferric versus base-to-acid ratio 
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A method has been shown whereby: 


1 The behavior of a slag in the tapping range can be defined 
by its temperature-viscosity relationship. Pree 

2 Viscosity limitations on a known boiler design can be estab- 
lished. 

3 Viscosity-temperature relationship of an ash can be de- 
lineated in the tapping range from the chemical analysis and 
fusion temperatures. 

4 An unsatisfactory ash may be altered and the flux require- 
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Plot of temperature for 250 poises viscosity at 20 per cent ferric versus base-to-acid ratio 


ments predicted to obtain the desired viscosity-temperature rela- 
tionship. 


Through these means, therefore, it is possible not only to match 
flow characteristics with ash analysis but also to modify these 
characteristics to suit the requirements of a given design, when- 
ever this proves to be economically feasible. 
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APPENDIX 
The viscometer used for this work is of the rotational type as 
illustrated in Fig. 13. The liquid under study is contained in a 
platinum-rhodium crucible, cylindrical in shape, and a cylindrical 
bob is rotated in the liquid at a constant speed through a cali- 
brated suspension wire. The torque or amount of twist produced 
The 


amount of twist is measured and recorded as the interval between 


in the suspension wire is proportional to the viscosity. 


impulses from light beams reflected from mirrors attached to the 
ends of the wire. The strip chart recording of the twist, as shown 
in Fig. 13, provides a convenient record and method of measure- 
ment. The suspension wires used are calibrated against viscosity- 
standard oils obtained from the Bureau of Standards and a 
glass measured and supplied by the Corning Glass Company. 
The electrically heated furnace is of the Globar-tube type with 
temperature regulation provided through a controlling type po- 
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tentiometer actuated by a thermocouple located in the furnace 
adjacent to the sample crucible. A thermocouple imbedded in 
the slag-crucible support serves to indicate sample temperature. 

The ash sample, 150-200 gram minimum, is prepared by 
standard laboratory methods and pelletized. Sufficient pellets 
are fed into the platinum crucible to obtain the desired slag level. 

Atmosphere Control. Since the oxidation level of the iron in the 
slag has a marked effect on viscosity it is necessary to control the 


















































CULATED TEMP, POR 50 POISES °F 
CALCULATED TEMP. FOR 100 POISES °F 











~~ 2800~«*«S 0] 


MEASURED TEMP, FOR 50 POISES °F MEASURED TEMP, FOR 100 POISES °F 





—y 


+ 


i 


N 








oH Ne 


| | 


iS 








CALCULATED TEMP, FOR 250 POISES °F 











2000K_ 
2000 2800 
MEASURED TEMP, FOR 250 POLSES °P 





Fig. 12 Plot showing correlation between measured viscosity and 
calculated viscosity at 20 ferric percentage for base-to-acid ratio method 
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ferric percentage in the slag. This is accomplished by holding a 
hydrogen-nitrogen atmosphere over the surface of the melt. The 
ratio is varied to give the desired ferric percentage. A mixed gas 
flow of 0.5 cfh containing 5 per cent hydrogen, after passing over 
heated copper to remove any oxygen and then being dried, is re- 
quired to hold about a 20 ferric-percentage level. 

Ferric percentage is determined periodically from samples 
taken from the melt by inserting a heavy platinum wire. The 
wire is withdrawn and plunged into water to prevent any oxida- 
tion and to break off the slag. 

The ferrous-iron and total-iron contents of the slag sample are 
determined by standard analytical methods and the ferric per- 
centage is calculated. 

Test Procedure. After the sample is heated to well above true 
liquid temperature and the desired ferric-percentage level ob- 
tained, the slag is cooled in controlled steps of 30 F over a 30-min 
period and viscosity determination made at each step. When 
changes in viscosity are very rapid, cooling may be in smaller 
steps. Several calibrated suspensions are required to cover the 
range of 0 to 10,000-poises viscosity. 

The normal viscosity determination including analysis requires 
about 80 man hours. 


DISCUSSION 
William T. Reid* 


Twenty years ago the Bureau of Mines provided the first real 
insight into the complicated relationship between slag viscosity 
and composition. Since then, these data have furnished the only 
practical means of predicting the flow characteristics of coal-ash 
slags from their chemical analysis. Such data have been helpful 
in the few cases where fluxes have been used to make tapping 
easier, they have been of assistance in some cases in selecting 
coals for mixing to get slags of desired flow characteristics, and 
they have provided basic understanding of the factors influencing 
the buildup of slag on heat-receiving surfaces. In general, how- 
ever, during these twenty years furnace designers learned how to 
improve the tapping ability of slag-tap furnaces, and_ boiler 
operators found by experience how to select coals whose ash 
produced easily tupped slag. 

More recently, new problems with slags in the cyclone furnace 
led the authors to reopen this study. Thanks to them, the 
moderately wide range of compositions studied originally by the 
Bureau of Mines has now been expanded to include coal ashes of 
widely varying silica:alumina ratio, and with FeO; contents 
ranging from 2.4 per cent to 41.2 per cent and with CaO + MgO 
from 1.2 per cent to 30.8 per cent. It is gratifying indeed to note 
that the general relationships determined in the Bureau of Mines 
studies have been confirmed by this work. Such confirmation 
seems to have been provided not only in the viscosity-temperature 
data but in the concept of the temperature of critical viscosity 
and the effect of ferric percentage as well. 

It is equally interesting to note that recent work in Great 
Britain, as reported by Shaw,> also confirms these results. 
Hence it appears now that the relationship between composition 
and viscosity above the temperature of critical viscosity is in good 
shape. A few unusual coals may still turn up which appear to be 
outside the predictable range. But, in general, the high-tem- 
perature viscometer can be considered to have served its purpose. 
Routine measurements of viscosity are no longer necessary when 
the composition of the ash is known. 

A problem area still exists below the temperature of critical 
The 


viscosity where the slag behaves as a thixotropic material. 


4 Assistant Technical Director, Battelle Memorial Institute, 
Columbus, Ohio. Mem. ASME. 
5 J. T. Shaw, ‘‘Coal-Ash Slags. Part II Flow Properties,’’ Mon. 
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development of yield stress in such slags certainly has considera- 
ble importance in slag deposits on wall tubes. Yet the rate at 
which such internal structure develops in coal-ash slags and the 
factors that affect it are still largely unsolved questions. This is 
a difficult problem area, but it well appears to be of enough im- 
portance to justify extensive work. 

Another problem is the rate and the mechanism by which slags 
form from coal ash. Here the mineralogical composition of the 
coal and ash may well be an important variable, and certainly 
the temperature-time history during combustion and subsequent 
heating will greatly affect the slagging reactions. Only a start 
has been made in studying this difficult and complex problem, yet 
the same status of knowledge on it as on slag viscosity would be of 
great value. 

Hence it now appears that the care and attention paid by the 
authors and others to the problem of slag viscosity has reliably 
solved the viscosity-composition problem for most coal-ash 
slags. With this problem out of the way, perhaps it is time to 
make an equally determined attack on the slagging mechanism 
and the behavior of slags below the temperature of critical vis- 
cosity. These are problems of sufficient magnitude to test any re- 
searcher’s skill. 


D. M. Wilkins® 


Investigations have been carried out in the BCURA for 
many years with combustion systems, and more recently with 
gasification systems, from which the coal ash is removed mainly 
as liquid slag. Many data have been accumulated on the chemi- 
cal composition of the ash of British coals of the type widely used 
for power generation.?. The composition of the slags formed by 
many of these coals when burned in an experimental cyclone 
combustion chamber and in industrial boiler installations has also 
been determined. During the last three years a viscometer has 
been in operation in connection with a program of research on 
gasification sponsored by the Ministry of Power, and a large 
number of viscosity determinations have been made on both 
fluxed and unfluxed coal ash slags.. We thus fully appreciate the 
effort which has been entailed in obtaining the results published 
in this paper. The authors are to be thanked for revealing a 
synopsis of the actual results obtained rather than, as is so often 
the case, presenting only an interpretation of the findings. 

The correlation between silica ratio and slag viscosity in the 
liquid range which was proposed by Reid and Cohen’ has in gen- 
eral been substantiated by the results of the viscosity determina- 
tions which we have carried out with coal ashes and coal ash 
slags. When, however, the slag or ash contained a high pro- 
portion of CaO or CaO + MgO, as occurs with fluxing, the dif- 
ference between the determined viscosity and that predicted 
from the silica ratio was found to be large. <A correlation between 
equiv. Fe.O; + CaO + MgO 

SiO, + ALO, 
in quite good agreement. This ratio is similar to that proposed 
by Sage and McIlroy except that we have found no worthwhile 
reduction in the scatter of results by including terms for the al- 
kali, titania, or Mn;O, contents. Incidentally, the base-acid ratio 
used by the authors seems to be unnecessarily complex. As is 
shown in Table 1 of the paper, the eight oxides represented in the 
ratio together make up all of the slag, so that the ratio might be 
; 100 — (SiO, + ALO; + TiO.) 
written equally well as ——— a 
(SiO. + ALO; + TiO.) 

6 The British Coal Utilisation Research Association, Leatherhead, 
Surrey, England. 

7H. R. Hoy, A. G. Roberts, D. M. Wilkins, Journal Institute of 
Fuel, vol. 31, 1958, p. 429. 

8 W. T. Reid and P. Cohen, ‘‘The Flow Characteristics of Coal-Ash 
Slags in the Solidification Range,’’ Trans. ASME, Furnace Per- 
formance Factors, vol. 66, 1944, pp. 83-97. 


viscosity and in such cases results 
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Table 2 


Range of 
ferric, 
per cent 
circa 20 


Source of 
information 
Fig. 6 
Table 1 
Table 1 
Table 1 


= ().927 AST 


TiO, content of the samples is small and reasonably constant, the 
eee W — (SiO. + Al.O;) 
expression may be further simplified to = 
(SiO. + AlsO;) 

with great saving of analytical effort. 
Of the coal ashes examined by us 
approximately 1 in 10 have had equivalent FeO; 


mainly from untreated 
small coal 
contents greater than 15 per cent and very few of these have ex- 
ceeded 20 per cent FesO;. The graphs relating temperature of 
critical viscosity and chemical composition which were published 
by Reid and Cohen* in their paper to the ASME in 1944 are there- 
fore not applicable to the majority of our ashes and slags without 
extrapolation. At present we determine the temperature of 
critical viscosity experimentally in the viscometer and supple- 
ment the information by the use of crystallization techniques. 
However, the ability to predict Tey from chemical composition 
and/or a simple physical test would be highly advantageous. 
My main interest in the author’s paper is, therefore, centered on 
this possibility. 

It is noted that about two-thirds of the authors’ results refer 
to ferrie percentages within the limits of 6 to 34 per cent which 
were stated to apply to slags tapped from industrial boilers in the 
U.S. A. 

Since it is normal practice to operate modern slagging boilers 


Do the remainder of the results refer to coal ashes? 


with the minimum of excess air (1.5 — 2.5 per cent O2), the upper 
limit of 34 per cent for the ferric percentage appears to be high. 
In our experience the ferric percentage of slags from combustion 
systems seldom exceeds 10 per cent and is occasionally zero. Do 
these slags with high ferric percentage originate from plant work- 
ing under realistic operating conditions? 

The authors claim that there is a usable relationship between 
temperature of critical viscosity and ash softening temperature 
in a semireducing atmosphere and that some of the discrepancy 
from the relation Tey AST + 200 F results from determination 
of 7. 
of ferric percentage in this connection we have derived the follow- 


it other than 20 per cent ferric. To assess the importance 
ing equations from statistical analysis of the information pre- 
sented in Table 2° and Fig. 6 of the paper. 

In Table 2, none of the confidence limits differ significantly 
from each other. 


of ferric percentage is increased, scatter from the mean line does 


Furthermore, it will be noted that, as the range 
not increase significantly. Although these expressions afford a 
better estimate of T.y from knowledge of AST than those pro- 
posed by the authors, they are still insufficiently precise for prac- 
tical purposes and consequently the value of determining the ash 
Would the authors indi- 
cate what advantage the relation T'cy AST + 200 F has for 
American coal ashes over the relationship between chemical 


softening temperature is not apparent. 


composition and 7'cy given by Reid and Cohen.$ 

In the example given in Fig. 5 of the paper, a marked change in 
Tey With differing levels of oxidation is shown at ferric percentage 
less than 49 to 59 for a slag having a high equivalent Fe.O; con- 
tent. 
the solid which separates on cooling does not contain iron the 
with differing levels of oxidation may be very 


For slags with lower equivalent Fe,O; content and where 


change in 7, 

* There is apparently disagreement between the temperatures of 
critical viscosity given in Table 1 and those plotted in Fig. 3, code 
Numbers 25, 15, 34, 49, and Fig. 5, code Number 11. Consequently, 
some incorrect data may have been used in calculating the correla- 


tions. 
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Regression equation 
‘ev = 0.750 AST + 840 F 
+ 415 F 
, = 0.842 AST + 608 F 
‘wv = 0.630 AST + 1177 F 


AST 

range, 

deg F 
2180-2770 
2030-2550 
1980-2660 
1980-2780 


95 per cent Degrees 
confidence of 

limits freedom 
+227 F 32 
+165 F 18 
+162 F 32 
+184 F 47 


much smaller. Nicholls and Reid” found that for equivalent 
Fe,0; contents less than 12 per cent, little change in Ty resulted 
from varying the degree of oxidation. Our own experimental re- 
sults appear to confirm this. While Fig. 5 could serve as a warning 
to the imprudent plant operator it would have added to the value 
of the paper if the authors had shown the effect of varying the 
ferric percentage over the lower range applicable to modern slag 
tap boiler practice. 

Referring again to Fig. 5 in the paper, I should be interested 
to know whether the range of ferric percentage given for each 
curve is due to temperature changes or to changes in atmosphere. 

The temperature of critical viscosity is normally taken to be* 
the temperature at which on cooling a liquid slag the concentra- 
tion of particles increases to such an extent that there is a sudden 
transition from liquid to plastic flow. In the plastic range, vis- 
cosities will be dependant upon rate of shear and flow will not 
take place until the internal yield stress is exceeded. If these 
concepts are accepted, the apparent viscosities obtained by the 
authors at temperatures below Tey are of questionable value in 
practice unless the relationships between apparent viscosity and 
both rate of shear and internal yield stress are given. For this 
reason the practical application of the curves given in Fig. 7 of 
the paper must surely be very limited. 

I think that the authors will agree that the alternative interpre- 
tations of some of the experimental data given in their paper ma) 
be feasible. However, what really matters is whether the data 
presented by the authors allow prediction of the behavior of slag 
in an industrial plant. I would therefore ask the authors whether 
they have collected sufficient information to enable them to 
assess the relative merits for using: 

(a) silica ratio alone, (b) silica ratio plus the 7'.v-ash softening 
temperature relationship proposed, and (c) the determined tem- 
perature viscosity characteristics of the slag, to predict the 
suitability and operating limits for a coal in a slag tap boiler. 


Authors’ Closure 

The authors wish to thank Messrs. Reid and Wilkins for their 
enlightening comments on this paper. 

Mr. Reid points out the need for additional research study on 
slag behavior in the temperature range below the temperature of 
In this range physical and chemical changes 
The authors’ Company 


critical viscosity. 
occur making the study quite complex. 
is continuing work in this field in the hope that it may shed some 
light on the problem of wall slagging. 

Comments by Mr. Wilkins indicate a rather careful scrutiny of 
the data and point out an error! missed by the authors during the 
final editing stage. 

In general we concur with Mr. Wilkins’ comments. It is true 
that, in the vast majority of coals currently being used, the 
titania content is sufficiently low so that its omission would not 
appreciably affect the results. In a majority of slags, the alkali 
contents are low and their omission or the use of an average value 
does not appreciably affect the results. However, NazO does have 

1 P, Nicholls and W. T. Reid, ‘‘Viscosity of Coal-Ash Slags,” 
Trans. ASME, vol. 62, 1940, pp. 141-153. 

11QOn Fig. 3 curves numbered 25, 34, and 49 should have been 
numbered lle, 33, and 48 to correspond to the analyses given in 
Table 1. 
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very good fluxing properties and its omission on high-alkali slags 
or in cases where sodium salts are used as a flux can lead to con- 
siderable error. On the other hand, we can see no advantage to 
expressing the base-to-acid ratio in the manner suggested by Mr. 
Wilkins, because if the sum of the SiO., AlsO;, TiO», Fe.O:, CaO, 
and MgO and alkalis is made equal to 100 per cent the manner 
chosen to express the ratio is merely a matter of mathematical 
expediency on the part of the user. 

The 7v (temperature of critical viscosity) in a slag is thought 
to result from the presence of iron in the slag. Viscosity studies 
by J. L. Machin, et al. [7] on mixtures containing the major in- 
gredients of iron-free coal ash do not show this phenomenon. 
Our studies of slags having a very low iron content have also 
failed to show a temperature of critical viscosity or, at most, one 
that is only faintly defined. 

We agree with Mr. Wilkins that the use of AST + 200 F is not 
a very precise means of determining the temperature of critical 
viscosity; however, it is much simpler to determine it in this way, 
and its accuracy is comparable to the other methods suggested. 

The data shown in Table 1 are largely from determinations 
made on samples obtained by laboratory ashing of a coal; how- 
ever, they do contain some ash or slag samples obtained from 
operating units. In a few cases, laboratory mixes of chemicals to 
give the same composition were checked against the ash samples. 
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It is the authors’ belief and the basis of this paper that a given 
chemical mixture, after being heated to a true liquid condition 
and controlled at the same level of oxidation, will behave in ex- 
actly the same manner on cooling, irrespective of the source of 
these chemicals. 

The ferric percentages shown for the slags in Table 1 are en- 
tirely dependent on the atmosphere within the viscometer and 
have no bearing on boiler operation. In this paper, the authors 
suggest that, when calculating the viscosity of an ash from its 
chemical composition, the ferric percentage be assumed to be 
20. This percentage will give results on the safe side. We feel 
that Mr. Wilkins’ figure of 0-10 per cent is on the low side for 
typical boiler practice over the range of operating conditions en- 
countered. Actually, calculated viscosities in the tapping range 
do not differ significantly whether a ferric percentage of 10 or 20 
is used. 

The procedures outlined in this paper are used by the authors’ 
Company for predicting the suitability of coals for use in liquid 
ash removal furnaces. Representative coal samples are ashed by 
standard procedure and ash constituents are determined by spec- 
trographic means. 

The silica ratio and Tv are used only in the modified Bureau of 
Mines method of calculation shown. If the calculated Tey occurs 
at very low viscosities (under 50 poises) we prefer to use the base- 
to-acid ratio for predicting viscosity. 
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Pressurization of Granular Solid Fuels 


M. N. AREF 


Project Supervisor, 
Research Division, 
Foster Wheeler Corporation, 
New York, N.Y. Mem. ASME 


The different methods of pressurizing powdered solid fuels which have been developed 
are reviewed briefly, including proposals as well as those in operation. The work has 
been devoted particularly to a study of the characteristics of these pressurizers, in rela- 
tion to their common problems of: (a) Power loss due to venting the pressurizing air, 


(b) wear of the surfaces in contact with the solids, and (c) effective sealing. The most 
practical methods have been compared on power loss, rather than power required in 
their cycles. It has been found that power loss due to air venting cannot be totally 
eliminated and that the simple lock hopper is, so far, the best available device which has 
the lower capital cost, maintenance requirements, power loss, and better reliability. 
Where space is of major importance, however, the lock hopper with a moving wall can be 
considered as an alternative, though it may suffer from an increased wear rate and a 
reduced volumetric efficiency. 


0. OF THE MAIN PROBLEMS which faces the designer 
of a solid-fuel system for a pressure-fired device is how to intro- 
duce the pulverized or granulated particles into the combustion 
space against the system pressure. By far, the most convenient 
method of achieving this is to employ a solid-fuel pump de- 
signed to do much the same duty as a liquid pump. The prime 
need is to pressurize the fuel particles, or more accurately to pres- 
surize the air, if it is a fuel-air mixture, or the water in case of a 
slurry, in which they are entrained. 

Until recent years, no such pressurizer had been required and 
very little work had been done on the pumping of solids, but the 
idvent of the newly developed coal-burning processes in combus- 
tion and gasification at elevated pressures provided a stimulus to 
research in that direction. Among these processes were the 
coal-fired gas turbine, pressure-fired boiler for the steam-gas 


cycle, and pressure gasification. In these processes, the solid-fuel 


may be pulverized, crushed, or sized in the case of coal; or sized 
in the case of other materials. Loose-packed bulk densities may 
range from 20 lb/ft® for sized coke, through 30 for pulverized 


Contributed by the Fuels Division and presented at the Annual 
Meeting, Atlantic City, N. J., November 29-December 4, 1959, of 
THE AMERICAN Society OF MECHANICAL ENGINEERS. 
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1959. Paper No. 59—A-150. 


coal, 40 for sized coal, 55 for crushed coal, 100 for stone, up to 200 
for rich ores. Capacity may vary from 2 tons/hr for a small gas- 
turbine plant to 100 tons/hr or even higher for a big gasification 
plant. In the latter case, however, multiple units may be used 
to reduce spares, so throughputs of 2 to 25 tons/hr cover the 
likely unit range. Delivery pressures may range from 5 to 20 atm 
and the powdered fuel size, for satisfactory pumping, may vary 
from 80 per cent through 200 mesh for pulverized down to 0-'/, in. 
or 0-'/, in. for crushed. 

In considering the design of a solid-fuel pump for a unit ca- 
pacity of, say, 25 tons per hr of crushed or pulverized coal at 20 
atm as the primary objective, the following basic methods deserve 
further study: 

1 Cyclic action group: 

Simple lock hopper 
Lock hopper with a moving wall 
Positive-displacement ram pump 
Solid extrusion 
Slurry pump 

Continuous-action group: 
Positive-displacement gear pump 
Multistage centrifugal compression with air 
Peristaltic-displacement (rubber) pump 


So far, established methods employ the lock-hopper principle, 





Nomenclature 
surface area, ft? 
constant 


diameter of discharging orifice of 
hopper, D, cylinder diameter, 


pressure differential across the 
top and discharging sides of the 
solid-fuel in a hopper 

mass flow rate of solid-fuel, 


tons/hr, Qasp positive aspirated 


power loss in the cycle in hp-hr/ 
ton of coal 


fraction of water content in the 
fuel slurry 


tangent of the angle of repose of 


ft 
mean particle diameter 
function 
acceleration due to gravity 


specific weight per ft* of solid or 
gas, VM. of coal, M, of air 


number of hoppers if more than 


one hopper is used 


pressure, P, atmospheric, P,, 
P,.. P, above atmospheric 
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flow, Qatm atmospheric gravity 
flow 

gas constant of air taken 53.3 ft 
lb/deg F-lb 

compression ratio 

weight ratio of coal to air 

weight ratio of coal to water in a 
slurry 

absolute temperature, deg R; 7, 
atmospheric, 7, due to com- 
pression 

specific volume of air at pressure 
P in ft3/ton of coal 


the granular solid 

adiabatic compression index 
equals 1.4 

efficiency, Nisotn due to isothermal 
compression 

dynamic coefficient of dry fric- 
tion 


density, p, absolute solid density, 
p, bulk density 


cone angle of hopper feeding 
orifice 
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two extreme forms being the simple lock hopper which may be 
bulky and slow but reliable and the small fast lock hopper with a 
moving wall which has a problem of achieving an efficient rubbing 
seal. 

The main features of most of these methods have been given 
elsewhere, but the purpose of this work is to study their charac- 
teristics more closely, particularly with regard to the common 
problems of a solid pressurizer, which are: (a) Power loss due to 
venting of the pressurizing air, (6) wear of the surfaces in contact 
with the solid, (c) effective and reliable sealing. 

In this work, the first two problems have been studied. 
ever, as a basis for planning development, such study may offer 
economical solutions for high-capacity, high-pressure duties. For 
convenience, solid fuel has sometimes been referred to as coal. 


How- 


Factors Governing the Rate of Emptying 

It is prerequisite, for proper design of a granular-solid pressur- 
izer, to be able to assess with sufficient accuracy the flow of a 
granular mass through an orifice. Deviating considerably from 
the principles of fluid flow, granular solid flow is quite different 
from that of a fluid even in its very viscous condition due to the 


following reasons: 


1 With liquids, flow is proportional to head, and pressure is 
equally transmitted in all directions, whereas with granular solids, 
the flow is almost independent of head, and the pressure is not 
uniform throughout the mass. 

2 The formula for fluid flow cannot be applied to solids be- 
cause the term of the particle size is absent. 

3 Angle of repose applies only to solids. 

t Bridging of solids due to cohesive forces between a particle 
and neighboring particles or wall. 

In the absence of an exact theory, several investigators have de- 
veloped empirical correlations between dimensionless groups for 
gravity flow from inverted cones or from a centrally situated 
circular orifice at the base of a tube. Assuming that the flow Q 
depends on the relevant quantities: D, d, p,, p,, 8, ®, and g, then 
according to Buckingham’s ‘“‘pi’’ theorem, five different dimen- 
sionless groups may be derived, yielding the following expression 


Q i 
p,D*5g9-5 > 


Deming and Mehring [1]! managed to evaluate experimentally 
the function on the right side of equation (1) and gave one ex- 
pression in d/D, 8, and ®. In order to plot equation (1), prefer- 
erence has been given to a certain type of bituminous coal of an 
absolute density of 78 lb/ft’, to find the variation of its bulk 
density with particle size. This has been determined experi- 
mentally, as shown in Fig. 1, for several samples of sized particles 


Fld/D, p/P. 8B, ®) (1) 


1 Numbers in brackets designate References at end of paper. 
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of this coal. However, since the bulk density is not a reproducible 
quantity, two extreme cases were measured; namely for loose 
gravity packing and compact shaken packing. It was found that 
the latter bulk density exceeded the first by approximately 8 per 
cent. The angle of repose has also been measured for different 
particle size and is shown in Fig. 2. By substitution from Figs. 1 
and 2 in equation (1), Fig. 3 has been produced for d/D varying 
from 0.05 to 0.2 only, since, for values beyond this range, equation 
(1) cannot be directly adopted without appropriate correction. 
This explains that other factors limit the validity of equation (1), 
such as excessive cohesion due to electrostatic forces in the case 
of very small particle sizes and bridging in the case of large sizes. 

Means of Speeding Up Solid Flow. In order to accelerate the 
gravity flow of the solid fuel from a hopper, positive aspiration 
can be applied on the solid fuel by establishing a positive pressure 
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differential across the top and discharging sides of the fuel. It is 
possible also to augment this pressure differential by lengthening 
the orifice constriction at the discharge below the hopper, due to 
solid-fuel plug in the constriction acting like a piston. This in- 
creases the gravity force on the plug and consequently results in 
a better flow of the granular solid. 

The theoretical relationship between the flow rate under posi- 
tive aspiration and that obtained by gravity has been derived and 
is given in Appendix 1. This has the following expression: 


Qusp = Qaim V1 + AP/Poem (2) 


Fig. 4 is a plot of equation (2) and illustrates the manner Qasp is 
increased over Qaim by applying positive aspiration, assuming 
that the bulk density remains the same and the mean particle 
diameter is constant in both flow cases. 
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Cyclic Action Group 


Simple Lock Hopper. 
izing granular solids is the lock-hopper system. This employs, as 


Probably the simplest method of pressur- 


shown diagrammatically in Fig. 5, two high-pressure closed hop- 
pers connected in series to an atmospheric hopper. The sequence 
of their cycle is as follows: 
1 Valve 1 is opened. Solid is introduced by gravity at atmos- 
pheric pressure into the lock hopper, valve 2 being closed. 
2 Valve 1 Auxiliary pressurizing air is applied 
through valve 3 to the lock hopper until it reaches the same pres- 


1s closed 


sure as that of the feeder hopper 
Solid is then displaced by gravity into 
Valve 4 serves 


3 Valve 2 is opened. 
the pressurized feeder hopper at system pressure. 
to equalize pressure in both hoppers. 

4 Valves 2, 3, 


pressurizing air to the atmosphere and cycle is repeated. 


and 1 are closed. Valve 5 is opened to vent 


Apparently, in each cycle, a certain volume of pressurized air 
must be vented to the atmosphere. This volume of vented air 
equals at least the difference between the volume of the lock hop- 
per and any dead space inside it, as for instance that occupied 
by valves and rods. Therefore, it follows that the energy required 
to pressurize this particular volume of air up to the working pres- 
sure of the system is wasted. 

Inasmuch as the pressurizing air displaces the whole bulk vol- 
ume of coal, then an equal volume of air per cycle has to be 
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vented. This represents power loss for each cycle according to 
the volume and pressure of air involved. In such case, since pres- 
surizing air is furnished by an auxiliary air supply external to the 
hopper, the minimum power loss assuming polytropic compres- 
sion is given in hp-hr/ton of coal by, 


Wmin = CRT, loge r/nisotn S (3) 
where 
C constant equals 1.01 x 10~¢ 
S = (pm/M,r) — 1/1 — (prv/Pre) (4) 

The weight-ratio of coal to air, S, at varying compression 
ratios, as given by equation (4), has been explained and derived 
in Appendix 2. Fig. 6 shows the variation of S with the com- 
pression ratio r for different final bulk densities of the coal-air 
mixturé as a parameter. The isothermal cycle efficiency in equa- 
tion (3) for multistage compression has been assumed constant at 
0.90. The minimum power loss in venting the pressurized air from 
the lock hopper per ton of coal has been calculated and shown 
plotted in Fig. 7 versus the compression ratio for different final 
bulk densities of the coal air mixtures. 

Means of Reducing Power Loss. As can be seen from Fig. 7, the 
power loss in venting the pressurizing air from the lock hopper is 
quite appreciable and deserves consideration. One of the methods 
for saving as much as possible of the waste power, is to pressurize 
a second hopper before venting, in which case two lock hoppers 
are used successively. 

In order to fully understand this concept, let us consider a 
system of hoppers, as shown schematically in Fig. 8, where for 
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simplicity two lock hoppers are employed instead of one. Assum- 

ing that the compression inside the hoppers is nearly isothermal in 

case of using coal, due to the continuous intercooling effect of the 

coal as explained in Appendix 3, and instead of venting the pres- 

surized air, adjacent hoppers can be consecutively pressurized in 
series. : 

Then, suppose hopper 1 is originally at a high pressure P; after 

i | 1 discharging coal into feeder hopper, while hopper 2 is charged at 

" _ sai atmospheric pressure Po. In hopper 2, air in the voids V2 has 

CHERESEED: REND £ to be pressurized by a part of the bulk volume of air V; in hopper 


Fig. 6 1, where 











V2 = 1/p,[1 — p/p) 


When the second hopper is pressurized from the first, the final 


LOCK HOPPER COMPRESSION LOSS IN sssure j ra P. s . > same “here . 
vents GneenumeEte aan pressure in both hoppe rs P, should be the same. Therefore, 
P,V, = P«~V, + V2) and hence 





ASSUMED 0.9 FOR MULTI- STAGE 
| COMPRESSION P,/P; = V/(V; + V2) = 1/[2 — p,/p,] (5) 


| / 
| / : z 
one, the final pressure P, in the hoppers becomes, 


g 
Nsotn. 








Similarly, if n number of hoppers are pressurized from the first 


P, = P,/[n — (n — 1)p,/p,) (6) 








Fig. 9 illustrates the ratio by which P,, is reduced in proportion 
to P; for 2, 3, and 4 equal hoppers pressurized in series. It is 
noticed that a considerable reduction in pressure to be vented 
can be obtained. If the hoppers are of equal volume, the pres- 
sure of air to be vented after pressurizing two, three, or four hop- 
pers approaches 0.5, 0.34, and 0.25, respectively, as the bulk 
density of the mixture is reduced. Evidently, this results in a 
marked economy in power loss if, of course, the excess in capital 
cost could then be justified. 

Referring to Fig. 7, at 20 atm, instead of having a power loss 
of 3.0 hp-hr/ton of coal at a bulk density of 30 lb/ft’, it will be 
only 1.9 when two hoppers are pressurized in series. The power 
economy per cycle will be 1.1 hp-hr/ton and this has been shown 
«4 clearly in Fig. 10. 

oe On the other hand, even by achieving more steady and frequent 
ae | coal flow in the circuit using more than one lock hopper to reduce 
oe the power loss due to venting, sufficient justification for compensa- 
tion of the excess capital charges may be obtained in a large 


resicv fT -apacity plant than in : ll i 
Pry" 2° capacity plant than in a smaller unit. 
Lock Hopper With a Moving Wall. A positive displacement or 
moving-wall lock hopper is another attractive method for reduc- 


20 30 ing power loss by venting. Instead of feeding from a pressurized 
COMPRESSION RATIO I ay F , : 

tank, the principle of a moving wall is to develop a continuous 

Fig. 7 positive-displacement pressurizer which conveys into a flowing 
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stream of high-pressure air whatever quantity of pulverized or 
crushed coal it may receive from a separate feeder. It is possible 
also to increase its volumetric efficiency by applying positive 
aspiration (vacuum suck) during filling and excess high pressure 
as it empties. 

In practice, the only two available types of solid-fuel pumps 
working on this method are the rotary pumps designed by Yellott 
in U.S. A. [2] and by the Incandescent Heat Company in England 
{3}. As shown schematically in Fig. 11, these two pumps are 
designed for applications where the work pressure need not exceed 
10 atm. The speed of the rotary disk carrying coal to the high- 
pressure zone can be anything up to 60 rpm and coal flow rates up 
to 2.5 ton/hr, while at full load the pressure drop through the 
pump is between 3 to 5 psi. 
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The size of this pump will be determined by the rate of flow 
through the pockets in the rotary disk, rather than by the hopper 
volume required. This also determines the size of piping for 
maximum gravity flow, since flow is proportional to (pipe 
diameter)’, 

The practical problems center upon, 

1 Disk wear due to abrasion caused by entrained fine particles 
of coal during its traverse between sealed pads or covers from at- 
mospheric hopper to high-pressure air duct. 

2 Perfect sealing between high-pressure zones, while main- 
taining at the same time efficient operation of the sealing 
medium upon surfaces in contact. 

3 To avoid venting of pressurizing air from the pockets be- 
fore reaching the atmospheric hopper each cycle. This can be 
overcome by momentarily closing the: pockets after the coal dis- 
charge, by means of mechanical or hydraulic thrust, if the pres- 
sure to be regained is not exceedingly high. 





A- POWER LOSS BY VENTING AT 20 ATMOS WHEN ONE 
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8- POWER LOSS WHEN TWO HOPPERS ARE USED 
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A similar rotary pump with displaced pockets, Fig. 12, is pro- 
posed to help avoid venting of pressurized air if the excess 
mechanism will prove cheaper than venting the pressurized air. 

In comparison to the power loss by venting from a simple lock 
hopper, the probable irrecoverable power required in a lock hop- 
per with a moving wall is utilized: 

1 To introduce a given volume of coal into the high-pressure 
air stream. 

2 To compensate for any additional mechanical or hydraulic 
power exerted after the coal discharge to prevent pressurized air 
in the pockets from being vented before passing over to the 
atmospheric zone. 

3 To overcome dry friction loss of the rubbing seals. 

The minimum power loss in a lock hopper with a moving wall is 
therefore divided into: 

1 Power required to introduce a given volume of solid fuel 
into the pressurizing air stream in hp-hr/ton of coal is given by 
(14.7)(144)(2000)(r — 1)/(550)(3600)p, 

= 2.14(r—1)/p, (7) 


W, = (P, — P,)/p, = 


2 Power to displace pressurized air from the pockets of the 
same bulk volume as that of the coal in hp-hr/ton coal is given by 


We = 2.14 (r — 1)/p, (8) 
3 Power loss in dry friction in hp-hr/ton is expressed by 


W; = peripheral friction load to seal pressure 
xX speed of travel/33000 


If D, is the diameter of the cylinder in ft and yu is the dynamic co- 
efficient of dry friction of the solid on steel and assuming that the 
bearing pressure on the seal is at least equal to the maximum 
pressure of air in front of it or r atmospheres. Therefore, if A, is 
the sealing projected area in sq in. 
Friction power loss = (14.7)(2000)(144)rA,u/(33000)(60)( 2/4) 
p.D? 
or W; = 2.72rA,u/p,D? (9) 
Adding up equations (7), (8), and (9) gives the total minimum 
power loss in terms of r, and is represented in Fig. 13 for coal in a 
cylinder diameter of 5 in. sealing projected area of 5 sq in. and p 
assumed 0.3. 
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Positive-Displacement Ram Pump. This type of pump has proved 
satisfactory in operation to pressurize solid fuel at higher pres- 
sures than 10 atm. Fig. 14(a) depicts a diagrammatic arrange- 
ment of a ram pump and its cycle. It operates by drawing air- 
laden solid-fuel dust into a cylinder through a nonreturn valve by 
means of the suction stroke of a reciprocating ram, which is con- 
nected to a conventional crankchain assembly. The ram is 
usually hollow and cooled during operation. Situated at a mid- 
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point in the cylinder wall is a cam-operated, nonreturn valve 
which admits high-pressure air to the cylinder. At the start of 
the compression stroke, the tappet valve opens to admit air hav- 
ing a pressure at least equal to that of the system being supplied 
with fuel. Throughout most of the compression stroke the con- 
tents of the cylinder are discharged into the receiving system 
through a nonreturn outlet valve. 

Unfortunately, the main drawback of a ram pump is that it 
suffers from a limited low filling density of the solid fuel in order 
to avoid briquetting. Therefore, careful control of initial or 
filling density must be assured, particularly at higher pressures, 
otherwise damage of the pump is likely to occur. This necessi- 
tates the use of a sensitive quantity governor for the solid-fuel 
charge. Referring to equations (19) and (20) of Appendix 2, a 
relation between the initial and final solid-air mixtures versus the 
compression ratio is represented in Fig. 15. In practice, it has 
been found that briquetting of the solid fuel inside the cylinder 
may occur at 20 atm, if the final mixture density exceeds 30 lb/ft’, 
which corresponds to about 2.5 lb/ft* initial density, as shown in 
Consequently, a very low solid-fuel 


I ig 15 by the broken line. 


volume of expected, which is 
This, on the other 


hand, may be compromised by increasing the pump speed or num- 


output per unit pump may be 


greatly reduced by the rise in temperature. 
ber of cylinders in the cycle. Nevertheless, low filling-mixture 
density results usually in an excellent fluidization of the solid- 
fuel particles and, wherever this is sought, the ram pump may 
prove very useful. 

Behaving similar to the lock hopper with a moving wall, the 
ram pump minimum power loss can be determined directly from 
Fig. 13 by the curve of final mixture density 30 lb/ft?, which 
represents the maximum practical limit for 20 atm to develop 
Apparently, for higher pressures than 20 atm, 
to be obtained in exactly the same 


briquetting 
other final densities have 
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way Fig. 15 is calculated and then substituted in equation (8). 
In order to avoid any probable damage due to briquetting 
which may occur as a result of uncontrolled pressure, the cylinder 
of a ram pump is sometimes equipped with a clearance volume 
slightly larger than the volume of the solid fuel. In this case, in- 
stead of introducing pressurized air into the cylinder, air in the 
voids of the solid is compressed by the ram up to the required 
pressure and scavenging air assists at the end of compression to 
vent the cylinder to the atmosphere as shown in Fig. 14(b). The 
power loss in venting this system in hp-hr/ton of coal will be 


W =CRT,(loger — (1 — 1/r)]/nisotn S 10) 
where 
C is a constant equals 1.01 « 10 
S is the coal-to-air ratio that gives a final mixture density of 


30 lb/ft® at 20 atm 


Transactions of the ASME 





Solid Extrusion. It is known from work by Koppers [4] that a 
coal plug 20 in. long extruded by means of a 2-in-diameter 
plunger coaxially located in a 2.5-in-diameter tube at a rate of 120 
strokes per min and an extrusion pressure of 700 atm can provide 
a satisfactory seal against 30-atm gas pressure, and that the 
plug can be repulverized very easily in a simple air-jet micronizer. 
This establishes the practicability of the solid extrusion method 
as illustrated in Fig. 16. 
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Substantial briquetting occurs, however, when the greater part 
of air in the voids is removed by the extrusion pressure. This air 
escapes past the cylinder wall when the coal plug advances therein 
under the action of the plunger. A small part of air in the voids 
remains in the form of more or less small compressed bubbles and 
this air together with any pressurized air from the high-pressure 
chamber which may have penetrated into the plug tends to ex- 
pand during the return stroke of the plunger. This tends to press 
coal particles against the cylinder wall, thus increasing the sealing 
pressure of the plug, particularly in that constricted passage 
nearer to the high-pressure chamber. 

Sufficient compression by the plunger is, therefore, required to: 


1 Maintain high friction-bearing pressure between coal plug 
and cylinder wall in order to hold it in place against the pressurized 
air from below, particularly in the return stroke of the plunger. 

2 Achieve satisfactory sealing pressure of plug, depending on 
type of coal and its fineness. 


It is possible to reduce the friction pressure, which contributes 
to excessive wear of the cylinder, if the diameter of the plug is 
decreased to a practical minimum while the plunger stroke and 
speed are correspondingly increased. Although this procedure 
undoubtedly reduces the plunger loading and may increase the 
sealing efficiency of the plug, yet it has to be determined experi- 
mentally whether or not it significantly reduces wear. 

From test work by Hock and Roenneke [5], it is also known 


that, when air-dried granular coal is compressed in a hydraulic 
press (tests were carried out in a 16-in-bore press at 1.2-in/min 
speed), it forms a coal briquette the specific gravity of which 
increases appreciably. For example, at 200 atm, the specific 
gravity of coal briquette has reached as high as 2.5 times its 
value before compression. This increase in specific gravity de- 
pends to a large extent on the coal size and type and also on the 
limit of compression, but the exact physical explanation of the 
phenomenon is not yet known. It also has been found that only 
part of the work done is utilized in the formation of coal logs by 
briquetting, which means that much reduced optimum extrusion 
pressures may be likely obtainable only by experiment, to form 
a satisfactory sealing plug and reduce plunger loading. Some 
experimental results from the work of Hock and Roenneke are 
quoted in Table 1 and shown plotted by a broken line in Fig. 17, 
in order to show the marked increase in specific gravity of gran- 
ular coal with applied pressure. 

The minimum power loss in extrusion has been calculated as 


follows: 


1 Power to extrude the coal plug at higher plunger pressure P2 
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Table 1 


Particle size Ppoured 


--Briquetting-test— 
data 


Poriaq Pressure, Power, 


Pbriq 
lb/ft’ Pabs atm hp-hr/ton 
.475 20 .22 
44 20 AY 
. 20 20 .59 
.38 100 0 
85 110 8 

200 9 
300 .35 
300 9.0 
220 .05 
300 .76 


Pabs 
lb/ft? 
120 


U.S. mesh lb/ft 
0-16 M 61 

0-60 M 45.5 120 
0-200 M 28. 120 
0-16 M 45. 90 
0-60 M 39. 90 
0-200 M 29 90 
0-16 M 30. 97 
0-150 M 22. 97 
0-16 M 37 107 
Pulver. 13. 137 


Type of coal 
Anthracite 
Anthracite 
Anthracite 
Anthracite 
Anthracite 
Anthracite 
W. Brown 
W. Brown 
M. Brown 
Graphite 


ee Ce ee 
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against lower delivery pressure P; in the chamber in hp-hr/ton of 
coal is given by 


Pbriq 
W ea 2.14( y= to Pobrigq at r2 
Pa 


where rz and r; are the compression ratios at extrusion and cham- 
ber pressures. Equation (11) is integrated on the basis that ppriq 
is a function of P2, which is shown by the broken line in Fig. 17. 

2 Power loss to overcome friction of coal plug against the 
cylinder wall in hp-hr/ton may be determined in the same way 
equation (9) was derived and is given by 


(11) 


We = 2.72(re — 11)A,M/Poria at r2 D.? (12) 


where Ppriq at r, iS the density of the coal briquette at the ex- 
.trusion pressure. Equation (12) has been calculated for bi- 
tuminous coal which bears a friction coefficient on steel of about 
0.5. D, and A, have been selected 2.5 in. and 12.5 sq in., respec- 
tively. 

3 Power loss in compressing the air in the voids may be ob- 
tained in a similar method to that of equation (3), yielding the 
expression 


= CRT, loge (T2 — 11)/ Nisoth S re—ri) (13) 


where 
constant equals 1.01 x 1078 
average coal-to-air ratio between r2 and r; 


These power losses are shown in Fig. 17 and, as can be seen, W2 
due to friction caused by extrusion is quite appreciable in this 
case, 

Slurry Pump. Usually coal-water slurry contains from 37 to 45 
per cent water by weight, depending on size of the coal particles 
and type of coal. 

The slurry may be considered saturated with water when the 
bulk density of its mixture approaches the absolute density of 
coal, which often occurs when the water content reaches approxi- 
mately 40 per cent by weight. Oversaturation, however, takes 
place when the water content exceeds that limit and conse- 
quently the water, unless prevented, starts to drain out of the 
slurry mixture. 

However, the bulk density of the coal particles in the slurry may 
vary considerably with the water content, but the density of the 
slurry as a whole remains between that of water and coal. At 
the saturation limit, the specific volume of air in the voids is 
diminished to practically nothing, as it is replaced by an equal 
volume of water. 

In order to be able to assess the power loss in pumping the coal 
slurry, it is important to know the coal-to-water ratio by weight. 
Therefore, assuming that z is the fraction of water content in the 
slurry, then 

Weight of water = zp, 
Weight of coal = (1 — z)p, 
the ratio of coal to water in slurry will be 


Hence, 


Si = (1 — 2)/z = (1/z) - 1 (14) 


Typical experimental results by Leven [6], which have been 
quoted in Table 2 on wet coal having an absolute density of 80 


lb/ft? and an initial coal bulk density of 54.0 lb/ft’, proves the 
validity of S; in equation (14). 

In calculating the power loss to pressurize the coal slurry in this 
vase, the friction loss due to the presence of water can be ne- 
glected. Therefore, the minimum power loss in hp-hr/ton of coal 
slurry is divided into 


1 Power loss to pressurize coal from P; to P2 per ton of coal 
and denoted by W,, is given by equation (7). 

2 Power loss to pump the water in the voids to the higher pres- 
sure P: in hp-hr/ton of coal, is given by 


We = (14.7)(144)(2000)(r — 1)/(62.4)(550)(3600)S, 


= 0.0343(r — 1)/S, (15) 


Adding W; and W; yields the minimum power loss to pump the 
slurry. This is plotted against the compression ratio in Fig. 18 

It can be noticed from Fig. 18 that hydraulic pressurizing of 
coal-water slurries by lock hoppers with moving walls or ram 
pumps, reduces the power to almost that required for positive 
displacement of the solids alone, but this gain is swamped by the 
thermal effect of increased moisture content after dewatering 
and by the serious complication of pressure-drying, even if the 
residual moisture can be utilized effectively in the system. This 
can be visualized in Fig. 19, illustrating a schematic design of a 
slurry pump which can be employed, for example, in hydraulic 
elevation of coal from pit-bottom. Nevertheless, this pump may 
have a place where the solid is in any case being transported hy- 
draulically, such as recycle slag, limestone, or coke after wet de- 
dusting; and where some residual moisture is welcome as the 
case of the top feed of a fixed-bed producer. In application, this 
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Table 2 


Bulk density 
of coal-water 
mixture, lb/ft? 


Weight of 
coal, 
lb/ft? 
45.5 
49.6 
51.0 


Water content, 
per cent 


20 56.9 


30 71.0 
40 85.0 
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Weight of Weight ratio of Calculated S;, 
water, coal/water from 
lb/ft? S equation (14) 

11.4 4. 4.0 
21.4 2.32 2.33 
34.0 1.5 1.5 
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type of pump may also entail considerable design work, including 
particularly compact arrangement of water separators or filters 


inside the pressure system. 


Mechanical Considerations 

One of the deciding factors in the design of an efficient solid-fuel 
pump is the wear problem of surfaces in moving contact. Ap- 
parently, this inevitable effect occurs mainly when a rubbing seal 
is employed. In the simple lock hopper, however, there is no 
rubbing seal involved, since pressurized air has to be vented at 
the end of each cycle. Consequently, wear is only limited to that 
of the valves and seats, and can be neglected. 

In the lock hopper with a moving wall, however, the case is 
different, as the coal charge is conveyed from a zone at at- 
mospheric pressure into another sealed zone at a higher pressure. 
This incorporates the use of a rubbing seal between the surfaces of 
disk or ram and embodying cylinder. With the action of con- 
tinuous filling and discharge of the coal, there is a great tendency 
to have leakage of fine coal dust past the surfaces in moving con- 
tact, causing tremendous abrasion which is further increased by 
speeding up the cycle. 

In order to determine the practicability of a solid-fuel pump, it 
may be reasonable to develop a criterion for estimating the possi- 
ble potential wear of its rubbing surfaces in the form of 

Area of machined sealing surface in contact with coal 


Pressurized coal charge/cycle operation 


This relation, denoted as ‘“‘sealing factor,’’ can only be deter- 
mined experimentally and plotted versus the time of operation 
before which leakage or seizure occurs. 


Continuous-Action Group 


Positive Displacement Gear Pump. In this type of pump, positive 
displacement can be approximated and leakage reduced by mesh- 
ing several gears together, three to four gears being sufficient. 
This is decided according to whether or not a maximum sealing 
distance is obtained due to a central zone of intermediate pres- 
sure. Undoubtedly, this is a practical variant of the Yellott de- 
sign to eliminate venting, but in this case it is probably the wear 
of the rotary-disk seals which limits the speed and capacity rather 
than rate of filling and:discharging. 

Since wear could be very high even with relatively thick gear 
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SUCTION 


Proposed design of hydraulic slurry pump 


wheels, this pump may be a more effective way to secure high 
capacities than attempting speeds far above 60 rpm, as is the case 
with the Yellott pump. Leakage around the sides is actually re- 
duced as speed is raised in the multigear system, because it is 
more readily taken up by the low-pressure voids. For very high 
delivery pressures, two or even three units might be superim- 
posed. 

Multistage Centrifugal Compression With Air. A small rubber- 
lined multistage shrouded centrifugal compressor offers another 
possible method for pressurizing finely pulverized materials, with 
an intake density below 2.3 lb/ft* for 20-atm delivery. This, how- 
ever, is the limit over which briquetting of the solid might occur 
and cause damage of the impeller vanes. This pump has the ad- 
vantage of keeping the solids fully entrained, and their cooling 
effect at pressurizing will contribute something to efficiency. 
Nevertheless, excessive wear of the casing is liable to occur, par- 
ticularly at higher impeller speeds and solid concentration. Obvi- 
ously, this limits the practicability of this pump. 

Peristaltic or Rubber Pump. A peristaltic pump based on the 
principle of squeezing rubber ducts by a set of cams mounted on 
the same shaft, provides a promising method for pressurizing 
solid fuels in the low pressure range up to 6 or 7 atm. In this 
case, wear, leakage, and air venting are eliminated, while the 
pump offers a steady continuous flow path for the coal, even so 
it has the disadvantage of a pressure limitation with the unsup- 
ported side of the rubber duct. This is illustrated in the design 
shown in Fig. 20, where it has been established that collapsing the 
rubber at intermediate points along its length is more feasible 
than by the longitudinal roller system. 


Conclusions 

1 Power loss in pressurizing solid fuels cannot be totally 
eliminated, but can only be reduced. Since most of the power in 
pressure systems is recovered (as in the open-cycle gas turbine 
or pressure gasifier or even in the supercharged boiler), it is better 
to consider the irrecoverable power loss in comparing the various 
methods of solid-fuel pumping. Based on this principle, the 
simple lock hopper has the maximum loss due to venting. By 
utilizing a second hopper, this loss can be substantially reduced to 
nearly 60 per cent at a final bulk density of 30 lb/ft’, approaching 
that of a lock hopper with a moving wall at 20 atm. 

2 Two conflicting factors seem to govern the choice of a solid- 
fuel pressurizer in small and large-capacity plants. In a fluidizing 
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Fig. 20 Solid fuel rubber pump 


process of comparatively small capacity, capital and main- 
tenance costs are the deciding factors in design. It is likely that 
any reduction in capital cost of a lock hopper with a moving wall 
due to simplicity and compactness may be balanced by an in- 
crease in maintenance cost owing to its speedy action and wear. 
On the other hand, the power loss is the major factor in the design 
of a Jurge plant having a capacity of, say, 100 tons/hr. In this 
case, development of better pressurizers with lower power loss 
promises significant savings up to 5 per cent, or even higher, on 
fuel costs. 

3 Practical reliability for longer continuous operation of a 
satisfactory pressurizer should be determined only by experiment. 
Wear, of course, is the major trouble, particularly when a pres- 
surizer has to be operated at higher pressures and speeds, as does 
the lock hopper with a moving wall. Development of a means to 
reduce wear to a minimum in this system would probably make it 
the best answer for pressurizing solid fuel. 

1 At high pressures, ram pumps suffer from the risk of briquet- 
ting and const quent damage, unless it is required otherwise for 
low filling densities of pulverized fuel, involving higher rates of 
fluidization 

5 Itis unlikely that a device will be available for pressurizing 


solid fuel that will have lower capital cost, maintenance require- 
and power loss, together with better reliability, than the 


ments, 
simple lock hopper. If compactness due to space limitation and 
high frequency of pressurization are of major importance, the lock 


hopper with a moving wall will be the best alternative. 


APPENDIX1 


Relation Between Flow Rates of Granular Solid by Positive 
Aspiration to That by Gravity 


Treating this relation as a dimensional-analysis problem, the 
probable differential equations of flow for a mass of equal spheri- 
cal particles may be expected to contain the following physical 


quantities 


Q = F(D, d, p, P, B, ®) (16) 
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Four dimensionless combinations from these seven quantities 
are possible; d/D, 8, ® being three easily picked out. If (DzpP# 
Q) is another group, then z, y, z must satisfy the conditions of 
units of mass, length, and time pertaining to these quantities. 
This yields the following relation from equation (16) 

Q-'D%p-6P0.5 = F(d/D, B, ®) 
Hence, 


(17) 


Q = D%p°-5 Po-5/ F(d ‘D, B, ®P) 


Assuming that the bulk density remains the same in both cases 
of flow by aspiration (suction) and by gravity, then from equation 
(17) 


= Qatm V Prop! Pam - Qatm Vi + AP; Pass \ 18 


APPENDIX 2 


Relation Between Weight Ratio of Coal to Air and Final 
Bulk Density at Different Compression Ratios 
Bulk density 


total weight of the individual components in mixture 


Qasp 


total space occupied 

Wet granular coal consists of coal, water, and air. It is as- 
sumed that compressed air is admitted to a hopper filled with 
granulated coal. Due to the accompanying rise in temperature of 
the air due to compression, the coal is dried and hence the moisture 
is eliminated. 

Assuming that p,, and p, are the initial and final bulk densities 
of the mixture, then 


Py = (M, + M,)/(1 + (M./p,,.)] (19) 


With the presence of coal acting as an intercooler, compression 
of air becomes very nearly isothermal and thus: 


pp = (M, + M,)/(1/r) 4 (20) 


P (M /Pa.) 
Since, 

J = M./M, 
then, 

1)/(1/M vr) + (S/p,. 
from which 


S = (p,/M,r) — 1/1 — (p,/p,.) 


This ratio is shown in Fig. 6. 


APPENDIX 3 
Condition of Isothermal Compression of Air-Coal Mixtures 


The presence of coal in the compressed mixture acts as an in- 
tercooler, so that even with adiabatic compression, the rise in 
temperature of the mixture is inappreciable. 

This is illustrated in the following example, assuming the worst 
limiting case of briquetting, where the final mixture density is 30 
lb/ft® and its corresponding coal-air ratio is 31. 

If the compression of air commences adiabatically from room 
temperature at 520 R then, 


atr = 20, 
Therefore, 


and the heat input 


AQ = 0.24 708 = 170 Btu/Ib of air 
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When the coal particles are present in the compressed stream of 
air, and taking a specific heat of coal of 0.3 Btu/Ib deg F, then the 
mean temperature differential A7’,, of the mixture is obtained 
from 

AQ = 0.24AT7,, + 31 X 0.3 x AT,, 


from which, 
AT,, ~ 18 F 


However, A7',, becomes smaller at low values of r and, there- 
fore, even at the briquetting limit which bears the least coal-to- 
air ratio, conditions are very nearly isothermal. 
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DISCUSSION 
Andrew W. Jenike® 


The problem of pressurization of powdered solids and the 
related problem of counterflow of solids and gas have not received 
the attention they deserve in view of their frequent occurrence 
not only in combustion but also in many other processes in the 
chemical and petroleum industries. Dr. Aref is, therefore, to be 
congratulated on his timely and comprehensive presentation. 

I would like to make two specific comments. The first com- 
ment refers to the parameters of equations (1) and (2). Both 
these equations contain 6, the tangent of the angle of repose of 
the solid. $8 is assumed constant for a given ratio d/D of particle 
size to orifice diameter. It does not appear to me that the angle 
of repose is representative of the properties of the material during 
flow through an orifice. Instead of the angle of repose the angle 
of internal friction or, better yet, the effective angle of friction 
should be used. 

The error introduced by using the angle of repose instead of the 
effective angle of friction may be significant because the latter, 
and for that matter the angle of internal friction, is not constant 
for a given size of material but varies with pressure, especially for 
low values of pressure such as likely prevail in an orifice. There is 
strong circumstantial evidence that these angles vary all the way 
to 90 deg at zero pressure. This would make coefficient 6 tend 
toward infinity. If the tangent of the effective angle of fric- 
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tion is used instead of the tangent of the angle of repose, it is 
possible that equation (1) will retain its validity below the 
limit d/D = 0.05 which is imposed in this paper. 

The second comment refers to the use of a vertical, slightly 
convergent spout as a pressurizer. Such a spout can be designed 
to give no air loss at all. This will occur when the downward 
velocity of the solid is greater than the upward velocity of gas 
through the solid. In practice, this may not be attainable 
because of the excessive height that would be required to pro- 
duce a sufficiently low pressure gradient in the gas. However, in 
many cases a spout will provide the cheapest practical pressurizer. 


John J. Yellott® 


The author has given a very interesting theoretical analysis 
of the problems involved in “pumping” granular materials into 
high pressure areas. This was one of the most difficult problems 
which had to be solved during the research and development 
program of the Locomotive Development Committee, Bitumi- 
nous Coal Research, Inc. 

Reference has been made by the author to the LDC coal pump, 
but he had apparently not had an opportunity to read the story 
of the development of this pump, which appeared in unpublished 
ASME paper No. 54—A-201, presented at the 1954 Annual Meet- 
ing of the Society by the writer and his colleagues. In this paper, 
a detailed description is given of the steps which were taken to 
develop a successful means of feeding crushed or pulverized coal 
into a 100-psi compressed air stream. In this program, each of 
the five processes listed by the author under the heading “Cyclic 
Action Group”’ was given an actual test under full scale condi- 
tions. The difficulties which were encountered are far too nu- 
merous to be listed, but, in brief, it can be pointed out that the 
discharge from lock-hoppers is very difficult to control with the 
accuracy which is needed in operating a direct-fired gas turbine. 
Screw type pressurizers, which produced what the author calls 
“solid extrusion,” proved to be feasible for relatively low pres- 
sures, and, with very long screws, crushed coal could be delivered 
against pressures as high as 70 psig Tapered screws, similar to 
those used in plastics extrusion, were able to deliver pulverized 
coal against pressures up to 100 psi, but the coal was so compacted 
in the process that it entered the feed line as cylinders, rather 
than as powder. 

Piston pumps were used with some success, but after one had 
been designed which operated for more than 100 hours, against 
140 psi, it was concluded that this line of attack, even if successful, 
would not do the job for the coal-burning gas turbine. 

Rotary pumps, which were originally conceived as very close- 
clearance star feeders, were next approached, and an encouraging 
degree of success was attained with the first model. The dif- 
ficulties which were actually encountered were somewhat dif- 
ferent from those predicted by the author, because it is the teeth 
of the pump which wear, rather than the disk, and the sealing 
between the high pressure discharge region and the low pressure 
inlet region was accomplished by a relatively simple face seal. 

It was originally thought that the coal pump would run at 
constant speed, delivering to the conveying air line pulverized 
or crushed coal which would be supplied by a controllable feeder. 
However, it proved to be more difficult to develop a satisfactory 
feeder than to use a coal pump, with a variable speed drive, to 
perform both functions. Using a pump with a 5-in. wide rotor, 
the rate of pulverized coal delivery varied almost linearly from 
0 to 4000 lb per hour as the pump speed was increased from 0 to 
52 rpm. 

The major problems which the author is likely to encounter 
with the ingenious design shown in Fig. 12 are connected with the 
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annoying manner in which finely pulverized coal insinuates it- 
self into narrow spaces, closing up clearances which should be 
more than adequate, and jamming mechanisms, such as the cams 
which the author proposes. 

For extremely large quantities of fuel, such as a coal-burning 
steam turbine-gas turbine plant would require, it is quite likely 
that the best solution would be found in the combination of a 
a coal pump to do the actual feeding into 
a conveying line. This would probably be satisfactory for 
a system which operated at constant pressure, although the LDC 
coal pump proved to be a better solution with the variable pres- 


lock-hopper with 


sure gas turbine. 


Author’s Closure 


The author wishes to thank Mr. Jenike and Mr. Yellott for 
their valuable comments. Mr. Jenike brought up the specific 
question that the angle of internal friction of the solid should be 
applied in equation (1) instead of the angle of repose, to which I 
Since it has been very difficult to determine the angle 
time the 


agree. 
of internal friction 
author was studying the characteristics of fuel pressurizers as 


with sufficient accuracy at the 
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a prime objective, it was believed that the angle of repose might 
give an idea of the inter-reiations between the common factors 
involved in the problem of solid flow through an orifice. How- 
ever, it would be interesting to find out the influence of high 
pressure predominating in the pressurized system on this angle 
of internal friction experimentally, and hence on the solid flow, 
Mr. Jenike’s second comment on spout type pressurizer would 
only be feasible for very small pressures, but would be impracti- 
cable for the type of pressurizers studied in this work. 

Mr. Yellott pointed out some very interesting conclusions of 
his long experience with coal pumps. With 
problems with lock hoppers, I have had the same difficulties 
when I was checking a lock hopper feeding a pressurized coal- 
fired cyclone combustor and it was almost impossible to control 
Using a coal pump, with a variable speed 


reference to his 


the feed so accurately. 
drive, as Mr. Yellott indicated, doing both pumping and feeding 
simultaneously is an interesting approach, especially that the 
feed rate could be easily calibrated by the change of speed at 
higher pressures. The author is aware of the possible difficul- 
ties encountered with the proposed design of pump, Fig. 12, and 
the only practical way to do it is to develop first efficient longi- 
tudinal seals with sufficient cleaning air bleed. 
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ists everywhere. Over 400 reviews, are 
published monthly, highlighted by an 
article which surveys and evaluates progress 
in a specialized area. 


Published monthly. Annual subscription price: $25 
to nonmembrs, $10 to ASME members. 
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BASIE ENGINEERING 


Series D of the ASME Transactions 


Provides authoritative information on de- 
velopments in the following areas of the 
mechanical engineering field: Lubrication; 
fluid mechanics; water hammer; cavitation; 
fundamental researches in hydrodynamics; 
instruments and automatic controls; new 
metals, metal-forming processes and testing 
procedures, and better engineering ap- 
plication of metals. Research reports 
on automatic regulation theory, effect of 
temperature on the properties of metals, 
mechanical pressure elements, fluid meters, 
prevention of fracture in metals, and plastic 
flow of metals will also be published in this 
quarterly. 


Published in March, June, September, and December. 
Annual subscription price: $10* to nonmembers, $5* 
to ASME members. 











APPLIED MECHANICS 


Series E of the ASME Transactions 


Here will be found fundamental analytical 
and experimental studies in the fields of 
elasticity, plasicity, dynamics, vibration, 
impact, flow and fracture. structures, stress 
analysis, mechanical properties of materials, 
hydraulics, aerodynamics, internal flow, 
turbulence, boundary layers, aeroelasticity, 
thermodynamics, heat transfer, lubrication 
and weer, and computing devices and 
methods. 

Where subjects overlap the fields covered 
by the other publications of this Series, the 
Journal of Applied Mechanics will deal 
with the basic laws and the more advanced 
methods of studying problems. 


Published in March, June, September, and D ber. 
Annual subscription price: $10* to nonmembers, $5* 
to ASME members. 
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APPLIED MATHEMATICS 
AND MECHANICS 
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This bimonthly contains the important 
findings of Russian scientists and engineers 
working in the fields of applied mathematics, 
fluid and solid mechanics, aeronautical 
sciencies, astronautics, rocketry and space 
flight. It is a cover-to-cover translation of 
the Journal entitled ‘‘Prikladnaya Mate- 
matika i Mekhanika’’ which has been for 
over two decades one of the leading Soviet 
periodicals in the fields. 

This translation series began with Volume 
22 (1958) and is being published bimonthly 
or corresponding to the number of original 
Russian issues per annum. All subscriptions 
are on a calendar year basis. 


Annual subscription price: $35 to nonmembers, $28 
to ASME members. 





